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NOTE OF TRANSMITTAL 


The SEASAT Economic Assessment was performed for the 
Special Programs Division, Office of Applications, National 
Aeronautics and Space Administration, under contract NASW~2558. 

The work described in this report began in February 1974 and 
was completed in August 1975. 

The economic studies were performed by a team con- 
sisting of Battelle Memorial Institute; the Canada Centre for 
Remote Sensing; ECON , Inc.; the Jet Propulsion Laboratory; and 
Ocean Data Systems, Inc. ECON, Inc. was responsible for the 
planning and management of the economic studies and for the de- 
velopment of the models used in the generalization of the re- 
sults. 

The SATIL 2 program has been developed to assist with 
the programmatic evaluation of alternative approaches to estab- 
lishing and maintaining a specified mix of operational sensors 
on spacecraft in an operational SEASAT system. The program en- 
ables the assessment of the effects of operational requirements 
and reliability (spacecraft buses, sensors, and transportation 
systems) on the time phased costs of alternative approaches. 

The SATIL 2 program is specifically designed to allow for the 
explicit consideration of reliability and cost uncertainties. 

In order to perform this evaluation, the launch systems and 
spacecraft (buses and sensors) are considered in detail from 
the points of view of the major sub-systems reliability and 
costs. All costs are treated as uncertainty variables where 
ranges of possible values are considered as well as subjective 
estimates pertaining to the form of the uncertainty (the proba- 
bility distribution) within the range. As a result, the program 
develops a set of probability distributions associated with costs 
and events (i.e., number of launch attempts, etc .) as functions 
of time and the probability distribution of the present value 
of total recurring costs. All of the results are based upon an 
optimization of the mix of sensors, spacecraft buses, and trans- 
portation systems and launch schedules in order to meet the sen- 
sor demand function . 

Although the SATIL 2 program was developed to assist 
in the evaluation of the costs of operational SEASAT system 
alternatives, the program is general and can be used for the 
evaluation of the costs of any proposed geocentric satellite 
system. : 

The SEASAT Users Working Group (now Ocean Dynamics 
Subcommittee) chaired by Dr. John Apel of the National Oceano- 
graphic and Atmospheric Administration, served as a valuable 
source of information and a forum for the review of these studies . 
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Mr. S.W. McCandless, the SEASAT Program Manage 
activities of the many organizations that part 
studies into the effective team that obtained 
scribed in this report. 

The SATIL 2 mathematical model, its 
computational procedures were developed by Joe 
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1. OVERVIEW OF THE ASSESSMENT 


This report/ consisting of ten volumes/ repre- 
sents the results of the SEASAT Economic Assessment/ as 
completed through August 31/ 1975. The individual volumes 


Summary and Conclusions 

The SEASAT System Description and 

Performance 

Offshore Oil and Natural Gas Industry - 
Case Study and Generalization 
Ocean Mining - Case Study and Generali- 
zation 

Coastal Zones - Case Study and Generali- 
zation 

Arctic Operations - Case Study and 
Generalization 

Marine Transportation - Case Study and 
Generalization 

Ocean Pishing - Case Study and Generali- 
zation 

Ports and Harbors - Case Study and Gen- 
eralization 

A Program for the Evaluation of Opera- 
tional SEASAT System Costs. 


s report 

are 

Volume 

I 

Volume 

II 

Volume 

III 

Volume 

IV 

Volume 

V 

Volume 

VI 

Volume 

VII 


Volume' VlII 
volume IX 
Volume X 


Each volume is self-contained and fully documents 
the results in the study area corresponding to the title. 
Table 1.1 describes the content of each volume to aid readers 
in the selection of material that is of specific interest. 

The SEASAT Economic Assessment began during Fis- 
cal Year 1975. The objectives of the preliminary economic 
assiessment/ conducted during Fiscal Year 1975, were to identi 
fy the uses and users of the data that ; could be produced by 
an operational SEASAT system and to provide preliminary esti- 
mates of the benefits produced by the applications of this 



Table l.l; 


CoiiLent and Orgaulzation of the Pinal Report 


Volume No. 


Ti tie 



Summary and Conclusions 


A summary of benefits and costs, and a statement of the 
major findings of the assessment. 


If The SEASAT System 

Description and Per- 
formance 


A discussion of user requirements, and the system concepts 
to satisfy these requirements are presented along with a 
preliminary analysis of the costs of tlicse systems. A 
description of tlie plan for the SEASAT data utility studies 
and a discussion of the preliminary results of the simula- 
tion experiments conducted with the objective of quantifying 
the effects of SEASAT data on numerical forecasting. 


in 


Offshore Oil and 
Natural Gas industry- 
case Study and Gener- 
alization 


The results of case studies which investigate the effects of 
forecast accuracy on offshore operations in the North Sea, 
tlie Celtic Sea, and the Gulf of Mexico are reported. A 
methodology for generalizing the results to otlier geographic 
regions of offshore oil and natural gas exploration and de- 
velopment is described along with an estimate of the world- 
wide benefits. 


IV 


Ocean Mining 
Study and Ge 
iza tion 


n 


- Case 
e r a 1 - 


Tlie results of a study of the weather sensitive features of 
the near shore and deep water ocean mining industries are 
described. Problems with tlie evaluation of economic benefits 
for the deep water ocean mining industry are attributed to 
the relative immaturity and highly proprietary nature of the 
industry. 






Table 1,1: Content and Organization of tlie Pinal Report 

(ootitiiiKed ) 


Volume No. 

Title 

Content: 

V 

Coastal Zones - Case 
Stiiily and General- 
ization 

The study and generalization deal witli tlie economic losses 
sustained in the U-S. coastal zones for the purpose of 
quantitatively establishing economic benefits as a conse- 
quence of improving the predictive quality of destructive 
phenomena in U.S. coastal zones. Improved prediction of 

hurricane landfall and improved experimental knowledge of 
hurricane seeding are discussed. 

VI 

Arctic Operations - case 
Study and General t Za t ion 

The hypothetical development and transportation of Arctic 
oil and other tesovrees by ice breaking super tanker to 
tlie continental East Coast a.re discussed. SEASAT data will 
contribute to a more effective transportation operation 
through the Arctic ice by reducing transportation costs as 
a consequence of reduced transit time pec voyage. 

VII 

Marine 'J’ranspor ta tion- 
Case Study and General- 
ization 

A discussion of the case studios of the potential use of 
SEASAT occa:i condition data in the improved routing of dry 
cargo ships and tankers. Resulting forecasts could be 
useful In routing ships around storms, thereby reducing 
adverse weather damage, time loss, related operations costs, 
and occasional catastrophic lo^.ses. 

VIII 

Ocean Pisliing - Case 
Study and Generaliz- 
ation 

The potential application of SEASAT data with regard to 
ocean fisheries is discussed in this case study. Tracking 
fish populations. Indirect assistance in forecasting expected 
populations and assistance to fishing fleets in avoiding 
costs incurred due to adverse weather through improved ocean 
conditions forecasts were investigated. 

IX 

Ports and Harbors - Case 
Study and Generalization 

The case study and generalization quantify benefits made 
possible through improved wcathor forecasting resulting 
from the integration of SEASAT data into local weather 
forecasts. The major source of avoidable economic losses 
from inadequate weather forecasting data was shown to be 
dependent on local precipitation forecasting. 

X 

A Program for the Evalu- 
ation of Operational 
SEASAT System Costs 

A discussion of tlic SA'i'lL. 2 ('rogr:iin which was developed to 
assist in the evaluation of the costs of operational SEASAT 
system al.ter natives . SATIL 2 enables the assessment of the 
effects of operational requirements, reliability, and time- 
phased costs of alternative appcoachcs. 
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dataV* The preliminary economic assessment identified large 
potential benefits from the use of SEAS AT~produced data in the 
areas of Arctic operations, marine transportation, and offshore 
oil and natural gas exploration and development. 

During Fiscal Year 1976, the effort was directed to- 
ward the confirmation of the benefit estimates in the three 
previously identified major areas of use of SEASAT data, as 
well as the estimation of benefits in additional application 
areas. The confirmation of the benefit estimates in the three 
major areas of application was accomplished by increasing both 
the extent of user involvement and the depth of each of the 
studies. Upon completion of this process of estimation, we have 
concluded that substantial, firm benefits from the use of oper- 
ational SEASAT data can be obtained in areas that are extensions 
of current operations such as marine transportation and offshore 
oil and natural: gas exploration and development. Very large 
potential benefits from the use of SEASAT data are possible in 
an area of operations that is now in the planning or conceptual 
stage> namely the transportation of oil, natural gas, and other 
resources by surface ship in the Arctic regions. In this case, 
the benefits are dependent upon the rate of development of the 
resources that are believed to be in the Arctic regions, and 
also dependent upon the choice of surface transportation over 
pipelines as the means of moving these resources to the lower 

* _ 

SEASAT Economic Assessment, ECON, Inc., October 1974. 
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latitudes. Our studies have also identified that large 
potential benefits may be possible from the use of SEASAT 
data in support of ocean fishing operations. However, in 
this case, the size of the sustainable yield of the ocean 
remains an unanswered question; thus, a conservative view- 
point concerning the size of the benefit should be adopted 
until the process of biological replenishment is more 
completely understood. 

With the completion of this second year of the 
SEASAT Economic Assessment, we conclude that the cumulative 
gross benefits that may be obtained through the use of data 


from an 

operational SEASAT 

system , 

to 

provide improved ocean 

condition and weather forecasts is 

in 

the range of $859 

mi 1 lion 

to $2,709 million 

($1975 at 

a 

10 percent discount 


rate) from civilian activities . These are gross benefits 
that are attributable exclusively to the use of SEASAT data 
products and do not include potential benefits from other 
possible sources of weather and ocean forecasting that may 
occur in the same period of time. The economic benefits 
to U.S. military activities from an operational SEASAT sys- 
tem are not included in these estimates. A separate study 
of U.S. Navy applications has been conducted under the 
sponsorship of the Navy Environmental Remote Sensing Coor- 
dinating and Advisory Committee. The purpose of this Navy 
study was to determine the stringency of satellite oceano- 
graphic measurements necessary to achieve improvements in 
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military mission effectiveness in areas where benefits are 
known to exist.* It is currently planned that the Navy 
will use SEASAT-A data to quantify benefits in military 
applications areas. A one-time military benefit of approx- 
imately $30 million will be obtained by SEASAT-A, by p.TO- 
viding a measurement capability in support of the Depart- 
ment of Defense Mapping, Charting and Geodesy Program. 

Preliminary estimates have been made of the costs 
of an operational SEASAT program that would be capable of 
producing the data needed to obtain these benefits. The 
hypothetical operational program used to model the costs of 
an operational SEASAT system includes SEASAT-A, followed by 
a number of developmental and operational demonstration 
flights, with full operational capability commencing in 
1985. The cost of the operational SEASAT system through 
2000 is estimated to be about $753 million ($1975, 0 per- 
cent discount rate) which is the equivalent of $272 million 
($1975) at a 10 percent discount rate. It should be noted 
that this cost does not include the costs of the program's 
unique ground data handling equipment needed to process, 
disseminate or utilize the information produced from SEASAT 
data. Figures 1.1 and 1.2 illustrate the net cumulative 
SEASAT exclusive benefit stream (benefits less costs) as a 

* 

"Specifications of Stringency of Satellite Oceano- 
graphic Measurements for Improvement of Navy Mission 
Effectiveness." (Draft Report.) Navy Remote Sensing 
Coordinating and Advisory Committee, May 1975. 
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Figure 1.1 SEASAT Program Net Benefits, 1975-2000 
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function of the discount rate. 


This volume describes a mathematical program called 
SATIL 2 which has been developed to aid in the programmatic 
evaluation of alternative approaches to establishing and 
maintaining a specified mix of operational sensors on space- 
craft in an operational SEASAT system. 



10 


2. INTRODUCTION 

The evaluation and comparison of alternatives which 
are based upon technologies in various stages of research, 
development, design, and use should take into account not only 
expected cash flows and present values, but should also take 
into account the associated levels of risk. Risk may be 
characterized as the probability that pertinent evaluation 
measures exceed specified values. Major contributors to risk 
are system unreliability, performance uncertainties, and cost 

uncertainties. Because of the normally complex interactions 

; “1 

between performance, reliability, and cost, it is desirable 
to consider these factors simultaneously when establishing 
the probability distributions of measures such as annual cash 
flows and present valuer of cash flow. 

This report is concerned with the documentation of 
the SATIL 2 program which computes the probability distribu- 
tions of events (i.e., number of launch attempts, number of 
spacecraft purchased, etc.), annual recurring cost, and 
present value of recurring cost. This is accomplished for 
the specific mission of placing a mix of sensors in orbit in 
an optimal fashion in order to satisfy a specified sensor 
demand function. The optimization is based upon selecting 
the mix of spacecraft buses (the platform consisting of the 
basic housekeeping and support functions upon which the 
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sensors are mounted), sensors and launch schedule which mini 
mizes the present value of satisfying the specified sensor 
demand (as a function of time) requirements. The SATIL 2 
program models the interactive effects of reliability, 
performance, and cost. 

This report consists of a brief overview of the 
methodology for evaluating alternatives which explicitly 
takes into account reliability and cost uncertainties and 
develops quantitative risk measures. The methodology makes 
possible the comparison of alternatives ’ with the explicit 
consideration of risk. The methodology also allows an 
economic value to be established for changes in requirements 
and/or capabilities. For example, the economic value of 
subsystem reliability improvements can be assessed and can 
be of considerable importance when making reliability 
allocations . ' 

SATIL 2 is intended to work with a data set which 
describes the operational SEASAT system requirements and 
alternatives. Because of time constraints, it was not 
possible to develop' a SEASAT data set during this phase of 
the work. The development of a SEASAT data set and the use 
of SATIL 2 to evaluate the costs oi operational SEASAT 
alternatives is planned for a subsequent effort during 
FY-1976. However, to demonst rate the methodology and the 
capability of the SATIL 2 program, the economic consequences 
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of unreliability are considered for a hypothetical space 
mission requiring the emplacement and maintenance {for an 
extended period of time) of a system of earth observation 
sensors in geocentric orbits. The demand for sensors may be 
satisfied by utilizing different spacecraft types, each with 

i ^ 

a different mix of on-board sensors and each requiring a 
different transportation system capability. The objective 
is to select that mix of spacecraft and launch schedule which 
provides at least the desired sensor demand and which mini- 
mizes the present value of the costs. Reliability plays an 
important role in this analysis since the number of sensors 
which must be placed in orbit each year depends upon those 
which were placed in orbit in previous years and are no 
longer operational. The hypothetical example illustrates an 
assessment of the effects of operational requirements, per- 
formance capability, and reliability (spacecraft support 
subsystems, sensors, and all pertinent operational aspects 
of the transportation systems) on the time phased costs of 
alternative approaches. It also illustrates how the impact 
of technology improvements, reliability variation, and their 
combined effects can be evaluated by considering the present 
value probability distributions. 

The value of reliability improvement and/or tech- 
nology improvements can be established as the change in the 
present value of recurring costs with and without the 
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improvements. The value of risk reduction through both 
reliability improvements and operational capability is also 
discussed. 

2 . 1 Uncertainty r Risk, and Decision Making ^ ^ 

A typical problem confronting decision makers is to 
determine the cost, present value, etc., of alternative ways 
of satisfying a mission requirement and, thence, to select 
that alternative which minimizes a performance measure such 
as total cost, present value, etc. Since, in the aerospace 
industry, the decision is usually concerned with selecting an 
alternative from among a Set of alternatives which utilize 
different technologies, are in various stages of planning, 
research, design and development, different levels of uncer- 
tainty, and hence, risk exists. The comparison and selec- 
tion of the various alternatives should consider the different 
levels of risk. The addition of the risk dimension is 
illustrated in Figure 2.1 and implies the consideration of 
reliability and cost and performance uncertainties. 

The economic analysis of space programs involves 
cash flow patterns that occur over several (n) years; thus, 
it is desirable to present the results of economic analyses 

I , 

(which cjombine the interactive effects of performance, cost, 
and mission requirements) in terms of the present value of 
costs. The pre sen t value , which explicitly takes into 




/> 
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account the magnitude and the timing of the cash flow patterns 
is defined as the summation of future annual costs discounted 
to the present and is given by 


[l+(r/100)]^ 

where PV=present value of cost, C^=cost in the i^^ year, and 
r=discount rate (%) or cost of capital. 

The costs entering into the above equation, however, 
are not deterministic quantities. Variations in the yearly 

i 

costs, occur because of the uncertainties in predicting 

future item (or per event) costs and due to the uncertainty 
both in the number of events necessary %o perform the desired 
program and the time of occurrence of these events. Thus, 
the present value of costs must also be characterized by a 
probability distribution. The probability of present value 
exceeding a specified level is the area under the probability 
density curve for all values greater than the specified 
level and is henceforth referred to as a risk profile. 

Typical risk profiles of present value are shown in Figure 2 . 2 
where the probability or chance, p, of exceeding the various 
levels of present value, PV, is indicated. In general, the 
steeper the curve , the lower the risk (or variability). When 
comparing alternatives, it is important to compare the 
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In the certainty situation/ it is generally desir- 
able to select the alternative which yields the minimum 

present value of costs when all alternatives are evaluated 

* 

on an equal capability basis. When present values are equal, 
the choice is immaterial. The selection process becomes more 
difficult when uncertainties are considered; tradeoffs must be 
made between alternatives possessing different expected pre- 
sent values and associated levels of risk. When the risk 
dimension is available, Alternatives A and B are found to be 
quite different. The risk associated with Alternative A is 

greater than that of Alternative B. Thus, a conservative 

* * 

decision maker would normally select Alternative B, provided 
that other unquantified pressures to select Alternative A do 
not exist. 

Generally, one can identify many alternatives that 
must be compared for selection of the best one. The problems 
of comparison are eased somewhat by the fact that the proba- 
bility distributions of the present value of costs are usually 
very nearly normal. Thus, the distributions can be fully 
characterized by their mean, m, and standard deviation, a, and 


•k 

Equal capability is defined as providing the same level 
of s e r vi ce . 

■k it 

Conservative in the sense that risk aversion preferences 
are evident . 
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each alternative can be represented by a point on the m-a 
plane. An example is illustrated in Figure 2.3. Here, 
Alternatives 1 and 2 have the same level of risk {i.e., 
cr^^Q^) / but the expected PV of the cost of Alternative 2 is 
greater than that of Alternative 1. Therefore, Alternative 1 
is preferable to Alternative 2. In a similar manner, it can 
be argued that Alternative 3 is preferable to Alternative 
Also in a similar manner. Alternative 3 is preferable to 

Alternative 1, since both have the same expected PV , but 

* 

Alternative 1 is riskier. This process can be continued 



Figure 2.3 General Problem of Decision Making 
Under Uncertainty 


This assumes a rational risk averse decision maker. 
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with all alternatives being considered. In the limit, it can 
be seen that a frontier of "best” alternatives can be estab- 
lished. Each of the points or alternatives represented by the 
frontier are different in the respect that the risk and ex- 
pected PV are different. The class of best alternatives has 
thus been obtained and the "best" alternative can then be 
selected based on the decision maker's risk judgment. That 
is, the decision maker must decide what the tradeoff is be- 
tween a reduction in expected PV of cost and an accompanying 
increase in risk. 

The risk judgment may be purely intuitive or it may 
be based on a quantified utility function [5 , 6 3 ■ In each case, 
the tradeoff between expected rewards and risk is made explic- 
itly clear to the decision maker. 

It is clear from Figure 2.3 that the economic value 
of pursuing Alternative 1 relative to Alternative 2 is m^-m^ 
and Alternative 3 relative to Alternative 4 is m^-m^. Thus, 
one could afford to spend up to an additional amount m^-m^ on 
non-recurring costs in order to achieve Alternative 1 and 
still be better off than with Alternative 2. Considering 
Alternatives 1 and 3 with equal expected present value of 
costs, the value of Alternative 3 relative to Alternative 1 
is clearly associated with risk reduction and can be assessed 
either intuitively or by developing utility functions. Again, 


1 


I 








» 


t 
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one can afford to spend an additional amount of up to this 
established value on non-recurring costs to achieve Alter- 
native 3 . 

The following paragraphs are concerned with the 
methodology incorporated into the SATIL 2 program for estab- 
lishing the annual cost probability distributions in terms of 
reliability and cost uncertainties and establishing the PV 
probability distribution. 

2 . 2 The SATIL 2 Methodology (Overview) 

Because of the complex interactions of performance, 
reliability and cost, the methodology [7,8,9] for evaluating 
the annual cost and present value probability distributions 
is based upon Monte Carlo simulation modeling techniques [10]. 
Specifically, a simulation model has been developed to assist 
with the programmatic evaluation of alternative approaches to 
establishing and maintaining a specified desired mix of 
operational sensors on spacecraft in geocentric orbits. The 
program enables the assessment of the effects of operational 
requirements and reliability (spaceciaft support subsystems, 
sensors, and transportation systems) on the time phased costs 
of alternative approaches to satisfying mission requirements. 
The program is specifically designed to allow for the explicit 
consideration of reliability and . cos t uncertainties . In order 
to perform this evaluation, the launch systems and spacecraft 


(support systems and sensors) are considered in detail from 
the points of view of reliability and cost. All costs are 
treated as uncertainty variables where ranges of possible 
values are considered as well as subjective estimates pertain 
ing to the form of the uncertainty (the probability distribu- 
tion) within the range. The input to the program consists 
primarily of a set of numbers which describes the demand for 
various operational sensors in orbit as a function of time, 
the mix of sensors available per spacecraft type, the trans- 

I . • • . ^ 

portation system to be used for each spacecraft type as a 
function of time, spacecraft subsystem, sensor and transpor- 
tation syst'em reliability, spacecraft and tr aQ.spp rtat ion 
system costs including explicit quantitative unceftainty 
assessments and cost learning rates. The output from the 
simulation program consists of a set of probability distri- 
butions associated with costs and events (i.e. , number of 
launch attempts, etc.) as functions of time and the proba- 
bility distribution of the present value of total recurring 
cost. 

The reliability, uncertainty and risk assessment 
capability embodied in the simulation model allows for; 


The group of support systems is frequently referred 
to as the spacecraft bus . 
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• the specification of the mix of operational 
sensors required in geocentric orbits as a 
function of time; 

• the consideration of multiple spacecraft which 
are defined in terms of the reliability of 
the major support subsystems, the mix of on- 
board sensors and their reliability and 
spacecraft cost; 

• the consideration of spacecraft subsystem and 
sensor failure models which allow for both 
random and wearout failures; 

• the specification and consideration of multiple 
transportation systems which may consist of 
current or new expendables or the space shuttle. 
The transportation systems may also include (as 
necessary) orbit- to-orbit shuttles or propulsion 
modules (for example, Agena, Centuar, Space Tug, 
etc.). The propulsion modules may be expendable 
or reusable and may be used for placing space- 
craft in orbit and retrieving spacecraft which 
fail and require replacement. The specification 
of the transportation systems include cost and 
reliability assessments. Reliability is 
considered at the major subsystem level; 
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• the specification of transportation systems to 
be utilized for placing different spacecraft 
into orbit as a function of time. Changing 
the specification of transportation system- 
spacecraft assignment as a function of time 
allows performance capability (such as 
allowable mission modes and reliability) 

and cost variations to be considered; 

• the explicit consideration of multiple time 
periods, thus allowing for annual costs 

to be established; 

• the consideration of cost learning curves 
an d 

• all costs to be treated as uncertainty 
variables . 

The simulation model, taking into account the 
required number of sensors as a function of time, number 
of operational sensors in orbit (as determined from space- 
craft subsystem and sensor reliability characteristics) and 
Spacecraft and launch costs, determines a near-optimal mix 
of spacecraft launches as a function of time. Since the 
simulation is based upon Monte Carlo te chnique s , it is 
possible to: establish the probability distributions of 
pertinent performance measures which allows alternatives 
to be compared by considering both expected values of 


performance measures and the chance of variation (i.e., the 
risk) of the values of the measures. Specifically, the 
simulation model establishes the probability distributions of 

• pertinent quantities by year (for example, 
number of launch attempts, number of space- 
craft required, number of propulsion 
modules required/ number of propulsion 
module ref urbishments , etc.); 

• launch, spacecraft, and total costs by space- 
craft type and by year and 

• present value of recurring costs. 

The basic structure of the simulation model is 

illustrated in Figure 2.4 and is typical of that required to 

* 

evaluate service-type of missions. The overall simulation 

i , ‘ 

model consists of a number of submodels, namely a Unit Cost 
Model, Sensor/Spacecraft Reliability Model, Demand/Scheduling 
Model, Operational Analysis Model, and an Annual Cost Model. 
Each of these are described briefly below. 

Unit Cost Model : The Unit Cost Model establishes 

the cost associated with each of the cost incurring events. 


* 

Those missions requiring an establishment and main- 
tenence of service (for example, remote sensing of 
earth resources, ocean surveilance, etc.). 


Input Data: 

Sensor Demand 

S pac ec rn f t /Sen so r Mix Specification 
Transportation System Definition 
Reliability Assessment 
Uncertainty Assessment 


input Data Verification 
C initialization of Variables 


Unit Cost Mode i 
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Figure 2.4 Basic Structure of The SATTL 2 Simulation 
Model 








That is, it establishes the initial cost per launch, propul- 
sion module refurbishment, spacecraft bus purchase, etc. The 
costs are established once for each Monte Carlo simulation run 
and are basically first unit costs. Future costs take into 
account learning effects and are discussed in the Annual Cost 
Model. 

Recurring cost uncertainties arise from the diffi- 
culties of predicting the cost of producing an item before it 
has been designed. An overall economic analysis must be 
concerned with these uncertainties. The problem is how to 
quantify uncertainty. This requires that informed estimates 
be made of the ranges: of uncertainty of key cost variables 
an.d. their probability distributions within the range. The 
estimates of uncertainty might be made, for example, at the 
CER (cost estimating relationship) level or they might be 
made at the unit cost level. The uncertainty assessments 
[4,11] can be made by individuals with the assistance of an 
experienced analyst or they can be made by an experienced 
group of analysts using Delphi-type techniques . The estimates 
are very subjective in nature and quantitatively express the 
attitudes regarding the uncertainties. The estimates reflect 
past experience with s imilar efforts , problems which have 
been encountered in the past , insights into problem areas 
which might develop, etc. 



A methodology for establishing the shape of the cost 
certainty profiles is illustrated in Figure 2.5 and hasi been 
employed in many previous risk analyses. The first step is 
to establish the range of unit cost uncertainty based on 
knowledgeable persons assessing what can go right and what 
can go wrong. The range is then divided into five equal 
intervals and a relative ranking of the likelihood of the 
cost variable falling into each of the intervals is performed. 
The general shape (skewed left, skewed right, central, etc.) 
of the uncertainty profile is thus established. The next 
step is to establish relative values of the chance of falling 
into each of the intervals. For example, in the illustration, 
the chance of falling into the first interval is estimated 
to be half as likely as falling into the second interval. 

This is repeated for each interval relative to the previously 
considered interval. The last step is to solve the illus- 
trated equation for the quantitative values by substituting 
the data from the previous step. 

In order to simplify this procedure, a large number 
of typical uncertainty profiles are stored in the computer. 

The evaluator may thus simply specify the range of Uncertainty 
(minimum and maximum values) and the name of the .uncertainty 
profile which reasonably represents liis feelings. If none 
of the stored profiles is suitable, then the previously 
described procedure may be followed. 



T 


j ^ 1000 _ 2000 

I d) Establish Quantitative Values 

. 

I, . 

i-"' 

I . . ■ . ■ 

[' Figure 2 . 5 Methodology for Establishing Shape of Cost 

[ V Uncertainty Profile 
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The Unit Cost Model is based upon the above method- 
ology. The Unit Cost Model employs minimum and maximum cost 
.estimates and the associated uncertainty profiles (probability 
distributions) and random sampling techniques to determine the 
unit costs which are to be used for each of the Monte Carlo 
simulation runs. 

Sensor/Spacecraft Reliability Model ; The purpose of 
the Sensor/Spacecraft Reliability Model is to determine , based 
upon estimated failure characteristics, which of the pre- 
viously-launched sensors are still operational in the time 
frame being considered. 

A service mission typically consists of two phases, 
namely the establishment of a desired level of service and the 
maintenance of that level of service. The establishment of 
the service is concerned with successfully placing a number 
of sensors in orbit over a period of time as determined by 
service needs. If it is assumed that both the need and the 
capability of the spacecraft are known, it follows that the 
number of spacecraft required to establish the service is 
deterministic. The maintenance of the service is concerned 
with maintaining the desired number of operational sensors in 
service. As a result of les S'^ than-perf ect reliability and 
variability associated with wearout phenomena, such as 
attitude control gas depletion, failures will occur in a 
random manner with the result that the traffic associated 
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with the maintenance of the service will be known only in a 
probabilistic sense. The probabilistic nature of the traffic 
due to spacecraft failures and replacements adds an important 
degree of uncertainty and risk to the service mission. It 
is, therefore, necessary to consider the impact of reliability 
on the overall mission. 

Spacecraft failures occur as the result of four 
causes: improper design or manufacture generally due to an 

incomplete understanding of all related physical principles, 
imperfect quality control during manufacture, uncertainties 
in the environment during storage, transportation and oper- 
ation, and because of design and technology limitations 
generally associated with lifetime. These failure causes 
generally manifest themselves in three different types of 
failures: early, chance, and wearout. In this paper, early 
spacecraft failures are considered in the Operational Analysis 
Model. Their consequences are dependent upon the Space 
Transportation System and the operational modes. 

The general failure model [8] utilized in this 
analysis considers both random and wearout failures and is 
given by 


R (n) 


exp ( -Xn ) 
( 2 -^) 



exp 


(a-M^) 


20 


da 
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This expresses the combined effects of the random and wearout 

phenomena as the probability of surviving through the n^ time 

period (years) and A is the failure rate (failures/year) or 

reciprocal of the mean-time between failures (MTBF ) , is 

the mean or expected wearout life (years) and is the 

f 

standard deviation of wearout life (years). The probability 
of a failure in time period n, F(n), given success to the 
beginning of time period n is thus 


r 


(n) = R (n-1) -R (n) 


R (n-1)/ 


The failure characteristics can be modeled by esti- 
mating, based on current design and/or past experience, the 
anticipated MTBF (1/A) for the random failures and the expec- 
ted wearout life and the standard deviation of wearout life. 
Although this failure model is a very simple one, it serves 
to illustrate the way in which reliabilities enter into the 
economic analysis. More detailed reliability models could 
be developed as necessary. 

The Sensor/Spacecraf t Reliability Model considers 
the reliability of each of the sensors and each of the space- 
craft support subsystems individually keeping track of when 
each is placed into orbit and determines by sampling tech- 
niques when each fails. Sensors are considered to be inde- 
pendent of each other, whereas sensor operation is dependent 
upon successful operation of all the support subsystems. 
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Demand/Scheduling Model -. The purpose of this Model 
is to determine which spacecraft and sensor mix should be 
launched each year taking into consideration the number and 
type of operational sensors at the start of each year and 
future sensor demand and unit costs. The objective is to 
minimize the present value of recurring cost by selecting the 
optimum mix of spacecraft buses, sensors, and launch systems 
and schedule such that the demand for sensors is always satis- 
fied even though sensors placed in orbit in previous years 
may no longer be operational. The problem is complicated by 
the fact that it is assumed that several different spacecraft 
buses (group of support subsystems) are available which are 
capable of supporting different mixes of sensors. 

The optimum solution can be obtained using an 
integer programming algorithm. However, because the Demand/ 
Scheduling Model is contained within a Monte Carlo loop, a 
linear programming algorithm, with its increased computational 
speed, has been utilized to find a very near optimum solution. 
The objective function and constraints for the linear program- 
ming algorithm are as follows : 



Minimi ze 


subject to 


N 


NO, + 
k , n 


SX . , X 
3 / k 


j = l 


I 


^ ND. 


i , k 


i=n 


for each i such that and I’^k^K, where 

n = start of the time horizon (years) , 

N = number of time periods (years) to be 
considered , 


PVF . = present value factor applicable to the 

i^h time period, 

J = maximum number of spacecraft bus types, 

th 


SC. . = CO sit. {$) of the j 


spacecraft (bus 


plus sensors) in the i 


th 


LC . . = cost ($) of launching the j 

in the ith year. 


year 
. th 


spacecraft 


K = maximum number of sensor types, 

NO = number of sensors of type k operational 

at the end of the n^h year, 


SX 


. th 


j f k 


ND . 


X , k 


= sensor types available on the j 
spacecraft (4> or 1), 

= number of sensors of type k required 
to be operational in the ith year. 


The above quantities are provided to the Demand/ 

Scheduling Model which thence determines NA . the number 

of spacecraft of type j to be placed successfully in orbit 

in the i^ year. NA. . thence serves as a primary input to 

1 A D 

the Operational Analysis Model. 



Operational Analysis Model ; The function of the 
Operational Analysis Model [9] is the quantification of the 
probability distributions of the numbers of recurring cost- 
associated events required to establish and maintain the 
desired level of service (as determined by the Demand/Sched- 
uling Model/ i.e., the value of NA. .). These events may 

^ » D 

include the number of spacecraft purchased, the number of 
launches, the number of spacecraft retrievals or refurbish- 
ments, and so on. Each of these events has associated with 
it a cost that, at least in the planning phase, includes 
some uncertainty. 

Operational analysis explicitly considers the 
possibility of failures in the space transportation systems, 
the chance of these failures occurring, and the consequences 
if the failures do indeed occur. More specifically, the 
analysis considers (1) the number, sequence, and complexity 
of operations to be performed, (2) the recovery modes, that 
is, given that a failure has occurred, the possible resulting 
sequences of events, and (3) the probability of successfully 
performing each of the required operations , both in the 
success and failure recovery sequences. 

The Operational Analysis Model is based upon a 
generalized mission scenario. The mission scenario is defined 
as the timeline sequence of all possible events (within the 
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desired level of detail) that can occur from start to comple- 
tion of a mission. The success -oriented path through the 
mission scenario is the mission profile. A typical mission 
profile based upon the Space Shuttle-Tug Transportation System 
is illustrated in Figure 2.6. Note that/ in the general sense, 
this represents a multi-stage launch system* As will be dis- 
cussed in following paragraphs , the multi-stage launch system 
may consist of various combinations of recoverable and re- 
usable and/or expendable stages. This is possible by providing 
appropriate values for the various recovery probabilities. 

The mission scenario includes both the success-oriented paths 
and all of the pertinent failure recovery paths as well. The 
mission scenario can be depicted as a series of nodes con- 
nected by branches. Each node is a decision point representing 
a group of events. Emanating from each node are branches for 
the success and failure recovery paths. The probability of 
departing the node on one branch or another depends upon the 
probability of success (or failure) of the events represented 
by the node. The failure recovery paths must ultimately pro- 
vide a route to mission completion ( i . e . , p laoi ng NA . . space- 

1 f J 

craft successfully into orbit) . This mission scenario is 
illustrated in Figure 2.1 as a logic flow diagram where the 
nodes are represented by the diamond-shaped boxes and the 
branches as lines with major cost-associated events given in 
the rectangular boxes . 







Figure 2.7 Space Shuttle - Space Tug Space- 
craft Placement Mission Scenario 
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It should be noted that the logic diagram of the 
mission scenario is based upon the Space Shuttle - Space Tug 
spacecraft placement mission. The reason for this is that 
a large class of space transportation systems and associated 
scenarios form a subset of that shown in Figure 2.7.' For 
example, the scenario for an expendable launch system is 
achieved simply by setting all of the appropriate probabi- 
lities (for example, probability of orbiter (upper stage) 
recovery) to zero. 

The Operational Analysis Model establishes the 
probability distributions of the number of events which occur 
as a result of the requirement to place NA^ ^ spacecraft 
successfully in orbit. This is accomplished through the use 
of Monte Carlo techniques and the associated random sampling 
to determine which of the branches are to be followed at each 
of the node s . 

Annual Cost Model ; The purpose of the Annual Cost 
Model is to establish the probability distribution of the 
annual recurring cost associated with the service-type of 
mission. The annual cost is established as the sum of the 
various unit costs multiplied by the associated number of 
events. The unit costs are adjusted individually on an 
annual basis to allow for learning effects as follows: 

C . = UC X i^ 

1 , '■■■'. ■ 

a = -1 + [log^^CAL-l. 699] /. 3(|)1 
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th 

where is the cost in the i year based upon a first unit 
cost of UC. a is the learning exponent and CAL is the cumu- 
lative average learning rate (%), When CAL=1 (|h|)/ there is no 
learning. When CAL=S4>, the cost is reduced to 8<j)% of its 
value every time the number of years doubles. 

The unit costs are the result of sampling of the 
unit cost uncertainty profiles and ranges of uncertainty and 
the number of events is the result of the operational analysis 
simulation. By repeating the process a large number of times, 
the annual cost probability distributions are established. 

The final step, after all time periods have been 
considered, is the establishment of the present value proba- 
bility distribution and the expected value and standard 
deviation of PV. 

2.3 An Application Example 

To illustrate the use of the SATIL 2 program and to 

demonstrate typical results obtained from the explicit con- 

■ . * 

sideration of reliability and cost uncertainties, a typical, 
though hypothetical, service mission is considered over a ten- 
year time period. It should be noted that this mission is 


Cost uncertainties have not been considered in the re- 
sults shown so that the impact of reliability consider 
ations can be clearly seen. 
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strictly hypothetical and is not intended to represent an 
operational SEASAT program. It is planned that the data set 
needed for the analysis of operational SEASAT alternatives 
will be compiled jointly by NASA, JPL, and ECON during FY-1976 
Estimates of costs using the SEASAT data set and the SATIL 2 
program will be prepared during FY-1976. 

The basic mission is described by the sensor demand 
function given in Table 2.1, where the assumed demand for five 
different sensors in geocentric orbits is illustrated for a 
10-year period. Table 2.2 illustrates the assumed sensor/ 
spacecraft mix based upon four spacecraft types. Spacecraft 1 
can accomodate Sensors 1 and 2, Spacecraft 2 can accommodate 
Sensors 2 and 3, and sc forth. The assumed reliability param- 
eters are mean-time to failure ? expected wearout life, and 
standard deviation of wearout life. 

A typical specification of launch systems techno logy 
is illustrated in Table 2.4, where the probability of success 
of performing major operations is shown for five launch tech- 
nologies. These technology specifications are representative 
of the Thor-Delta, Atlas -Cen taur , Space Shuttle and Interim 
Upper Stage, Space Shuttle and Improved Interim Upper Stage, 
and Space Shuttle and Reusable Space Tug Transportation 
Systems, respectively. Table 2.5 describes the applicability 
of these launch technologies to placing each of the four 
spacecraft into orbit as a function of time. For example. 




Table 2.1 Assumed Sensor Demand 


Year 


0 0 1 2 2 2 2 
0 0 1 1 2 2 2 


Table 2.2 Sensor/Spacecraft Mix 


Spacecraft 


Sensor 



Sensor Reliability Assumptions 


Reliability Item 

Sensor 

1 

2 

3 

4 

5 

Mean Time to Failure (yrs) 

5 

7 

7 

5 

5 

Expected Wearout Life (yrs) 

5 

5 

5 

5 

5 

Std. Dev. of Wearout (yrs) 

1 

1 

1 

1 

1 
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Table 2.4 Reliability Soecif ications of 

Launch Technology 




Launch Technology 

.* 


Probability of Success for 

1 

2 

3 

4 

5 

Booster burn 


.95 

.95 

1.0 

1.0 

1.0 

' Booster recovery 

given success 


.0 

.0 

.0 

. 5 

Booster recovery 

given failure 

.91 

.0 

.0 

.0 

.0 

Or biter achieves 

low earth orbit 

.95 

.95 

1.0 

1.0 

1.0 

Orbiter recovery 

given abort to orbit 

.95 

; .0 

1.0 

1.0 

1.0 

Orb iter recovery 

given booster failure 


.0 

f — 1 

1.0 

1.0 

Orbiter recovery 

from low earth orbit 

.95 

.0 

1.0 

1.0 

1.0 

Upper stage (i.e. , space tug) checks 
out successfully in low earth orbit 

1.0 

1.0 

.95 

.95 

. 95 

Upper stage transfers to s/c placement 
orbit 

.95 

.92 

.95 

.95 

. 95 

Upper stage transfers from terminal 
orbit to low earth orbit and rendez- 
vous with orbiter 

.0 

.0 

.0 

. 9 

.95 

Upper stage reacquiring s/c which does 
does check-out in terminal orbit 

.0 

.0 

.0 

.0 

. 9 


.... 


s/c 

type 






1 2 

4 

4 


s/c checks-out in 

terminal orbit 

• 

95 .9 

1.0 

1.0 


it 

Probabilities of 

1.0 imply greater than 

0.99 

chance 

of success . 
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it is assumed that Spacecraft 1 will be placed in orbit by the 
Thor-Delta during the first three years, by the Space Shuttle 
and Interim Upper stage during the next three years, etc. 

Table 2.6 illustrates the assumed costs for sensors, 
launch technologies, and spacecraft buses. The model allows 


Table 2.5 

Launch Technology Applicable 

per 

Spacecraft 


Spacecraft 





Time 

(years ) 




D 

2 

3 

B 

5 

6 

7 

■ -'8 

,9 

10 

1 

1 

1 

1 

3 

3 

3 

4 

,,, 

5 

5 

2 

2 

2 

2 

3 

3 

3 

4 

4 

5 

5 

3 

2 

2 

2 

3 

3 

3 

4 

4 

5 

5 

4 

2 

2 

2 

3 

3 

3 

4 

4 

5 

5 



SSS9B 

6 Cost Estimates* 

C$mi 

llions) 


Sensor 


Cost 

Launch 

C os t 


Bus 

: Cost 



C$M) 

Technology 




($M) 

1 


15 

1 

10 


1 

15 

2 


7 

2 

15 


2 

30 

3 


7 

3 

B 

1 



4 


5 

4 

H 

1 



5 


40 

5 

15 


. ..t. 


* 

Costs 

assumed 

known with certainty 
















all of these costs to be treated as uncertainty variables 
(ranges plus probability distributions). However, since 
emphasis is being placed upon the effects of unreliability, 
all costs have been assumed to be known with certainty (i.e., 
single values ) . 

The Mathematical Model develops, based upon the 
Monte Carlo methodology, probability distributions for all 
pertinent events and costs and present value of cost. Typical 
event statistics (based upon the data presented in the previous 
tables) are shown in Table 2.7, where the probability density 
function of launch attempts as a function of time is presented. 
This data encompasses the effects of both transportation sys- 
system and spacecraft (bus and sensors) unreliability. Table 
2.8 illustrates the summary statistics, in terms of expected 
value and standard deviation of annual launch attempts, when 
all random failure rates (bus and sensors) are zero and all 
launch reliabilities are unity. Only bus wearout failures are 
considered (sensor wearout life is assumed in excess of bus 
wearout life). In effect. Table 2 . 8 i 1 lus trates the number of 
launch attempts that would be considered when allowing for 
"planned" replacement. The effect of random failures is 
immediately apparent by comparing the summary statistics in 
Tables 2 . 7 and 2 . 8 . 

Table 2.9 illustrates the probability density func- 
tion of total annual recurring costs. The variability of the 













Table 

2.7 Probability of 

Indicated 

Number 

of 

Launch Attempts 


Quantity 










10 

.0 

.0 .0 

. 0 

. 0 

. 0 

.01 

. 0 

. 0 

. 01 

9 

. 0 

.0 ,02 

. 0 

, 0 

, 0 

. 01 

. 0 

.01 

. 0 

8 

.0 

.0 .03 

. 0 

.03 

. 0 

. 02 

.01 

.02 

.02 

7 

- 0 

.0 .06 

. 0 

.05 

.07 

.03 

.03 

.02 

.05 

6 

.0 

.0 .07 

. 0 

, 10 

m O / 

.07 

. 10 

.04 

.07 

5 

. 0 

.01 . 25 

. 0 2 

. 16 

. 11 

. 15 

. 11 

. 08 

. 11 

4 

. 0 

.02 .34 

.03 

, 33 

. 16 

.18 

.13 

.27 

.17 

3 

. 01 

.10 .23 

.13 

. 30 

.19 

. 2 6 

. 17 

. 19 

.25 

2 

,15 

.30 .0 

. 26 

.03 

. 19 

.22 

, 3 0 

.23 

.22 

1 

,84 

.57 .0 

. 50 

. 0 

. 16 

. 04 

. 13 

. 12 

.06 

0 

, 0 

.0 ,0 

. 06 

. 0 

. 05 

.01 

.02 

.02 

.04 


1 

2 3 

4 

5 

6 

7 

8 

9 

10 





Year 






Expected 
No . 

1.17 

1.60 4.56 

1 . 63 

4.27 3 

. 21 3 

.71 

3 .20 

3 . 26 

3.51 

Std . 

Dev - 

. 40 

. 82 1.40 

1.0 2 

1.36 1 

.89 1 

. 79 

1 . 78 

1.72 

1.89 









Table 2.9 

Probabili-ty of 

Indicated Total Launch & Spacecraft Costs 


Cost Range 
( Millions 
of Dollars ) 











480.00 500.00 

0. 0 

0.0 

0.0 

0,0 

0.020 

0.0 

0,0 

0.0 

0.0 

0.0 

460.00 480.00 

0,0 

0. 0 

0 . c 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.010 

440. 00 460.00 

0.0 

0.0 

0.0 

0.0 

0.020 

0.0 

0.0 

0.0 

0.010 

0.0 

420.00 440.00 

0.0 

0.0 

0.0 

0.0 

0.010 

0.010 

0 .020 

0.0 

0.010 

0.010 

400. 00 420.00 

0.0 

0.0 

0 . 0 

0 . 0 

0.0 

0.010 

0 .010 

0.010 

0.0 

0.010 

380 . 00 --- 400. 00 

0.0 

0.0 

0.020 

0.0 

0.020 

0.030 

0.020 

0.010 

0.0 

0.020 

360 . 00 --- 380.00 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 020 

0.020 

0,030 

0.020 

0.020 

340 . 00 --- 360. 00 

0.0 

0,0 

0.030 

0.0 

0.020 

0.050 

0,020 

0.040 

0.030 

0.0 

320.00 340. 00 

0.0 

0. 0 

0.010 

0.0 

0.040 

0.020 

0,020 

0.010 

0.010 

0. 050 

300 . 00 320.00 

0.0 

0.0 

0.050 

0.0 

0.060 

0.030 

0.050 

0 . 020 

0.030 

0.060 

280 . 00 300 . 00 

0.0 

0.0 

0.0 

0.010 

0.080 

0.070 

0.080 

0. 110 

0,030 

0.040 

260.00 280.00 

0.0 

0.0 

0.070 

0 . 0 

0.03 0 

0.010 

0.030 

0 .050 

0.080 

0.050 

240 . OQ 260.00 

0.0 

0.0 

0.0 

0.010 

0.14 0 

0.080 

0.050 

0.020 

0 . 0 10 

0.030 

220 . 00 ; 240. 00 

0.0 

0.0 

0.150 

0.020 

0,130 

0.060 

0 . 100 

0 . 080 

0.120 

0.060 

200 . 00 --- 220.00 

0.0 

0 . 0 

0.100 

0 . 0 

0.020 

0.0 40 

0.010 

0.010 

0.050 

0.010 

180.00 200.00 

0.0 

0.010 

0.0 

0.050 

0 . 220 

0.120 

0.150 

. 0.050 

0.090 

0. 160 

160.00 180.00 

0.0 

0.020 

0.340 

0.0 

0.020 

0.030 

0.050 

0.030 

0.030 

0.020 

140 . 00 - — 160. 00 

_ jO . GI O 

0.0 

0.0 

0 .070 

0.140 

0.050 

0.120 

0.030 

0.080 

0.120 

120. 0 0--- 140 ,. Q ' 3 ' ' 

0 . 0 

0.040 

0. 230 

0.030 

0.0 

0.050 

0.060 

0.110 

0. 130 

0.080 

100 . QG --- 120. 00 

0 . 0 

0.060 

0.0 

0.040 

0.0 

0.0 

0.030 

0.040 

0.040 

0.060 

-■"■ 80 . 00 --- 100. 00 

0.150 

0.150 

0.0 

0 . 200 

0 . 030 

0.110 

0 . 100 

0.190 

0.110 

0.090 

60 . 00 - — 80.00 

0.0 

0. 150 

0.0 

0 . 020 

0 . 0 

0.030 

0,010 

0.010 

0.010 

0.0 

40. 00 60.00 

0 . 840 - 

0 . 0 

0 . 0 

0 .490 

0.0 

0.130 

0.040 

0.130 

0.090 

0.060 

20 . 00 — - 40.00 

0. 0 

0. 570 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 20.00 

0.0 

0 . 0 

0.0 

0 . 060 

0 . 0 

0.050 

0.010 

0.020 

0.020 

0.040 

Year 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Expected Costs 

55.0 

61.8 

202 . 0 

82.8 

241.1 

187.8 

203.7 

176.7 

180.8 

192.7 

Std . Dev . 

18.9 

34.1 

62.2 

56.9 

80.1 . 

110.8 

93.7 

103 . 7 

95.1 

103.4 
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costs and quantities (with cost uncertainties set equal to 
zero) is due to system and subsystem unreliability. 

The risk profile of present value of recurring cost 
(i.e., probability of PV exceeding specified values) is illus 
trated in Figure 2.8. Since this risk profile is approxi- 
mately normal/ the expected, most likely and median values, 
are the same, i.e., $1037 million with a standard deviation 

of $132 million. If unreliability of launch systems and 
random spacecraft and bus failures are not considered (as in 
Table 2.8) / the expected present value is $545 million with a 
standard deviation of $35 million. 

The importance and economic value of unreliability 

is addressed in Figure 2,9 which indicates the effect of MTBF 

and expected wearout time (o^ held constant) . The value of 

achieving different MTBFs and apparent. Moving from 

an expected wearout time of M to M_ , and on MTBF to 

^ / Ji- ^ 

MTBF^ , results in an expected economic value of m^-m^ and a 
reduction in ris;k of a^-a^ . Thus, a development program 
should be undertaken which would move the technology from 
Point A to Point B if the expected present value of the cost 
of the development program were less than m^-m^ . More 
precisely, it is necessary to develop the probability distri- 
bution of the present value of the development program and 
compare it with the probability distribution of the present 
value difference between technologies indicated as points A 
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and B. Whether or not the development program to move from 
A to B should be undertaken depends upon the chance that the 
present value of development costs will exceed the reduction 
in the present value of the recurring costs and the decision 
makers risk aversion preferences. To illustrate the above 
quantitatively/ m^-m^ is equal to :$24 million and cr^~g3 is equal 
to $12 million when bus expected wearout time is increased from 
f our- to-eight years with a constant MTFB of 15 years, and the 
other parameters are as given in the previous tables. 

The importance and economic value of unreliability 

is addressed in Figure 2.9 which indicates the effect of MTBF 

and expected wearout time {O ^ held constant). The value of 

achieving different MTBFs and M^s is apparent. Moving from 

an expected wearout time of M_ , to _ and MTBF, to MTBF^ 

'f/1 f,3 1 3 

results in an expected Economic value of ^^d a reduction 

in risk of • Thus ,, a development program should be under- 

taken which would move the technology from point A to point B 
if the expected present value of the cost of the development 
program were less than m^-m^. More precise ly , it is necessary 
to develop the probability distribution of the present value 
of the development program and compare it with the probability 
distribution of the present value difference between technolo- 
gies indicated as points A and B. Whether or not the develop- 
ment program to move from A to B should be undertaken depends 
upon the chance that the present value of development costs 
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will exceed the redu.ction in the present value of the recur- 
ring costs and the decision makers risk aversion preferences. 
To illustrate the above quantitatively, ™l”'^3 equal to $23 
million and is equal to $12 million when bus expected 

wearout time is increased from f our- to-eight years with a 
constant MTBF of 15 years and the other parameters are as 
given in the previous tables. 
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3. INPUT DATA REQUIREMENTS, DEFINITIONS 
AND DATA SHEETS 

The SATIL 2 data requirements, definitions, and input 
data sheets are presented in the following pages. The Fortran 
name of the input data variable is presented along with the 
definition of the variable. Each variable has a reference 
number associated with it. This reference number appears on 
both the input data specification sheet and the input data 
sheet, thus allowing for easy cross-referencing. 

A range of valid data is given for each of the input 
variables. If the input data does not lie within the valid 
range, an error message is printed and the program execution 
terminates. In order to simplify the inputting of data, all 
variables are initialized to "default" values by the program. 
The default values are indicated on the input data specifi- 
cation sheet. If input data is omitted for a variable, the 
program will use the default value. Therefore,,; if the input 
data for a variable corresponds to the default value for that 
variable, it is not necessary to input the data to the program. 

The input data sheets contain the Fortran name of 
each of the input variables and space for the input data. 

The input data system utilizes the Fortran "NAMELIST" feature. 
This feature requires that the name of each variable be input 
along with an equal sign and then the data. Each data 
term must be followed by a comma. This is all taken care of 
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by the input data sheet as well as the variations for integer 

and real variables. It is only necessary to place the input 
J * 

data in the appropriate spaces on the data sheets and thence 
to keypunch one card (and continuations if necessary) for 
each row shown on the data sheets. Everything in the ruled 
boxes is to be keypunched. 

I, 



I 



w 



VARIABLE 

NAME 

1. MAXN 

2 . MAXLSC 

3. MAXES 

4. MAXNB 

5 . MAXLB 






'W 


INPUT DATA SPECIFICATION (SATIL 2) 


DEFINITION 


RANGE 


DEFAULT 

VALUE 


Number of time periods (years) to 1-15 

be considered. 

Number of spacecraft types to be 1-5 

considered. A spacecraft consists 
of a bus, made up of a number of 
subsystems , and one or more sensors . 

Number of sensor types to be con- 1-10 

sidered. This is the total number 
of different sensors used by the 
MAXLSC spacecraft. 

Maximum number of spacecraft buses 1-5 

to be considered. A bus consists 

of a number of support subsystems 

(i.e., AVCS , TTC, power subsystem, 

etc.) all of which must function 

properly for the sensors (i.e., the 

payload) to be operational. 

Maximum number of spacecraft subsys- 1~5 

terns which comprise the bus and whose 
reliability are to be explicitly 
considered. 

Maximum number of launch system con- 1-10 

figurations to be considered. A 
launch system consists of a set of 
cost and reliability data in terms of 
spacecraft type. 


1 

1 


1 


1 




6 . 


MAX IE 


1 


cn 

4 ^ 


VARIABLE 

NAME 


DEFINITION 


RANGE 


DEFAULT 

VALUE 


7 . 

MAXR 


Number of Monte-Carlo simulation 
runs to be performed (even multiples 
of 100) . When only a single run is 
to be performed, enter MAXR = 0- 

Unlimited 

1 

8. 

MAX I - 


Maximum number of rows to be printed 
in "Probability of Indicated Quantity" 
output table. 

1-40 

25 

9. 

NO I NT 


Number of rows to be printed in 
"Probability of Indicated Costs" and 
"Probability of Present Value of 
Cost" output tables. 

1-40 

25 

10. 

NSEN__ 

(N) 

Number of operating sensors of type 
LS required in orbit in year N. 

Values must be specified for N < MAXN. 

0-10 

for 

0 

all N 

11. 

IMIX 

(LSC) 

Specification of sensors of type* LS 
per spacecraft type (LSC) . When equal 

0 or 1 

for 

0 

all LSC 


to 1, the LS sensor is on-boar<3 the 
LSC spacecraft. When a sensor is not 
on-boaifd a specific LSC spacecraft type, 

the corresponding value of IMIX is 0 . 

Values must be specified for 1 < LS < 

MAXLS and 1 < LSC < MAXLSC. 

12. IBMIX (LSC) Identity of the spacecraft bus* (i.e., 1-5 1 

value of NB) to be used to form the for all LSC 

LSC spacecraft. Values must be speci- 
fied for 1 < LSC < MAXLSC. 


* The general subscript nomenclature is illustrated in Figure 3.1. 



Figure 3,1 General Subscript Nomenclature 
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VARIABLE 

NAME 


DEFINITION 


RANGE 


13 . AMTBF (NB) 


Mean time before failure (years) 
for random or chance failures of 
the LB subsystem which comprises 
the NB bus. Specify values for 
1 < LB < MAXLB and 1 < NB < MAXNB 


14. STAND (NB) 


Standard deviation (years) of the 
wearout life of the LB subsystem 
of the NB bus. Specify values for 
1 < LB < MAXLB and 1 < NB < MAXNB. 


15. AMEAN (NB) 


Expected or mean wearout life (years) 
of the LB subsystem of the NB bus. 
Specify values for 1 < LB < MAXLB 
and 1 < NB < MAXNB. 


16 


AMTBS (LS) 


Mean time before failure (years) for 
random or chance failures of the LS 
sensor. Specify for 1 < LS < MAXLS. 


17. STANDS (LS) 


Standard deviation (years) of the 
wearout life of the LS sensor. 
Specify for 1 < LS < MAXLS . 


18. AMEANS (LS) 


Expected or mean wearout life (years) 
of the LS sensor. Specify for 
1 < LS < MAXLS . 


. 01 - 100 . 


.01-100 


. 01 - 100 . 


. 01 - 100 . 


. 01 - 100 . 


. 01 - 100 . 


DEFAULT 

VALUE 


5.0 

for all NB 


0.5 

for all NB 


5.0 

for all NB 


5.0 

for all LS 


5.0 

for all LS 


5.0 

for all LS 






















w 
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VARIABLE 

NAME 


DEFINITION 


RANGE 


19. AMAXC (LS) 


Maximum estimated unit cost (10 ?) 
of the LS sensor. This establishes 
the upper limit for the rai(ge of 
cost uncertainty. Specify for 
1 < LS < MAXLS. 


. 0 — 1000 . 


20 . AMINC (LS) 


Minimum estimated unit cost (10 $) 
of the LS sensor. This establishes 
the lower limit for the range of 
cost uncertainty. Specify for 
1 < LS < MAXLS . 


.0-1000 


21. IPT (LS) 


Name of the uncertainty profile* 

(i.e. , probability density function) 
to be associated with the range of 
uncertainty established by AMAXC 
and AMINC for the LS sensor. This 
establishes the subjective probability 
distribution of cost. Typical uncer- 
tainty profiles are illustrated (see 
Figure 2-2) and may be modified by 
changing values of PPP . Specify for 
1 < LS < MAXLS. 


1-30 


22. CALC (LS) 


Cumulative average learning rate (%) 
for the LS sensor. Assumes that costs 
are reduced by 100. -CALC percent every 
time the number of years from first 
use doubles. Specify for 1 < LS < MAXLS 


01-100 


The default uncertainty profiles are illustrated in Figure 2-2 


DEFAULT 

VALUE 


.0 

for all LS 


.0 

for all LS 


16 

for all LS 


90 

foi all LS 


cn 















VARIABLE 

NAME 


23. AMAXB 


24. AMINB 


25. IPTB 


26. CALB 


DEFINITION RANGE 


Maximum estimated unit cost (10 $) .0-1000. 

of the NB spacecraft bus. This 
establishes the upper limit for 
the range of cost uncertainty. 

Specify for 1 < NB < MAXNB. 

g 

Minimum estimated unit cost (10 $) .0-1000. 

of the NB spacecraft bus. This 
establishes the lower limit for 
the range of cost uncertainty. 

Specify for 1 < NB < MAXNB. 

Name of the uncertainty profile 1-30 

(i.e., probability density func- 
tion) to be associated with the 
range of uncertainty by AMAXB and 
AMINB for the NB spacecraft bus. 

Specify for 1 < NB < MAXNB. (See 
IPT (LS) ) . 

(NB) Cumulative average learhihg rate .01-100. 

(%) for the NB spacecraft bus. 

Assumes that costs are reduced by 
100, -CALB percent every time the 
number of years from first use 
doubles. Specify for 1 < NB < MAXNB. 


(NB) 


(NB) 


DEFAULT 

VALUE 


.0 

for all NB 


.0 

for all NB 


16 

for all NB 


90. 

for all NB 


VARI?i 

BLE lil'-yiE 

DEFINITION 

RANGE 

DEFAULT 

VALUE 

27. 

MiAXM 

Maximum estimated spacecraft re- 
furbishment and maintenance cost 
expressed as a percentage of total 
spacecraft (including bus and sensors) 
cost. This estabishes the upper limit 
for the range of cost uncertainty. 


.0 

28. 

AWINM 

Minimum estimated spacecraft refurbish- 
ment and maintenance cost expressed 
as a percentage of total spacecraft 
(including bus and sensors) cost. This 
establishes the lower limit for the 
range of cost uncertainty. 

.0-100. 

•0 

29. 

IPTB 

Name of the uncertainty profile (i.e. 
probability density function) to be 
associated with the range of uncer- 
tainty by AMAXM and AMINM. (see IPT 

(LS)). 

1-30 

16 

30. 

CALK 

Cumulative average learning rate (%) 
for the spacecraft refurbishment and..- 
maintenance cost. Assumes that costs 
are reduced by Ij20.-CALR percent every 
time the number of years from first 
use doubles . 

.01-100. 

90. 



VARI ARLE 
NAME 


DEE IN IT ION 


RANGE 


31. AMXLC_(IE) Max^muiri estimated launch cost ,0-1000. 

(10 $) for the LSG spacecraft 
with the IE launch system. This 
establishes the upper limit for 
the range of cost uncertainty. 

Specify for 1 £ IE £. MAXIE. 

32. AMNLC_(IE) Minimum estimated launch cost .0-1000. 

(10 $) for the L3C spacecraft 
v/ith the IE launch systen^. 

Specify for 1 5 IE < MAXIE. 

33. IPTLC_(IE) Name of the uncertainty profile 1-30 

(i.e., probability density func- 
tion) to be associated with the 
range of uncertainty by AMXLC_ 
and AMNLC_ for the IE launch 
system. Specify 1 £ IE < MAXIE. 

34. CALLC_(IE) Cumulative average learning rate .01-100. 

{%) for launching the LSC space- 
craft v;ith the IE launch system. 

Assumes that costs are reduced by 
100.~CALLC_ pe“cent every time the 
number of years from first use 
doubles. Specify for 1 £ IE £ MAXIE- 

35. LNCHS_(N) The name of the launch system (value 1-10 

of IE) to be used in year N for 
launching the LCS spacecraft. 

SpGcify for 1 £ N £ MAXN. 


DEFAULT 

VALUE 


.0 

for all IE 


.0 

for all IE 


16 

for all IE 


90. 

for all NB 


1 

for all N 
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VARIABLE 

NAME 


DEFINITION 


RANGE 


DEFAULT 

VALUE 


36. ORBCAp (LSC) Indicator of Orbiter capability Y or N N 

for placing the LSC payload directly for all LSC 

into the desired orbit. When ORBCAP 
= Y it is implied that the orbiter 
is capable of placing the payload 




into the desired orbit. When ORBCAP 
7 ^ Y it is implied that a "space tug" 
or upper stage is required. Specify 
for 1 ^ LSC ^ MAXLSC. 




37 . 

PPS (LSC) 

Probability of the LSC spacecraft 
functioning properly when placed 
in final desired orbit (wear-in 
period) . Specify for 1 < LSC < 

MAXLSC. 

.0-1.00 

1.00 
for all 

LSC 

38 . 

PBS (IE) 

Probability of a booster success 
per launch for the IE launch system. 
Specify for 1 < IE < MAXIE. 

.0-1.00 

1.00 
for all 

IE 

3 9 . 

PBRS (IE) 

Probability of a booster recovery 
given a booster launch success. 
Specify for 1 < IE < MAXIE. 

.0-1.00 

1.00 
for all 

IE 

40. 

PBRF (IE) 

Probability of a booster recovery 
given a booster failure or malfunc- 
tion. Specify' for 1 < IE < MAXIE. 

.0-1.00 

1.00 
for all 

IE 

41 . 

POS (IE) 

Probability of an orbiter success 
(i.e., no abort prior to placing 
spacecraft and/or upper stage and 
spacecraft in orbit) . Specify for 

.0-1.00 

1.00 
for all 

IE 


1 < IE < MAXIE. 


o 



VARIABLE 
NAME - 


DEFINITION 


RANGE 


42. POROF (IE) Probability of an orbiter recovery , 0 - 1.00 

given an orbiter abort. Specify 
for 1 < IE„ £ MAXIE. 

43. PORBF (IE) Probability of an orbiter recovery . 0 - 1.00 

given a booster launch failure. 

Specify for 1 £ IE _< MAXIE. 

44. POR (IE) Probability of an orbiter recovery . 0 - 1.00 

given an otherwise successful 
flight. Specify for 1 £ IE £ 

MAXIE. 

45. PPMO (IE) Probability of propulsion module . 0 - 1.00 

(space tug or other upper stage) 
checking out successfully (wear-in 
period) prior to leaving vicinity 
of orbiter. Specify for 1 £ IE £ 

MAXIE. ~ ~ 

46- PPMP (IE) Probability of propulsion modules . 0 - 1.00 

(space tugi or other upper stage) 
transferring to and placing space- 
craft in desired orbit. Specify for 
1 £ IE £ MAXIE. 

47. PPMRP (IE) Probability of propulsion module .0-1.00 

(space tug or other upper stage) 
successfully returning and rendez- 
vousing with the orbiter. Specify 
for 1 £ IE £ MAXIE. 

48. PPLRR (IE) Probability that a spacecraft which .0-1.00 

does not check out successfully in 
terminal orbit is reacquired by the 
propulsion module (space tug or other 
upper stage) for return to low earth 
orbit. Specify for 1 £ IE £ MAXIE. 


DEFAULT 

VALUE 


1.00 

for all IE 


1.00 

for all IE 


1.00 

for all IE 


1.00 

for all IE 


1.00 

for all IE 


1.00 

for all IE 


1.00 

for all IE 



RANGE 


DEFAULT 

VALUE 
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VARIABLE 

NAME 


DEFINITION 


49. 

MAXM 

(LSC) 

Maximum number of spacecraft of 
type LSC which will be operational 
at any point in time. 

1-20 

10 

50. 

NN 


Number of time periods to be con- 
sidered for the optimization of 
spacecraft launch schedule. 

1-MAXN 

5 

51. 

DR 


Discount rate, expressed as a 
percentage, used in the present 
value of cost computations. 

.0-100. 

10.0 

52 . 

DETC 


If it is desired to have detailed 
cost reports (i.e., costs by 
spacecraft type) printed then 
DETC = ' Y ' . 

Y or N 

N 

53 . 

PPA__ 

(lA) 

Uncertainty profile data; i.e., 
probability associated with the 
lA interval of profile IB where 
1 1 1 5 and 1 < IB ^ 20. 

.0-1.00 

See 

Table 3 


to 



INPUT DATA SHEET 


LSC 




1 

2 

3 

4 

5 

36 . 

ORBCAP = 

1 1 

$ 

r r 

9 

1 t 

9 

f t 

9 

1 1 

9 

37. 

PPS 

• / 

• / 

• / 

• / 

• / 


IE 

123456789 10 


38. 

PBS 


• 9 

• 9 


9 


9 

.. *. 

/ 


/ 


9 


9 

• 9 

^ 

□ 

39. 

PBRS 


• 9 

• / 

• 

9 


9 


/ 


9 

'•{ 

9 


9 

• 9 

• / 

40. 

PBRF 


• 9 

• / 


9 


9 


/ 


9 


9 


9 

• / 

• 9 

41. 

POS 


• 9 


• 

9 


9 


9 


9 


9 


9 

• 9 

• / 

42. 

POROF 

_ 

• 9 

_ 

9 


9 

• 

9 


9 


9 

• 

9 

• 

9 

• 9 

. • ./ _ 

43 . 

PORBF 

— 

• / 

• 9 


9 

• 

9 

. 

9 


9 


9 


9 

• 9 

• 9 

44. 

POR 

_ 

m 


■ 

■ 

■ 

■ 

■ 

1 

1 

■ 

■ 

■ 

■ 

1 

H 

m 

45. 

PPMO 

— 

m 

m 

■ 

■ 

■ 

■ 

■ 

1 

1 

■ 

■ 

■ 


1 

■ 

M 

46. 

PPMP 

_ 

• 9 

* 9 


9 


/ 




/ 


/ 


9 

• 9 

• 9 

47. 

PPMRP 


• 9 

• 9 


9 


t 




t 


9 


9 

• / 

• / 

48. 

PPLRR 


• r 

• 9 

. 

9 

. 

9 

. 


. 


4 

9 

V 

t 

• 9 

• 9 


LSC 

1 2 3 4 5 


49. 

riAXM 

1 

J 

9 

. 9 

/ 

9 

50. 

NN 

S 


i 

i . 

I 



51. 

DR 

m 


i • 




52. 


DETC 
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Table 3.1 Initialization (Default) Values 

of Uncertainty Profiles 

PPPXX(IA) 


XX 

1 ■ 

2 

lA 

3 

4 

5 


01 

. 50 

.25 

. 15 

.07 

.03 

w 

02 

.30 

.25 

.20 

.15 

.10 


03 

.30 

.30 

.20 

.13 

.07 

: 

04 

.35 

.40 

.15 

.07 

.03 

; 

05 

.21 

.32 

. 27 

.15 

.05 

f- 

06 

.23 

.30 

.23 

.16 

CO 


07 

.25 

.35 

. 2 5 

.10 

.05 


08 

. 16 

.49 

.24 

.09 

.02 


09 

.12 

.32 

. 3 2 

.17 : 

.07 


10 

. 15 

.34 

.37 

.12 

.02 


11 

.20 

.20 

.20 

.20 

. 20 


12 

. 15 

.22 

.26 

.22 

.15 


13 

. 10 

.25 

.30 

.25 

.10 


14 

.08 

.25 

.34 

. 25 

. 08 


15 

.05 

.25 

.40 

.25 

.05 


16 

.10 

.20 

. 40 

.20 

.10 


17 

.03 

.30 

.34 

.30 

.03 

# 

18 

.05 

.20 

: .50 

.20 

.05 


19 

.03 

.20 

. 54 

. 20 

.03 


20 

.03 

.07 

. 80 

.07 

.03 


t 




w 
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MATHEMATICAL STRUCTURE 


4.1 


Subscript Notations 

N Time (years) 

LS Sensor type 

LSC Spacecraft type 
LB Bus subsystem identifier 

IE Launch system technology 

NB Bus type 

J Column of output data 


4 . 2 Some Notes on Input Data 

The following input data (see Chapter 3) are input 
in the form of linear arrays and form the indicated two-dime 
sional arrays. *, 


Two-Dimensional 


Input Variable 

Array 



Comment 

10 .' 

NSENXX (N) 

NGSEN (LS , N) 

XX 

i s 

value of LS 

11 . 

IMIXXX(N) 

ISMIX (LS ,LSC) 

XX 

II 

" " LS 

13 . 

AMTBF2 (NB) 

AMTBF (LB , NB) 

2 

If 

" " LB 

14 . 

STANDS (NB) 

STANDL (LB ,NB) 

2 

II 

" " LB 

15 . 

AMEAN2 (NB) 

AMEANL (LB , NB ) 

2 

11 

" LB 

31 . 

AMXLCy ( IE) 

AMAXLC (LSC , IE) 

Y 

II 

" " LSC 

3 2 . 

AMNLCY (IE) 

AMINLC (LSC, IE) 

Y 

11 

" ” LSC 

3 3 . 

IPTLCY ( IE) 

IPTLC (LSC , IE) 

Y 

It 

” ” LSC 

3 4 . 

CALLCY (IE) 

CALLC (LSC , IE) 

Y 

If 

" LSC 

3 5 . 

LNCHSY(N) 

LNCHS (LSC ,N) 

Y 


" " LSC 

5 2 . 

PPPXX ( lA) 

PPP (lA, IB) 

XX 

If 

« M Jg 



4. 3 


List of Variables 


The following is a list of the variables used in 
the mathematical formulation of the SATIL 2 program. The * 
indicates input data variables, and the + indicates integer 
variables. The computed variables are defined by the equa- 
tions and operations indicated in the computation functional 
flow given in Section 4.5. 



RELS (LS ,N) 

* 

AMAXLC (LSC , IE) 

+ 

NOOP(LS,N) 


PFS (LS , N) 

* 

AMINLC (LSC , IE) 

+ 

NOSA(LSC,N) 


RELB (NB , N) 

+ * 

IPTLC (LSC , IE) 


IS (LSC ,M) 

+ *' 

MAXLB 


LCELC (LSC , IE) 

+ 

ISS (LS ,LSC ,M) 


AMTBF (LB) 

★ 

CALLC (LSC , IE) 

+ 

lYEAR (LSC ,M) 

*: 

STANDL (LB) 


RANGLC (LSC , IE) 

+ 

NFAIL (LSC) 


AMEANL (LB) 


UCPLC ( I , LSC , IE) 

+ 

ILA (LSC,N) 

*1 

AMTBS (LS ) 


SCCST (LSC ,N) 

* 

PBS ( IE) 

★ 

STANDS (LS) 


UCLC (LSC , IE) 

* 

PBRF ( IE) 

★ 

AMEANS (LS) 

+ 

ITS(LS) 

+ 

IBR(L'SC,N) 

+ * 

MAXLS 


PVF (N) 

+ 

IB (LSC , N) 

+* 

MAXN 

* 

DR 

+ 

IPL(LSC,N) 


LCE(LS) 

+ 

IX (LSC) 

* 

POROF(IE) 


CALC(LS) 

+ 

ITL 

+ 

lOR (LSC ,N) 

+ *; 

MAXIE 

+ * 

LNCHS (LSC , N) 

+ 

10 (LSC, N) 

+ *i 

NN 

+ 

IXB(NB) 

* 

ORBCAB (LSC ) 


PFB(NB,N) 

+ 

IN 

* 

IPM (LSC , N) 

* 

AMAXC(LS) 


LCST(LSC,N) 

* 

PBRS(IE) 

* 

AMINC(LS) 


LCER 

* 

POS(IE) 


RANGE (LS) 

* 

CALR 

* 

POROF ( IE) 

+ * 

IPT(LS) 


LCEB (NB) ; 

+ 1 

IPM (LSC , N) 


UCP ( I , LS) 

* 

CALB (NB) 

* 

PPS(IE) 

* 

PPP(IA,IB) 

+ 

NOSSC(LSC^ : 


POR(IE) 

+ * 

MAXB 

* 

AMAXB (NB) 

+ 

IPLR (LSC ,N) 

+ * 

MAXNB 

* 

AMINB(NB) 

* 

PPMO(IE) 


ACC 


RA.NGEB(NB) 

+ : 

IPMR (LSC ,N) 


RN 

+ * 

IPTB(NB) 


PPMP(IE) 


UC(LS) 


UCPB (I,NB) 

* 

PPLRR ( IE ) 

+ *\ 

MAXR 

* 

AMAXM 

* 

PPMRP(IE) 

+ *\ 

MAXI 

* 

AMINM 


C(K,LSC^N,NR) 

+ * 

NOINT 

+ * 

IPTM 


SUMC (K, LSC , N) 

+ * 

ISMIX(LS, LSC) 


UCPM(I) 


SUMSC (K,LSC,N) 

+ * 

MAXLSC 


UCB(NB) 


CT(K,N,NR) 

+ * 

NOSEN (LS,N) 


UCM 


SUMCT (K, N) 
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SUMJ ( J ,N) 
SUMSJ{J,N) +* 

ANPV (NR) 

SUMLA (LSC ,N) 

SUMPL (LSC , N) 
SUMSPL (LSC , N) 

PPL (II ,LSC ,N) 

* DEIC 

AMX ( K, LSC) 

AMN (K,LSC) 
CSTI(K,LSC) 

P (K,LSC ,L ,N) 

AMXT (K) 

AMT (K) 


RANGEM 
IBMIX (LSC) 

SUMPV 

SUMSLA (LSC ,N) 
CSTIT (K) 

PT (K, L, N) 

PVX 

PVN 

PVI 

PPV(L) 

MCST (K,LSC ,N) 
STDCST (K,LSC ,N) 
ME ( J , N) 

SIDE ( J , N) 


SUMSCT (K,N) 

PE (II, J ,N) 
SUMSPV 

PLA (II .LSC .N) 
MT(K,N) 

STDT (K,N) 
MLA(LSC ,N) 
STDLA (LSC ,N) 
MPL (LSC ,N) 
STDPL (LSC ,N) 
MPV 
STDPV 

+ NOS(LSC,N) 






* Integer 
+ Input Data 




I 


7 9- 


ij 


i t 




t 


B 






i 




4 . 4 Initialization Data 

Initialization of the following variables to the 
indicated values takes place both inside and outside of the 
Monte Carlo loop. List 1 refers to data initialized prior 
to entering the Monte Carlo loop, and List 2 refers to in- 
itialization after entering the Monte Carlo loop. 


List 1 


List 2 


RELS(LS,(Z5) 

= 

1.0 

for 

RELB (NB ,9T) 

= 

1.0 

for 

PFS(LS,N) 

= 

.0 

for 

PFB (NB , N) 

=■ 

.0 

for 

SUMC (K,LSC ,N) 

= 

.0 

for 

SUMSC (K,LSC , N) 

= 

.0 

for 

MLA(LSC,N) 

= 

.0 

for 

MPL 

= 

.0 


SUMJ ( J ,N) 

= 

.0 

for 

SUMSJ(J,N) 

= 

.0 

f or 

PE(II,J,N) 

= 

.0 

for 

SUMLA (LSC ,N) 


.0 

for 

SUMSLA (LSC ,.N) 

= 

.0 

for 

PLA (II, LSC ,N) 


.0 

for 

SDMPL 


.0 

for 

SUMSPL 

= 

.0 

for 

PPL 

= 

.0 

for 

SUMPV 

= 

. 0 


SUMSPV 


.0 


P (K,LSC,L,N) 

= 

.0 

for 

PT (K, L ,N) 

= 

.0 

for 

PPV(L) 

= 

.0 

for 

MI 

_ 

.0 


SCCST(LSC,N) ' 


1 rtf 

for 

LCST(LSC,N) 

= 

10^^ 

for 

I S ( L S C , M ) 

= 

.0 

for 

NOOP (LS , N) 

= 

.0 

for 

ILA(LSC,N) 

s= 

.0 

for 

IB (LSC ,N) 

=s 

.0 

for 

IBR(LSC,N) 

= 

.0 

for 

IPL(LSC,N) 

= 

.0 

f or 

10 (LSC, N) 

= 

.0 

for 

IOR(LSC, N) 

= 

.0 

for 

IPM(LSC,N) 

= 

.0 

for 

IPLR(LSC,N) 

= 

.0 

for 

IPMR(LSC,N) 

S! 

.0 

for 

ISS ( LS , LSC ,M) 

= 

.0 

for 

NFAIL (LSC ,M) 

= 

.0 

for 

ANPV(NR) 

= 

.0 

for 


all LS 
all NB 

all LS and N 

all NB and N 

all K,LSC,&N {K=l,2,3) 

all K,LSC,&N 

all LSC and N 

all J and N 
all J and N 
all II ,J and N 
all LSC and N 
all LSC and N 
all II, LSC,&N 
all LSC and N 
all LSC and N 
all II , LSC, &N 


all K,LSC,L,N 
all K,L,N 
all L 


all LSC and N ( 1<N<MAXN + NN ) 

all LSC and N (1<N<MAXN+NN ) 

all LSC and M 

all LS and N 

all LSC and N 

all LSC and N 

all LSC and N 

all LSC and N ; 

all LSC and N i 

all LSC and N 

all LSC and N 

all LSC and N 

all LSC and N 

all LS , LSC,M 
all LSG,M 
all NR 
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4 . 5 Computation Functional Flow Description 















Sensor Reliability: RELS(LS,N) Assumption is 'that once a sensor is placed in orbit 

it goes Into an operational (i.e., non-dormant) mode 


(2) R E L S ( L S . N ) = ^ ^ ^ ^ - -^ - ^§ .1]-^ 
/TTJv -i STANDS(LS) 


oo 

/exp j"- 

y L 


(x-AMEANS (LS) ) ^ 

2.0 * (STANDS (L5)) 


d X 


t 
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for l<LSC<MAXLSC 

and l<i£<HAXIE ■ 

LCELC (LSC .IE)=- 1 

.+ jlog (CALLC (LSC .IE) ) 


-1.699j/.301 


Launch cost 
learning exponents 



RANGE (LS)= [amAXC (LS) -AMINC{LS)] *, 2 


1=1 



1 = 1 


UCP(I ,LS)=PPP(I ,IPT(LS) ) 


set up probability 
d i s t . of sensor 
cos t s 


UCP( I,LS;)=PPP(6-I. IPT(LS) - 20 


1 = 1+1 
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for HAXN<rf<MAXl! + 'IN and 
for I<lS<HAXLS 
IIOSEN (LS ,Hi=M0bEN (IS .MAXH) 
















Start of Monte-Carlo loop. 

fjR = NR+l - — (7 




' 

f J r. = N K + 1 


Initialize as Per 
List ^2 


AMAXC (LS) 

« 

AMitIC (LS) 


JL-J UC ( L S ) -AMAXC ( LS ) 


]aCC=ACC+UCP (I ,LS) 


ACC > RN 


Sensor cost 
compu tat ion 


juC(LS)=AMlNC (LS)+RANGE(LS)i'.- [I- I] 

^[rn-acc+ucp(i,ls)> 


LS=HAXLS 


*Note: RN implies the generation 

of a random number. 


OBiGWAb PAGiS 
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^ IX(LSC)S»«1AXN ISLSCSMAXLSC , f r « / tum / 1 «r > » 

mrST(LSC.N)-UCB(lBMTir{LSC) )*(N-ITL^iy^^® ^ 

MAXL S 

♦ ISMIX{LS.LSC)*UC(L$) 


LS-I *(N-ITS(LS)*I) 


LC£(LS) 


Sp«cecr«rt 

cott In year N 
(but tentort) 


for MAXN<N<NAXN^NN and I< 1 SC^MAX L SC 
SCCST(ISC.N)»SCCST(ISC .MAXN) 



N ^«S¥ 

<C N-HAXN ^ 


LSC-MAXLSC 


LSC-uSC^ I 


for ISLSCSMAXLSC 
and MAXNcNsMAXN+NN 
lT?T(LSC ,N)-LCST(LSC .MAXN) 
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reproducibility of the 
ORIGINAL PAGE IB POOR 


















I 
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IS(LSC ,M)-0 


lSS(LS,LSC.M)-afor ISIS'S MAXLS 






(See S/C 
launch sched u I c 
notes) 


For l^LSC^hAXlSC 
£oinpu t a t i on of 
NOSAUSC .N) 


0 


NOS{LSC,N)-NOSA(LSC.N) 



N 


Y 
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Notes on Dc t c rn i na t i on of S/C Launch Schedule 

(TTosaTlsc , N~)T 



N-niN 

MAXLSC 


Minimize 

JL, 

PVF(l)y ^(SCCST(LSC , D+LCST (LSC , I) )*nOSA(LSC . I) 



I-N 

LSC- 1 



I 


Subject to: 


M /UIS Q 


N+iiJN 


N00P(LS,N)+> ISMIX (LS ,LSC)* 7 NOSA(LSC.l) 


I 


N 


> NOSEN (LS. I) 


for each I such that N<I<N + NN 
and for 1<LS<HAXLS 


I 


> 


» 


Inputs 


N 

NN 

PVF(I) 

MAXLSC 

SCCST(LSC ,1) 
LCST(LSC ,1) 
N00P(LS,N) 

ISMIX(LS.LSC) 
NOSEN (LS ,I) 


start of the time horizon (years) 
number of time periods (years) to be considered 
present value factor 
maximum number of S/C types 
cost of the LSC S/C In year 1 
cost of the launch of the LSC S/C in year I 
number of sensors of type LS operating successfully 
at the end of year N 

sensor types available on the LSC spacecraft (0 or I) 
number of sensors of type LS required in operation 
I n yea r I . 


Output 

N0SA(LSC,I) number of S/C of type LSC to be placed successfully 
In orbit in year I. 


» 


I 
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n-tNCHSUsc .»i) 


UAU$C.N)-ILA(LSC.N)»I 


T booster 

r«ni 

AN<?8S(II ) >- 


booster 
recovered 
V _ 


RN<?CRF (II ) 


II BR(LSC.N)-: BR(LSC .N) »l 


^ b90Ster not 
recovered 


IB(L$C .»)• IBUSC ,N)*I 


booster not 
recovered 


orbiter not 
recovered 


RN<PBRS(I1) 


IPL(LSC.N)-1PI(LSC.N)+I 


RM<P0R0F(II) 


booster 

recovered 


BR(L5C.ri)-ISR(LSC.N) + l 


IB(LSC.N)-IB(LSC.N)«I 


orbiter Y 

recovered X 

I I0R(LSC.N)-IOR(LSC.ri)*l 



Rfi<P0S(ll) 


No orbiter 
eoor t 


or b i ter 
ebor t 


ORBCAP(LSC)-V 


IPM(LSC.fi)*lPfl(LSC ,N)*I 


orbi ter not 
r ecove red 


orbiter V 

recovered 1 

lORUSC .N)«I0R(1.SC ,N) ♦! 


IPL(LSC.M)-IPL(LSC,U)>I 


IO(LSC.H)»lO(LfC.IJ)*l 


ORBCAPdSC) 
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i 


I 


I 






I 



> 




I 
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I i 


L) 


IJ 


I 




i> 
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Pay load functions O.K, 
Y ^ 


RN<PPS(LSC) 


“n Payload does 
not function 


Payload not 
reacquired 



RN<PPiiRR (IX) >— K- 


Y Payload 
r eacqu i red 


RN<PPHRP(II ) 


PM does not return 
PL to orbiter 

V — -Jm- 


Orbiter 

recovered 


Orbiter not 
recovered 


RN<P0R (II ) 


RN<P0R(ll) 


Orbiter no 
r ecove red 
V. N . 


recovered 

I0R{LSC,N)=I0R(LSC. H)>l 

V PL was not reacquired 


I0(LSC,«)-I 

Q.(LSC ,N) + I 




I0R(LSC,N)-I 

0R{:L3Cj,fi) +1 

' 


IPL (LSC,N)»I 

PL (LSC .N) *1 

. f ^ i ■ 


IPM(LSC,N)=IFM(LSC,N)+| 





PM does not return 
to orb I te r 



^ orb I ter 


















2 


C (1 ,LSC ,N ,NR) = ILA(LSC ,N)*LCST (LSC ,N) 


C ( 2 . L S C , N , N R ) = IP L ( L S C , N ) * S C C ST ( L SC , N ) + , . _ „ 

IPLR (LSC , N) ••UCM* . 0 1 *SCCST (L SC , N ) * (N-ITL + 1 ) 


for K=1 ,2,3 

SUMC (K,LSC , M) =SUMC (K, LSC ,N)+C (K, LSC , N , MR) „ 
SUMSC (K.LSC ,N)=SUMSC (K, LSC ,N)+ [C ( K , L SC , N , NR )] 

I 


LSC=MAXLSC 


for K= 1 , 2 , 3 

MAXI SC 

CT{K,N,NR) = N '' C(K, LSC,N,NR) 

LSC=l 

M AXLS C 

SUHCT(K,N)=SUMCT-(K,N) + ^^^C (K, LSC ,N) 

LSC=I 

maxi s c 

SUHSCT (K, N) = SUMSCT (K, N)+\ C.(K, LSC ,N) 

lsT^i 
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ESPBODUCIBILITy OF THE 
pUMiNAL PAGE IB POOR 








































6 






















I 
































DETC = Y 


For all K.LSC and N^TL 
" MCST (K,L5..C.,M) = SUMC (K,LSC ,N)/HAXR 

STDCST(K,L5C,N) = [SUMSC (X,LSC ,N)/f'^XR- 
(MCST (K, LSC , N) ) ^ 


i po r ail J and fy^TL 

HE (J , N)=SUMJ (J , N) /MAXR , 

SIDE ( J , N)= [SUHST (J , N)/HAXR- (ME (J ,N) ) 

For a 1 1 K and N^jTTL 

MT(K,M)-SUMCT(K,N)/HAXR . , 

STDT (K , N)= [SUHSCT ( K. N) /MAXR- (HT (K, H) ) ]^ 

FoiP all LSC and N^TL 

:MLA(LSC-,N) = SUHLA(LSC,fl)/MAXR , 

|STDLA(L5C, N)= [sUM S L A ( L S C , N ) /MAX R- (MLA(LSC,N) ) ‘ 

Fcjr all LSC and N>iTL 

IMPL (LSC , N)=SUMPL (LSC , N) /MAXR , 

|STDPL(LSC.N>[SUrVSPL(LSC ,N)/MAXR - (MPL(LSC,M))]“ 

lMPV-=SUMPV/MAXR „ , 

'STDPV=[SUM5PV/MAXR-(HPV) 3^ 


Print 

Reports 


Fo| N<ITL^al l 
oKpec t cd values 
and s t ii ti d 0 r d devia- 
tions arc equal to 
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4.6 Output Document Description 

The following pages summarize the output documen- 


tation . 



ESTA SIIMENT & MAINTENANCE 
OF A SYSTEM 'J;‘ EAKTII OKUITING SPACECRAFT 


PROBABILITY OF INDICATED QUANTITY 

QUANTITY 


XXX 

X . XXX 

x.xxx 

X . XXX 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

X . XXX 

XXX 

X . XXX 


. 


. 

. 

. 

, 

XXX 

• 


• 

. 

. 

. 

. 

. 

♦ 

• 


• 

• 

• 

. 

• 

. 

• 

. 


. 

. 

. 

, 

. 

. 

XXX 

x.xxx 

X.XXX 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

XXX 

X . XXX 

x.xxx 

X . XXX 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

X.XXX 


LAUNCH 

ADDIT. 

ADDIT. 

ADDIT. 

SPACECRAFT 

BOOSTERS 

ORBITERS 

PM • s 

ATTEMPTS 

BOOSTERS 

ORBITERS 

PM • s 

REQUIRED 

REFURB. 

REFURB. 

REFURB 


REQUIRED REQUIRED REQUIRED 


EXPECTED 
NUMBER XXX. XX 

STD. 

DEV. XXX. XX 


XXX. XX XXX. XX XXX. XX 

XXX. XX XXX. XX XXX. XX 


XXX .XX 
XXX . XX 


XXX . XX 
XXX. XX 


XXX .XX 
XXX .XX 


XXX .XX 
XXX .XX 


X.XXX 


x.xxx 

X.XXX 


SPACECRAFT 

REFURB. 

XXX . XX 
XXX .XX 


iiimill 




a 


9k 
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ESTABLISHMENT & MAINTENANCE 
OF A SYSTEM OF EARTH ORBITING SPACECRAFT 


QUANTITY 

MI 

MI- 1 
MI-2 


3 

2 

1 

0 


EXPECTED 

NUMBER 

STD. DEV. 


PROBABILITY OF INDICATED QUANTITY 


YEAR: N 


PE (MI+1 
PE (MI , 1 
PE (MI-1 


1 , N) /MAXR 
N) /MAXR 
1 , N) /MAXR 


PE (MI+1 , 9 , N) /MAXR 



PE ( II , J , N) 


PE ( 1 , 1 , N) /MAXR 


PE (1 ,9,N)/MAXR 


LAUNCH 

ATTEMPTS 


SPACECRAFT 

REFURB 


ME (1 ,N) ME (9,N) 

STDE ( 1 , N), STDE (9 ,N) 


KJ 


COST 

range 
(MILLIONS 
OF DOLLARS) 

XXXX .XX--XXXX- XX 
XXXX . XX— XXXX . XX 
XXXX . XX — XXXX .XX 


XXXX .XX--XXXX. XX 
XXXX. XX--XXXX.XX 

YEAR 

EXPECTED COST 
STD. DEV. 


ESTABLISHMENT & MAINTENANCE 
OF A SYSTEM OF EARTH ORBITING SPACECRAFT 

I TOTAL LAUNCH S SPACECRAFT COSTS 
TOTAL LAUNCH COSTS 
TOTAL SPACECRAFT COSTS 

These three sets 
of reports are 
to be printed on 
request only, 
i.e., DETC=Y. 


PROBABILITY OF INDICATED 


LAUNCH & SPACECRAFT COSTS 
LAUNCH COSTS (LSC=X) 
SPACECRAFT COSTS (LSC=X) 


(LSC=X)' 


X .XXX 
X . XXX 
X .XXX 

x.xxx 

X .XXX 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx * 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

x.xxx 

X . XXX 

x.xxx 

1 

2 

3 

4 ....... 

MAXN 

xxxx.x 

xxxx.x 

XXXX . X 

xxxx.x • • • . . . . 

XXXX . X 

XXXX .X 

xxxx.x 

xxxx.x 

xxxx.x • 

xxxx.x 
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ESTABLISHMENT f. MAINTENANCE 
OF A SYSTEM OF EARTH ORBITING SPACECRAFT 

{ TOTAL LAUNCH & SPACECRAFT COSTS k=3 
TOTAL LAUNCH COSTS K=1 

TOTAL SPACECRAFT COSTS K=2 

COST 
RANGE 
(MILLIONS 
OF DOLLARS) 


AMT(K) + [NOINT-li *CST1T (k ) --AMT ( K) +NOINT*CSTIT ( K) PT ( K , NOINT , 1 ) /MAX R . - . PT ( K , NOINT , N ) /MAXR 


AMT (K) T2*CSTtT (K) --AMT (K) +3 *CSTIT (K) • 

AMT (K) FCSTIT (K) -- AMT ( K ) +2 *CSTIT ( K ) 

AMT(K) --AMT(K) +CSTIT(K) PT ( K , 1 , 1 ) /MAXR • ♦ • • PT (K , 1 , N) /MAXR 

YEAR 1 .......... N 

ESPECTED COST MT{K,1). . • . . MT(K,N) 

STD. DEV. STDT(K,1) STDT(K,N) 
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MU' 


ESTABLISHMENT (* MAINTENANCE 
OF A system of earth ORBITING SPACECRAFT 


PROBABILITY OF INDICATED 


{ LAUNCH & S 
LAUNCH COS 
SPACECRAFT 


& SPACECRAFT COSTS (LSC-LSC) K»3 
COSTS (LSC=LSC) K“1 
AFT COSTS (LSC=LSC) K*2 


COST 
RANGE 
{MILLIONS 
OF DOLLARS) 

AMN(K,LSC) + [NOINT-11 *CSTI (K,LSC) — AMN(K.LSC) +NOINT*CSTI (K.LSC) P < K , LSC , NOINT, 1 ) /MAXR. . . P (K, LSC , NOINT , N) /MAXR 


AMN (K,LSC)+CSTI(K,LSC) 
ASH(K,LSC) 


-AMN{K,LSC)+2*CSTI(K, LSC) 
-AMN(K,LSC)+CSTI (K,LSC) 


P(K,LSC,1,1, )/MAXR. ,P(K,LSC,1,N)/MAXR 


EXPECTED COST 


MCST{X,LSC.1) .... . HCST{K,LSC,N) 

STDCST(K,LSC,1) • • • * STDCST (K,LSC , N) 
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m 




ESTABLISHMENT & MAINTENANCE 
OF A SYSTEM OF EARTH ORBITING SPACECRAFT 


SPACECRAFT TYPE: XX 


PROBABILITY OF INDICATED 


LAUNCH ATTEMPTS 
SPACECRAFT REQUIRED 


QUANTITY 


XX 

X . XXX 

X . XXX 

X.XXX • 

X.XXX 

XX 

X.XXX 



# 

XX 



• 

• 

XX 

• 

• 

• 


• 

. 

• 

• 


• 


• 

• 


XX 

X.XXX 

X .XXX 

X.XXX 

X.XXX 

YEAR 

1 

2 

3 . • 

• . • • » MAXN 

EXPECTED 

NUMBER 

XX 

XX 

XX • . 

• • • * • XX 

STV . DEV . 

XX 

XX 

XX • ‘ 
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ESTABLISHMENT & MAINTENANCE 
OF A SYSTEM OF EARTH ORBITING SPACECRAFT 

SPACECRAFT TYPE; LSC 

PROBABILITY OP INDICATED LAUNCH ATTEMPTS 


QUANTITY 


MI 

PLA (MI+l ,LSC,1) /MAXR 


MI-1 

PLA(MI,LSC,1)/MAXR 

• 

• • 

MI-2 

PLA (MI-1 ,LSC ,1) /MAXR 

• » 
• • 

• 


: • • 
• • 

PLA (II /LSC , N) 

3 

2 

• 

• 

• 

• 

• 

I 

• 

• 

• 

• 

0 

PLA(1,LSC , D/MAXR- • 


YEAR 

1 . . . . . . 


EXPECTED 

NUMBER 

MLA(LSC,1) . . . * . 



STD .DEV 


STDLA (LSC ,1) 


STDLA (LSC>MAXN) 




XXXXX .XX^ — XXXXX .XX 
XXXXX .XX — XXXXX. XX 
XXXXX .XX — XXXXX .XX 

• • 

# • 

XXXXX .XX--XXXXX .XX 

EXPECTED PV 
STD. DEV. 


X .XXX 
X . XXX 
X .XXX 

• 

x.xxx 

XXXXX .XX 
XXXXX .XX 
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ESTABLISHMENT & MAINTENANCE 


OF A SYSTEM OF EARTH 


PRESENT VALUE 
(MILLIONS OF DOLLARS) 


PVN+ [NOINT-l] *PVI — PVN+NOINT*PVI 



PVN+2*PVI — PVN+3*PVI 
PVN+PVI --PVN+2*PVI 
PVN — PVN + PVI 

EXPECTED PV 
STD. DEV. 


ORBITING SATELLITES 


PROBABILITY OF INDICATED 
PRESENT VALUE OF COST 
(DISCOUNT RATE=DR) 


PPV (NOINT) /MAXR 


PPV (1) /MAXR 


MPV 


STDPV 


12 0 


4.7 Optimum Scheduling Formulation 

The problera is to satisfy the yearly demand for var- 
ious sensor types in every year by launching spacecraft which 
transport a certain subset of the sensor types so that the 
overall program cost is kept to a minimum. If there were no 
random elements such as sensor failures/ then the problem 
could be solved by considering the entire time frame of the 
program. In this analysis, however, the random element is 
incorporated and the decision as to which spacecraft to launch 
must be on a year-to-year basis (each year looking ahead to 
the future but taking into account what has happened in the 
past which will effect the future) . 

The variables for the problem consist of the number 
of new sensors which must be launched each year, the sensor 
mix which is carried on each spacecraft, and the space- 
craft associated cost each year. In order to achieve a bet- 
ter conception of the potential launches, a "look-ahead" 
period beyond the current year is considered so that the costs 
will .be minimized with respect to the "look-ahead" period 

and .not just with respect to the current year. That is to 

say that it may be economically preferred to launch excess 
sensors in the current year to satisfy later year demands. 

Formally stated as a linear integer program, the 
problem become s : 
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I 


I 


I 


I 


I 


I 


I 


Minimize : 

Total Cost 


L m 

EE (Cost ) ^ jX 
i=l j»l 


( Launch ) . 

^ f3 


Subject to 


and 


(Launch), .is integer valued for i<=l,2...,L; j *=1 , 2 , 


• • m • f 


m 


E 


( Sen) 


k , i 


X 


(Launch). .+(Excess), 

1,3 K 


(Excess) ^ 




Whe re 

i 

j 


k 

L 


m 

n 

Cost. . 

1 . 3 

Launch . 

1 »3 


Sen. 
k , 1 

Excess 

Dem 


k, 


k , j 



for 


i=l , 2 
k = l , 2 


L 
n ; 


index relating to spacecraft type 

index relating to year of the "look-ahead” 
period (current=l) 

index relating to sensor type 

number of spacecraft types 

number of years in the "look-ahead" period 
number of sensor types 

Cost to launch spacecraft i in year j 

Number of spacecraft of type i to launch in 
year j 

i0 if spacecraft i does not carry sensor k 
(1 if spacecraft i carries sensor k 

number of excess sensor k in orbit in year j 
demand for new launches of sensor k in year j 


f 
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As the launch optimization occurs within a Monte 
Carlo simulation an efficient algorithm must be employed. 
Conventional linear integer programming algorithms are too 
time consuming to be considered, and a suboptimal choice 
procedure is employed. The problem is treated as a con- 
tinuous linear program and a simple upwards rounding of any 
non-integer first year launches is employed. Although not 
always optimal, the rounded continuous solution will corre- 
spond to the best integer solution the vast majority of the 
time due to the nature of the setup of this specific problem. 
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5. SATIL 2 PROGRAM DESCRIPTION 


The FORTRAN coding for SATIL 2 follows the flow- 
chart (Section 4.5.) very closely. The following is a de- 
scription of the program implementation. 

5.1 Organizat ion 

The computational variables, i.e., counters and sums, 
are initialized in one of three ways. 

1. data statements 

2 . assignment statements 

3 . the subroutine ALI5T1 

The nominal values of the input variables are established in 
subroutine DFLTIN ♦ The periodic reinitialization of vari- 
ables within the Mohte -Carlo loop is done by subroutine 
ALIST2 . Many of the NAMELIST input variables are arrays 
inputted as vectors in row rather than column order. The 
transposition of those arrays is accomplished by subroutine 
TRAMS P . The reliability equations are implemented in sub- 
routine RELIAB . The linear programming algorithm is imple- 
mented in subroutine NQSAT which calls two assembly language 
routines. The output of the program is done by subroutine 


OUTPUT . 

Those variables dimensioned in the flowchart by the 
number of Monte Carlo runs have been written out to disk on 
standard FORTRAN logical unit 2, and are then read back in 
after the Monte Carlo runs are completed. 
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5 . 2 Control Cards and Input 

Standard FORTRAN compile, llnk-edit, and go cards 
are used if the source cards are being run and standard link- 
edit and go cards are used if object decks are used* A 
U.F4STAT SYSLIB card must be included in the link-edit step 
to locate and link subroutine TRAMSP . Inputs via NAMELIST 
should be entered in the appropriate place in the back of the 


decks. 
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6. SATIL 


PROGRAM LISTING 
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sat: 1.0 2 LAUNCi; Sr^^HLATrON 
2>:TliHf’AL riy.ur 
v.£hl LCE?v,K?V 
INTEGER pFS'JG ,OETT 

FKAL *P F.AT?,!}: 

* 2 I VDEXF^INOiX: 




Dr!Ot!SION i'!ATF,I;«' ( 1n,2.R) (16? /I NDEXC (23) , COZ? (5,23), 

D'^TrNrlON ??P01 (3) ,???02 (3) ,Pr,‘?03 (5) ,?PP04 (5) ,P?P05(5) ,rpP06 (3) 
DI”iS NSTON ? >?07 (3) ,???0o (3) , PPPO") (5) , ??i?10 (5) , ???1 1 (3) ,rp?12 (5) 


DiMRtlSId’ 

F 

P P 1 3 ( 

5) , 

? P P 1 4 

(3) , PPP1 

5 (5) ,??P 16 (5) 

,P? 

P17 (5) ,PPP1d 

DI.*^F 

NSIOK 

? ? 2 1 4 ( 

5) , 

?2?20 

(5) 





IHTE 

SIR 

r 

I 

I.C(3, 

10) 


(10) , 

TV 

(3), lnc:-13 (5, 

25) 


T?;R£ 

g:-:~ 

T 

T* 

{^) , 

nc 

SC (5)^ 

, IPIB 

(5), I3MSX(5) 



IN'"E 

j ' r 

i* 

0 (10, 

25)“ 

r !vCF 

A (5,25 


NSW333 (10, 25) 



I .'I 

Gs:': 

I 

pr^ 

(10) , 

IS!! 

(10 

, p) ,:ic 

SH 

N (10,25) ,MA2>i 

(5) 

, NFAIL (5,20) 

*? IS 

■I I ^ 

T 

I. A 

(5,23 

) 

d? (3 

, 23) , 

IB 

(5,25) 



T 

G>:? 

T 

PI, 

(5,25) , 

lOA (5 

,25) , 

10 

(3,2 5) 



TIME 

GFR 


?L 

^ (5,2 

3) , 

1 ? ’1 R 

(5,25) 

/ 

NOS (5, 25) 



IN'??: 

5FS 


S ( 

3,2 0) 

r 1 3 

S (10, 

5,20) , 

lY 

.:A?.(5,20) 



4EAL 

C(1 

f 

5/ 

2 3 

) , 

:t (3, 

4 5 

) , 

PE (2 5,9,25) 



KSAI, 

r(3 

9 

r>, 

4 0,25 

) ,9 

T (3,4 

0,25} , 

PP 

v'{U0) 



REAL 

PELS 

(1 

0 , 25) 

0 ^ 

9 

S (1C, 

2 5 ) , R 2 

LF- 

(5 ,2 5) , AE- DF (10, 

5) 


REAL 3TA::dI (1C, 3) ,A:-:EA!:L (1C,5) ,A;;:3E (10) , stands (10) ^A«iiMiS(lO) 
REAL LCrnC) , CALC (10) ,??? (3,25) ,AMA:xC (10> 

PEAL ARlRCdO) ,rAi:G?nO) ,UCP (3,1)) ,?2r(3,20) 

PEAL UC (1 C) , A:-:a:*LC (5, 10} , AMI'.?LC(5, 10) ,LCEIC{3,10) 

REAL CAI.LC (3, 1C) , PANGLC (5 , 10) , UCc'LC(5,5, 1C) , SCCST ( 5 , 2 5) 

REAL UCLC (3, 1C) : (25) ,LCST (5,23) ,LC.EL1 (3) , CALS (5) 

PEAL AEAEE (5) ,ArLN5 (5) , RANGES (3) ,i)C?n (5, 3) ,UC?N (5) 

OSr'iRFSLON :.'iO);C(5,25,9) 

IIITiJGEP C =LCAP (5) 

REAL rjCS(.3) , ?SR (d) ( 10) ,?CR02? (10) 

REAL P3R? (10) ,RDS(1C) , J?S( 5) , roR ( 10 ) , Pt\'lG ( 10) 

IMdGK? I?': (5,25) 


HEAL 


d (10) , 

PPLFL (10) , 

p?;:?p (10) , sa 

MC(3,5,2 

51 



REAL 

PCk 

dFd'^) 


, ?LA (23, 3,25) , PPL (25 

,5,25) 



REAL 

so: 

dC(3,5 

,2 3) ,30^'C'7 

(3 ,25) ^SUMSC 

7 (3,25) , 

•3UMu (9 ,2 5) 



RFA.L 

SI). ISO (0,2 

3) , SUM L A (5 

,25) ,SUMSIA( 

5,25) ,SU 

NPI (5,25) 



REAL 

S'J‘ 

^3?L(5, 

23) ,AM>: (3, 

5) ,AMN (3,5) , 

CSC7 (3,5 

) 



" F AI. 

k ' 

■7 (3) ,A 

.d (3) ,C5TI 

7 (3) ,XCS7 (3, 

5 ,25) , MS 

(1,25) 



71 *1 T 

-* »» fix. 

SIDE (3,25) ,MT(3,25) 

,3707 (3,25), 

MLA (3,23 

\ 

t 



HEAL 

STI 

1C 37 (31, 

5 22 5) , A M 13 

X: (5 ,10), 37A 

NP7 (5, 10) , CALLCI (10 

, 5) 


PEAL 

AM' 

lK LT ( 10 

,5) ,.\MA>:L'’ 

(10,5), AM FAN 

7 (5 ,10) 




ir: ?E 

4 F R 

i:a m '' XT 

(3,10) , 003 

'•‘■w (25,1c) ,i 

PTLCI (10 

,5) ,l;ichst (2 5,5) 


REAL 

C T I 

5LA (5,2 

5} , nPL (5,2 

5) ,Si DFL(5,2 5) 




:nte 

GER 

N S F M 0 1 

(2 5) , N5EN0 

2 (25) , N SEN 03 

(25) ,US3 

Id 4 (2 5) ,N3R 

i:05 (2 5) 

IK IF 

GFE 

ys 3N'‘'6 (25) , NS EKO 

7 d5) , NSEKOd 

(25) ,N5E 

NO 9 (2 5 ) , NS ?; 

u 10 (: 

:5) 

T>* 

OFF 

i;:7:-'0i 

(5) 

(3) ,IM7'd J (5 

) ,IMIX04 

(5) ,lMl:v05 { 

5) , I MIX 06 ( 

:) ,10 

I NO' 

7 (3 ) 

i:i;03 (3): ,t;i 

I • C9 (5) ,L'1IX 

10(5) 




s;Hir 

V ?• L i 

e:;cf(im 

1X01 , 

T (1,1) ) , ( I XI 

/, 02, 1 3 MIX! ( 1,2) ) 



REAL 

AhEAN 1 (5) 

, A mean 2 (3) 

,AXFAli3(5) ,A 

mean 4 (5) 

, A M FAN 5(5) 



RrAL 

^ M *, 

CLC 1( 1 9) , A.1XL7 2 (1 0) ,A:'!X1CJ ( 10) , AM:>:LC4 (1 0) , AMXLC5 

(10) 


REAL 

k K ' 

■II.C1 ( L^) , a:-:n'lc 2 (i 

0 ! , A M 'IS LC 3(10), A M N L C 4 ( 1 0 ) , A X N L C 5 ( 1 0) 




I p-LC 1 (d) , IPTLC 

2 (10) , IP--LCJ 

(10) , IPT 

LC4 (10) ,I?T 

LC 5 (10)/ 

REAL 

CALLCI (1'"' 

) , rALI.C 2 (1 

C) ,CALLC3 ( 10) ,CALLC4 

(10) , CAL1C5 

(10) 


TNCb 

GEK 

L NC:iS1 

(2 3) , LlvCHF 

2 (25) , l::cus3 

( 25 ) , L PC 

HF 4 (25) ,L’iCHS5 (: 

>5) 

REAL 

3Ti 

'.N D 1 ( S) 

,3TAND2 (5) 

,37AND3 {:.) ,3 

I AKu4 (5) 

,S7AND5(5) 



REAL 

A “in F I ( 3) 

, A M7FF2 (5) 

, AXT3FK5) ,A 

M7flr'4 (5) 

, A MTQ F5 (5) 



E'.)nr VAL; 

•:1CE (“'1 

7X03, 13 MIX 

7(1,3) ) , (-M7 

XD4, loXi: 

::7(i,4)) 



EO 01 

VAIK 

^NCM'Y 

IR05, JSXRX 

T (1,5) ) , (IMI 

XOb j 1 3 Ml 

XT (1,6)) 




EEPRODUCIBrLITY OF THE 
ORIGiNAL PAGE IS POOR 


SATL0100 


SAIL0200 

SATL0300 

Sn LLCSOO 
5ATL0600 
SATL0700 
SATLOoOO 
SATL0500 

S/ITLIIOO 
SATI.1200 
SAIL 1300 
OATL14U0 
SATL1500 
SAIL 1600 
SATL1700 
SATLIbOO 
3 A T L 1 9 0 0 


SAIL 22 00 

SATL2500 

SATL2600 

3ATL2700 


SATL3000 
SATLdOO 
SAj:L 3200 
SATL3300 
SAIL 340 0 
S A I L 3 5 0 0 
£S Ai’LioOO 
SML370U 
SATLJfaOO 
SAIL3900 
SATL4000 
SATL4V00 
SATL4200 
SAIL 4 3.00 
SAIL 44 00 
SAIL 4 500 
SAIL 4 600 
SATL4700 
SAIL 4 a 00 
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VALENCE (iJCA ND"", STANDI) » (ANAMLC, AHAXLT) # (A .1 EANL, AflEANT) 
r.on: VALENCE (^sr- k.ise:/.:) , (lhcis.lncust) 

Kvjui vALErcFC'*:::''?, isr.ixT (i , 7 ) ) , cni:: 03 , isntxt ( i, 8 ) ) 

z^ur valence r •!:x''9,:3mi>:t (i, }) ) , isiax*: (i. iC) ) 

E.)UivAL-.':c2 (rrr:J0i,NO3E:J-^ (i,i)) r risEN02 , :jcss::c (1,2)) 

(NFEI'D NOSErC ( 1 , 3 )) 

(ESEI.C^.NCiCN'^ ( 1 / 5 )) 

(N^rNOT^NOEENT ( 1 , 7 )) 

(NSENO^,NOSEr’ ( 1 , 9 )) 

OZEAN 1 , A 'EAN" ( 1 , 1 )) 


zyUE vale:: 

FQUI VALEVCE 

eOu: vALn::c*" 

equi VAI E"CE 
sy u: vALzi’ct 


EO'JI VALENCE 

(a::fa;J 3 , 

A A LAN I 

ecu: VALE'' 

V. i. 

(A;-*L AN 5 , 

A.MEAN' 

EO'JX VALuN 

Cr 

(A.vnci, 

A V AXL' 

tear. VALJ:: 

CE 

(A Cl, 

A.UXL' 

EQUI VALEN 

CE 

(a:*.nlc 5 . 

A NAZI 

Z 3 U: VAIE*’CE 

(A’^NLCI, 


EQU: VAL'CN 

C - 

(AMNLCl, 

AMPL' 

?Q UI VALE’-' 

C- 

(A'*M C 5 , 

a:: IN I' 

EOUi VALor: 

r" £ 

(A.NINI;, 


ECUI VAL"’I 

CE 

(AMCIi "T, A 

m:df) 

EQU: VALEN 

cz 

(NTEZN, NCSZHT) 

F.'iUI VALE'* 

CE 

cruel , 

ircLC 

EQU: VALKV 

cz 

(TPTLC 3 , 

I 2 TLC 

FC'H VALr’: 

cz 

(TPTIC 5 , 

:?TLC 

ZQ'JT VALEN 

CE 

CPTLC,: 

P’LCT 

EC uc vale:-' 

CE 

(CALLCl, 

CALLC 

ecu: valcn 

CE 

(CALLC 3 , 

C A. LLC 

zou: VALt’N 

' C E 

(CAa- LC 5 , 

CALLC 

Egu: VALS -1 

ICE 

(ca:lc,calic: 


(S 3 CN 0 9 , NCSE?-: ( 1 , 9 )) 
(VS3N36,.;OShNI (1 ,b) ) 
(:i 3 Z!iOf 1 ,NOSEliT ( 1 , 6 )) 
(NS£Nio,NosE:r: (i , i 0 ) ) 
(A;isa'! 2 , A.iEA:;i (i , 2 ) ) 

(a;isa 3 4 , a:-.ean: (i , 9 ) ) 


1 ) ) , (A-.XLC 2 ,A« AyLT ( 1 , 2 )) 
. 1 ) ) , (a;u:c 9 ,anax-t (i , «) ) 
5 )) 

1)) , (A.NNLC2,h*'. 2S1C (1,2)) 
-(1,3)) , (A;iNLC4,ArliriC (1,9)) 


(:?CLC 2 ,IPTLCC ( 1 , 2 ) ) 
(:?rLC 9 ,i?iLCT ( 1 , 9 )) 


( 1 , 1 )), (CALLC 2 ,CALLCr ( 1 , 2 )) 
(1 , 3 )) , (CALLC 9 ,CALLCC ( 1 , 9 )) 


EQMT VAI 5*’C£ 
t’QUr VAICN Ci 
Z<.)UT V .MZr Cc 
EOU; VALfr C7. 

valence 
eod: v.me::cz 

EC ’J: VALENCE 
EQl’I VALZ::CE . 

Fy'r vALJ::cE (AiiC jee, a;iT 3?': (i ,5) ) 

rcuc vr.L-::;cF. (?3t> ( i , i) , ???oi ( i) ) 

£Qu: VAL-Ncr (??r ( 1 ,'i2) ,pp:72(1) ) 
E?UI VALt'i'C? (P^='(1,03) ,P??03(1)) 
ecu: VALE!: CE(?P?( 1,04) ,??P0U(1) ) 
EuU: Va1E?'C6 (F?? (1, 05) , PPP05 ( 1) ) 

KQDI VALr;::Ct (PP? ( 1 ,Oh) , PPPOo ( 1 ) ) 

equi v?.le::ce (ppl ( i ,07) ,p??07 ( i) ) 

EUUI VALENCE (PFP ( 1 ,09) , PPPOP ( 1) ) 
EQU’' valence (’?^ ( 1,09) ,???.V»(1) ) 
VALE!' Ct’ (PPi (1,10) , El'F 10(1) ) 
RQUT VALENCE (P?r ( 1,11) ,?PP1 1 (1) ) 
EO'J: VALENCE (PPP( 1,1 2 ) ,PP 0 12 ( 1 ) ) 
ecu: VALr.!‘CK (?PP ( 1 , 13) ,PPri3(l) ) 
sou: v.MLNCE (f'p: ( 1 , 19 ) ,r??i 9 (i) ) 
K 3 UI VALENCE (PPIM 1 , 15) ,?r?15(1) ) 
ecu: vALE"ci; (p?r{ i, i6) ,?prifa(i) ) 
EOUT V/.Lc,::CE (?^?(1,17) ,???17(1)) 
ECU' valenc'p (PP? ( 1 , 1'’) / PPr 10(1) ) 
S'CUI VALEI.CL (por (1,10) ,PP? 1 9 ( 1) ) 
ecu: VALEIICE (I'PIM 1,20) ,PPP20(1) ) 


(l:iche2,lncmc: (i , 2) ) 

(LMCHS9,LeCH3I (1,9) ) 


(L.\c!-si,: NCI'S : (1 ,i)) 

(L:;c:iS3,LNCH5: (1 ,3)) 

(LNC!.S5,1':CHET (1,5)) 

(^CANDl ,nTA:iDT (1,1)) 

(SCA ND3, oIANU" (1,3)) 

(3T.V'i)".,STANDT (1,5)) 

(A-: 11 -'1 , A'^Cli (1,1) ) , (AETB’2, A.1T3:T (1,1)) 
(; -!Tl-.E3, AMTPE" (1,3)), ( AM rB7 9 , AMCbE ^ ( 1 , 9) ) 
(AiiT J, 


( 3 CA‘’D 2 ,SrANDC ( 1 , 2 ) ) 
(S 7 a::d 9 ,sThi;i)t (i , 9 ) ) 


SATL 9900 
3 ATL 3000 
SATL 5100 
SATL 5200 
SATL 5300 
SACL 5900 
SATL 5500 
SATL 5600 
3 ATL 5700 
SATLSbOO 
SAIL 5900 
SATL 6000 
SAIL 6100 
SATL 620 J 
SATLbjOO 
SATL 6900 
S ATLO 500 
SATLooO'J 
SATLO 700 
MBS 

SATL 6300 
SATL 6900 
SATL 7000 
SATL 7 1 00 
SATL 7200 
SATL 7300 
SATL 79 C 0 
SACL 7500 
SAJLToOO 
SATL 7700 
3 ATL 7 O 00 
SACL 7900 
SATL 6000 
SATL 5100 
SAIL 820 0 
SATLd 300 
SACL 8900 
J A 1 L 350 0 
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CO «• T 5/ !1I , F E , X A y rw 1 r , SI D ii , A « T , !l 0 1!.'T , C5 Tir , ? , T , ST DT , KC 5 : , 

1 5TDCS?, PLA,IlA:-:N,:'.Lh,STDLA, FVS, ?VI, P?V,K?V,ST0?V,IT1,DETC 

2 ,P7,PPL, «FL,STCVL,ArK,CST:, KAXLSC, OS 
c'*.x'*ON/?r :.'T/D*;ruii 

‘ cos>'*ov'^o’’T/:iA, :n,:o, ib?.,:os,iP'i5,i?LE 
cor-o’:/L:::fc.\r./cc2"’,Li,!.?v;,sj,*.MSs:x,:4E'«sEN 
C0S'«0 'J/ImU''CH/K£AT( 5) , IC05I 
CC O’j/Ol S K/C , CT , A N ? V 
sou: VAlEr*C*E(IHO!IT (1, 1 , 1) ,ILA (1 ,1) ) 


AL?{ '() »-1 . ♦ (MOG V' (••:) -1.69?)/. 301 

0EV(x,y) =5osr(r/A:'AS?-y**2) 


DATA UCLC , 1 C ^LCV-Ci / I 1 ?»C ./ 

DATA 'lOSA , ••03/23C* '/ 

DATA ITS, i/.p,:;:/ 15*0, 

DATA '•A>'I/25/ 

DATA Ar.~3 r/50*5./ 

DA~A :rDS‘*:r, t'T'i".i/Ji4*c/ 

DATA fATr.IM/a'46-0. TOO/ 

data P?rO 1/. 5, .23, . 15, .07, .03/ 
DMA ?’?0 2/. 3, . 25, . 2, . 15, . 1/ 
DATA ?r?0 3/. 3,.3,.2,. n,.07/ 
DMA PprO'-i/. 15,. u, .15, .07, .0 3/ 
DATA ^PF0 5/.2 1,. .12 , . 2? , . 1 5 , . 0 5/ 
DA-^A P7P05/. 23,. 1,. 23 , .16, .05/ 
DAT.A P??'> 7/. 25,. 35,. 25,. 1 , .05/ 
DAT* PF?0 •]/. 1 0, . 4 0, .24 ,. 0 } ,. 02/ 

DA“1 r??0y/. 12,. .12, .3 2,.!*’,. 37/ 

DA-"*, ?r?VV. 13,. 34,. .17 ,. 12,. 0 2/ 
DATA ???1 1/.2,.2,.2, .2,.2/ 

DATA ???12/. 15,. 22,.25,.22,. 15/ 
Data r?pi 3/. i , .2 5, . 3 , . 25, . 1 / 
D.».TA “??1 4/. '' = , . 25, . ?4,. 2 5,. 08/ 
DATA ?PL‘1 5/. OS,. 25, . 4 , .25, .05/ 
DAT'A PPP15/. 1,.2,.4, . 2,. 1/ 

DATA 7P?1 7/.C.1,. 3, . 14 , ..1, .03/ 
DATA PPP1 3/. 05,. 2 ,. 5, . 2, . 05/ 
O.ATA V. 7 3,.:,. 5 4 ,. 2, . 03/ 

DA*1 P?r'2 0/.0 3,. 07, . 5 , .07 , .0 3/ 


NA:*- LTST 
2 f?Dl-.IT, 
?:5E*:06 
IKIX'^7, 
A«TR?5, 
A;1 FAN’ 4 
ATT'JR, : 
A« -aC4 
I?"LC3 
I. NT !IS2, 
PDKBf , 
? ?P01 
2 ? ? 1 1 


/ : '’ ? j '• / .A ' N- , i A ;• 1 s c , :i ? r L s , .1 A X n , .t a :c l n , :: a :c i £ , a 2 , :i ax:, 

!: 'I’ E T • 0 1 , »• 5 E *: } 2 , N' : S C 3 , N 3 E 0 4 , N S 5 ;j 0 5 , 5 £ 11 0 o , r S t X 0 7 , 

, MS£*:o >, MSTin o,:mixoi , 1:1: ~02 , :/.ixo3,:r.i;:04,::’:AOb,::iii0b , 
I :: r •: 0 i , : 1 1 /. 0 9 , 1 'i I : ■ 1 0 , 1 p M : , a M T s. ? 1 , : p £ 2 , a n t k f 3 , a .1 t 5 f 4 , 

S r A * ' D 1 , 5 T A N D 2 , 3 T A : : D 3 , 3 : A : D 4 , £ T A X 0 5 , A '1 : A 1; 1 , A M : A 1 1' 2 , A P. £ A 3 , 

, /tTEAVS, Ar.r33 , 3TA”D3 , AyEA.;.’?, A PA/.C , AP I’’ C ,I ?I , C AlC , APAXi , 

? ” u , C A ! P , .A .1 .1 ; : 1 , A p : N P , I F : P , C A !. 3 , A V A L C 1 , A .4 C 2 , A P X L C 3 , 

, AMir:, rp‘i: C1 ,AP'?1C2, .\:IMLC3 , .AP‘:iC4,A;;NL25,Ii>iL21 ,:9:iC2 , 
, IFTLTU, • ?TLC5,CALLC1 , TALL02 , CAL LC .1 , Cn L LC 4 , C«LLC 5 , 111 CHS 1 , 
Li:CHSJ,LVCH.5 4, LXcnS5,Oa?CAP, ?FS, F5S,t>b ?.3,?bSF, POP ,2050?, 

? D r. , ? *-0 , ? ? M , r ?y ' p , r ? Ls ? , p A < p , II N , D • , D F rc , 

, Pl-P0.3,r?P0 l,??PC4,p?r05, ??Pl6,Pr'?C7,???09,PF?05 ,???io, 

, ?P?12,.3??1 3, ?.>Pl4,?.^e 15, ?P?16, ?F?17,???1 o, ???19, PPP20 


3 , D E2 U 3 


NOb 

.IDS 


SATLOwOO 
3 aTL'95»i< 0 
SATL9o00 
3 AIL i» 700 
3AIL9bwO 
SAxLOjOO 
.lATLOOJO 
SATLOICO 
SkZlO^ 06 
SAILOiOO 


initialize input data to DiFAULT VALUES 


SAIL0500 

S.AlLObOO 
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.11*0 

!iO:*IT»2S 
DO 5910 K*1,3 
DO 5 913 '1*1,25 

C7 (:< 

s'lic r (K,r) -* 0 . 

5913 SOIS c: (<, V) *0 . 

Aixr (K)*-993. 

A.1T( K) «1'». KIO 
DO 5 9 IP LSC»1,5 

5910 Ar.'-’C K,LaC) *1 0. n P 
pV7»-9Q}. 
pv:;* 10. El 0 
D.^D'JGoO 
;;3*3 

CALL D7L^:i:( P?^P, 

1 V 0 s I K , : r : r , : J I c , -1 7 E E , r A :i 3 L , A 1 £ Ai: L , A.1 1 B 3 , s : A ;; j s , A is A !.’ s , A :i A X c , 

2 Am 7 c , I 3 T , C MC , A « A •• B , A : N 3 , 1 = r 3 , C.\ L B , A X .1 , A LL' M , I i 7 1: , C H I. .=1 , 

3 A '. A.*LC, A/.::;iC, :p:LC , CALLC, L NCHS , os JCA? , PPS ,? P.,, PBP.S , PuPP , P G3 , 

3 ?.9 POf , PC SB 7, pen, PP.IC, PPrfl-.?, PPL."?."? ,1 AXI , U:; , DP , L£IC) 

c 

C TKPijr 

f::ad (5,::; POT) 

•.'^IT h (6,If:?rj71 

W?.ir F (o,70C) 

■•00 ?0FrA’(lHl) 

W?I:e (A,32«) 

329 •'0?’'. A~ (1C-* , ' 3p?' ) 

DO 3 2<> ’=1,2*> 

32^ W-:-E(6,32T) (??? (IT,!) ,11,1^5) 

327 .^CP.IAT (3;':,5r7.2) 

W^IT r. (6,7 00) 

C VALTD.'.TF 

^ * <t * * * * w* <t A » fc), >» » * A 

CAiL 7 ':a 

CA LL 7.-:.;?r F ? ( a • "B ."T , S , l O , iT Aa 4 ) 

CALL "'?A::£?(L::CH3T,25,5,CTAE4) 

CAi:. 7 ?.A;;s?( 3:.V:3^,.S, lO,37A.7a) 

CA LL TAf S P ( CA LLCT , r: , 5 , 3 1 AF 4) 

CAIL :^AN'.3?( : >71 :t, 10, 5. 3:aF.4) 

CALL 7F.A^3?(\>;:i:L'", 1C, 5,5TAF4) 

CALI. T?.AN’SP(A.'’.a;{LT, n,5,STA?.4) 

CALL TP.A;:S?(A.-.CA^T,5, 1C,?"'A2U) 

CALL A'!'?(:s:::;--,3, io, .9CA.-.4) 

CALL r.PAjrs? ('r*.9r:;c,25 , lo ,3TA74) 

C 

«?.T7E (•',3 4 3) 

34 3 B‘Or'tA7( • I9 mi:'<) 

344 1=1,10 

34 4 V7 I- V (6 , 3 (I 5) (I S.'* IX (I , J) , J = 1 , 5) 

34 5 .=‘o?.r; A'’ (3::,5:a) 
c Lizn a.*? tt;]- 7 ,i 

w 

, ... ’^S PI , P , PPL , re , ?p V , an- p V, SU.IP V , s u r.c p V , 

, /?PD,3U.1LA,SUK.SLA,rLA, 

4. ,3J*'v-,5U1.5C,i'L:\,JUXJ,SUK3J,PC) 

c 


I 
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at 


<w> 


f 


i 


I 


* 


% 


» 


> 


I 


COMPUTATION OF f;US SFIIALULITY r, ^?0B. 07 PAIL UTF/YEAF. 

CALL PZLIAE (MA:*Io, 

1 «AXVB,VAyv^Pj,y^:.jj;,3^^^^^,^STANDL,A.\SA:;L,-AlL5,p73,F.LLS,AKTBS, 

2 S*^ AND?, AMEA'lf) ' 

•fr IT E (b, me) 

34 6 rOF''AT(/' TFh') 

DO 346 1*1,5 

346 4MTE(6,3U7) {?7U (I , J) , J«1 , 10) 

34*» FOM AT (10/, 10.-( . 3) 

• ? II £ (6, 3 4 ^ ) 

34 9 ?0 ?',:.*(/• ?F7') 

DO 350 T«i,iO 

350 W/IIF (6, 347) (I,J) ,J»1, 1C) 

I® (0 -DUO. LC.C) GC TO 3 51 

DO 3 52 1* 1,25 
DO 353 J*1,5 

353 PyQ(J,I)*0. 
no 3 5 4 .1* 1,10 

354 P?.5( J,T)*0. 

352 CO’;: 

351 COr.TIN'Jc 


SATI.3400 

SAIL3500 

SAIL 3 7 00 
SATLJeOO 


C J.EAr.NIXJ CUrVE EaPCNEMIS FCF SI'ISCPE 

C 

DO «^02C IF O , P?-.;*LS 
e020 LCE(L?)*Al? (CALC (LS) ) 

LCEF -u? (c.nr-) 

DC 3 021 %n=i ,:i.v-:.;3 

902 1 LC2J (VU) =AL? (CAL 3 (N3) )-* 

DO <5 022 LFc* 1 ,;-;a.-:lsc 
DO 6 02 3 I 5*1, MAX IS 

9023 N^SSC (LFC) *';OF5C (irC) ♦ ISKi:-; ( LS , LSC) 

9 02 2 CONTINUE' 


c 

DO 10 29 LS*1 , .ma;'ls 
1029 irj(LF.)*0 
ITL*0 

DO 1 0 30 r = 1,«A:.T 
DO 1031 L£ = 1 , r./i'-CLF 
IF C;OSE”(LS, 7) .L t. O) GOTO 1031 
I' (riL.^’u. 0) 

iF(::s(LS) .Eg.C) :rs(is)=N 

1031 con:i*;u" 

10 30 c.' NT IV ns 

u 

C COMPiJrS PISCOU’iT FACTOR 

r* 

7A>'V ’rvz:KAvN>!H 
DO 1 040 !.*= 1 , *<A V?; s'N 
104 0 ?Vr(N) =1.0/(1.0*. ''1*DF) **(%•- 1) 

C 

- START 0? ACTIVITY FOa SPAC5C3AFT AND DUSES 

DO 1 04 4 L0C*1,.TA\'LSC 

I>: (LSC) =-5 

DO 1 045 L<5 = 1 ,.'1AvlS 


SATLJ90C 


SATL4300 

3ATL440C 

oA4A44I>00 

SA ri4bOC 
SATL4700 
SAXi-**sOO 
SATL4900 
SATL5OJ0 
SAIL5500 
SATL56D0 
SATL5700 
SATL6C00 
SATL6 1 Ou 


SATL7000 

SAILTIOO 

SATL7200 

SATL7300 


* 


t 



) 


.1 


» 







I 


I 


I 


» 



it 


I 


I 


» 


- 

... 
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— 

*IF.AD»ITS(IS) ♦TSf'TX(LS, LSC) 

-- 

SATL7400 


If (:r-\D.G*.i '(isc)):>: (lsc) *:?.a3 


SATL7500 


cr-\T IN -It; 


SATL7900 

104U 

CON? I ‘it’?. 


SATLooOO 


00 10u7 V^=1 , rAXNB 


SA7L9100 

1047 

i*:o{ i:n) *i'l 


SATLc400 


C 

c 


C 

3 

SAIL3500 

0 • • • • 4 

LAUNCH COS? LEARNING EXPONZN? 

SAILSoOO 


00 1C50 :l‘ci ^.nAxiz 

SATL5700 


DO 1 0 50 LSC= 1 ,r.A’*LSC 

SATLBBOO 

1050 

C 

LCLLC(LSC,:S) =ALP (CALIC (L^C, IS) ) 
W?.I7 r (6,3 55) LCL 


3S5 

•:0F/.\?(/* L.'--' , 10E12.4) 
£ (6,356) LCI?. 


356 

?<^R.‘'AT(/' LCIr=* , 1i12. 4) 
WRIT t (6,357) LCSD 


357 

.?K A? (/ ' IC ih = • , 53 1 2 . 4) 
W? i: L (6, J5H) 


.359 

70F1A’'(/' LCilC) 

-’Rir £(6,.H7) ((LC?LC(T,J) ,1=1,5) ,J=1, 10) 


359 

FCR I A? (3X, 10512. 4) 
Wr.I’” E (6,360) TOL 


36.0 

roR:iAT(/' i:i=',ii4) 
£ (6, 361) I".5 


361 

70?.«\'!(/' IT 7= 1014) 

XRT7 £ (6,36 2) I" 


362 

70r.5A"(/' 1.5 = ', 514) 

W^rCE (6,36 3) lyo 


363 

30H‘'AT(/' :••:?= ',514) 

w: :,’F (6,364) ( JVT (:;) ,r = 1, 10) 


364 

7074;.<7(/i rvf=» , 10312 .4) 


w • • • • 4 

, . .-D 4 

DO 1 059 13=1 , ;'.Ay LS 

I? (.A ;5.AXC(L3) . 50. A.1TN: (L3) ) GOTO 1059 

SATL9100 


pan; E (17) = (A'* a:: c(is)-a.i:‘' c(is ))*.2 

SATL9300 


I? (I (IS) .Lr:.2C) GO^.l 1051 

SATL3400 


07 1 062 1=1,5 

SATL?500 

1052 

!tc?( I, LS) =PP? (e-:,lPT (L3) -20) 

SATL9o00 


GOTO 1053 

SATL9700 

1051 

DO 1 0 54 1 = 1,5 

Sat 19600 

1C 54 

ucp( I, L.*;) =???(?,:?? (IS) ) 

SATL590O 

1053 

con: IN'tF. 

SA'LCOOO 

1059 

c 

CONTINUE 


c 

i . • . • i« 3 A 

3ATL0700 

c 

4A 

D.O 3 5 NU=1,NA'*N3 

IF (Ay.A.>:B(ni) . FQ. A.'IIND (NO) ) GO TO 38 

TAf:;En(NP) = (a.v\aU(nd) -amsb(nd)) *.2 

SAIL0900 

c 

. ...SE*" up P’^OS.A 3I1ITY DI3TN 0? E'JS COSTS 

SATLUOO 


IF (I PTl’ (M.) . LG.2 -7) 3010 noi 

SATL1300 


DO 1 062 ■^=1,5 

SAIL1400 

106? 

NCPP (",M9) =•??? (6-1, TP-'B (N 3) -20) 

SATL1500 


GOTO 1063 

SATLloOO 

1061 

3<'’ 10 64 :=1,5 

SAT11700 
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1C64 

•jcrj (:,N'B) (kh) > 

SAIL 1800 

1063 

CON*: INtIM 

SAIL 1900 

18 

C 

COKI 

SATL2J00 

c 

c 

p 

AT 3 D 



I? (AKAX:?. FO. AMI':/) IDTO 39 

SATL2700 


RA NG SK= (A A>: .1-A INII) • . 2 

SAIL2CG0 

%* • • • • 

...31? UP DISTN. OF 3/C •'.AIN?. 5 aSPAIP. C0ST3 

5ATL29O0 

C • • • • 

. . . ( AS 1^ C • S/C COST) 

S A i. L 3 C 0 0 


IF (I ?TN.L£.20) GO-’O 1070 

SAIL 3 100 


DO 1071 1x1,5 

SArL3200 

1071 

(T) =P??(6-I, l»'-a-20) 

SAIL3300 


10"*? 1072 

SAIL34CC 

1070 

no 1 o*^'4 1, 5 

SAIIJSSO 

10”»4 

ucP'-' r:)*2??(i,i??:i) 

SAIL360: 

1072 

n 

cox: i.'-uF. 

SAIL 37 OC 


A? 3C 


c 



39 

00 34 12* 1 ,«AXTR 



.3 6 LSC = 1 , AI L sc 



IF (A •'.A'CLC (is:, ::.) . 20 . .A MlN LC (ISC, I 2 ) ) GC70 36 

SAIL460C 


SA*nLC(LSC,:r) = (M1AVLC (LS:,I£) -AXINLC (13C,I2) ) * .2 

3AIL470C 

C • • • • 


S AT 14 9^0 


IF (•" P'^lCdSC, T7) .12. 20) 30T0 1080 

SATL49CC 


DO 1 OH 1 1*1,5 

3 AIL 5000 

1081 

UC?LC C,LSC, r^) = ???(6-I,:?Ti:(LSC, IE) -20) 

3ATL510C 


;oT> 10=^2 

SAII5200 

1080 

DC 109 3 1x1,5 

SATLSiOO 

1033 

•j:r: c (i,i.3c, is) = ???(i ,i?TLC(LSc,ii)) 

S AILSuO j 

1092 

CON^ IMIE 

3ATL5500 

36 

CO *r" I •III 2 


3U 

c 

co%'7i:in£ 



KF *1= ■'» >:‘i+ «!| 

SA7L6500 


00 10P8 NN 

3 AIL bo 00 


DO 103d is = i , 

SAIL6900 

10«3 

::o r>" x ( is , n ) = :i o 3 e (is , x am « ) 

3AIL70CO 


MCI' E (6,37?) 


377 

- JCv.-: .a: (/• •10.5?::=') 



WPi: 2(6,3 78) ((’tOSS'Id ,j) ,J = 1,2 5) ,1 = 1, 10) 


37 8 

FOr:VA7(3\',25I3) 

. 






c 


c 

c 

C nri;?" 07 M0'^?7 CAP.LO 100? 

c 

I C 5 ■? 

C "D 5 

DO <? 0C1 '!?=1, 

C ir’TTI.Mr?.£ AC ?£?. LIST »2 


SAIL? 100 
SAXL7200 
SATL7J00 
SAIL7U0C 
3ATL7O00 
SAIL770U 

3ATL7bOO 


CALL ALIS72{ C02T, 

1 '*A>: LSC,''.AXL3,.“AX' ,r^, 3CC3T, LCST, ILA , 13 , IB? , Il'L ,10 , ICR ,I I PL3, 
^ NFAIL,I£,Ii?S) 




ft 
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li 

' # 

I 

I 


r *******•**•*••*** 4.* ***.»*Aii* 

DO 5010 K=1, 3 
DD 0010 "=1,25 
CT (K , *J) *0 . 
on 0010 L5C=1 
9010 ?(K, :oo,:.) x^. 

Ar;pv=c. 

c 

DO 10 '3 1 l;?=i,.<axls 

ir(A?;.\::c(L5) a.i:-ic (ls) ) jotc ioo2 

oc (L ?) =A-. nxc (L.') 

10 ro 10 01 

1092 ACCxC-.n 

^"=0 NDY5 ( Df'M.IV) 

DO 1093 ”=1,5 
ACC= AOr + 't'"? ( I , LS ) 

(A CC.LI . FV) OOrO 1**93 


?!iV30T COo” CO'iO'JTATIO: 


*^*''"-<">*<--"*(SS-ACC.acP(I,It)).tA!.-G 2 (ls, 

iOT“ 1091 
100 3 COTI'J'JE 
10 01 corri'jnj 

C '£ 9 

DO 2 0^1 *’ii = i 

ZT (A !* AXu ('.iH) . '’E. A’' :N*n (XB) ) GOTO 2002 
nos( rp) =A '‘Avj (VB) 
i10~0 2C01 

2002 AC 3=0.0 

R.i-?. NDYS ( DtJM vvj 
DO 2003 1 = 1,5* 

A :c= ACC«'”C ?P (T 0 ) 

1’=' (ACC.L". v;CT'^ 20*0} 

v'a?'r;’::rar'“‘'''’ < '3) • c- U M^s-accucpd »d, , n.A.N'3cii (n3) 

SOTO 2001 

2003 CDNOrN'fJp 

2001 CON":;: OB 

r 

c 

p ’^^^fofoishmsnt and kai:;t. cost CO.^PHTA'^ON 

c 

I? (A:VA7'1. ME. a:'.IV.'1) ooTO 2006' 

UO'*= A^AZ/; 

GC'"‘0 2 009 

2006 ACC=0.0 
i?y=:<NDY5 (DrjvMY) 

DO 2 00'^ Tat 1^ 5 
ACC= ACC + fJC?r (I) 

IF (A cc. r.T . h:j) goto 2dd7 

) *aA«=E,VUC?.". (I, 

2007 co>n:iK:JE 
2009 COlJ.'irUB 


10 
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DO 2 020 

DO 2 021 L 3 C» 1 ,KAy.ISC 

T 7 (.\«A:ac (V^.zry) rtii::LC(LSC,: 2 )) jcto 2022 

flCLC (I?C,IF) aC(LSC,It) 

;o ro 202 1 


2022 ACC= 0 . 

2 K=r iiDY 5 (Dl'!“.?tY) 

D? 2023 I»l ,5 

ACC» r.CC + uCPLC {I, L 3 C, Z£) 

T? (ACC. 3 K.HN) 30 70 202 U 


2023 C0trT!.”JE 
^ 0 ”* 4 ^0 * 1 ^ 7 N U i 

03L: {L5C,I2) =A:a::LC(LSC,lE)^?A::3LC(LSCME) *?LDA7(I-1) 

2 ♦ (“V -Accrue n.C(I»LbCr72)) *?AN3LC(L37/ li) /'JC ?Lt {1 , lSZ, It) 


2021 COVTZ'rJZ 
2020 co;rr:rus 


c ^ n 

S/C COS'" IS’ Y£A'» (3US<-3S:JS0R) 

DO 203^ LSC=1,yAXLSC 
NQ=r. ;( (Lsc) 

DO 2 031 >=:iQ, 

DO 2 0 12 IS = 1 .XAy.LS 

!>' (I TS (L*5) . '.I. S) 50 ro 2032 

y=i!- ir.1 (IS) >1 

f:jx“ snx+yLCA'O (r.F;ii;c (ls »lsc) ) *uc (l 3) ♦?»*lcb (ls) 

203 2 COy^ISMt 

17(171.07.':) 30 70 203 1 

7=1- :”i.+ i 

3CC37 (L5C ,::) =oc5 (ij:ir.y. (LFC) ) * {? ) **izz 3 (Ib.'iIa (lsc) ) ♦su^ 

203 1 coNiryot 

203 0 co.iriyuc 

c 

MK !J= 'lA;'?!* !**! 

DO 2 03 3 •; = yAX':,i:';*; 

DO 2 033 LSC= 1 , KAyiFC 
2033 SCC37 (L3C,K) -sees’" (LSC, MA::H) 


DO 2 034 LSC=1 ,KAy.LSC 
;IQ=TX (L3C) 

DO 2 03 5 :: = NQ,xAr*: 

ir (! l<~!IS(LSC,n .MS.LIICMS (LSC,::3)) 30 re 2036 
-IX (LSC) ♦ 1 
GO ro 20 3 7 


2036 FM =1 

i:c(Lsc)=y 

20 37 ICS" (LSC, N') ='JCI.C (L3C, LSCHS (L3C ,'.J) ) *? 'A** ( LCELC ( LSC , LMCSIS ( LSC , N) ) ) 
203 5 CONri’^'JI 
20 34 co:m!niF, 

DO 2 0 3ft s' = .':.\ys’,::NN 
DO 2 033 LKC =1 ,fiA;-:LSC 
203 ft Lcn^ (LSC, y) -LCrr (LSC, ’'AXO 
DO 3 00.1 ’IS’:; 


pvrx =pv? (N) 

DO 3 00 ft LSC» 1 ,XAXL 3 C 

90 Cft cor.r (I.FC, *.) =pvr; * (sees? (lsc, n) ♦lcoI(Lsc, ::) ) 


non •! n n 
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Ir (r. A7F.VE. 1) TO 5900 
WPTT2 (6,365) 

:0:.’UT(/' SCC5*?') 

DO 366 1=1,5 

V5r?5: (6,357) (SCC-T(T,J) ,J*1, 10) 
A* (31 , 1075. 1 ) 

WR-TTr: (6,3 59) 

?0F:«AT(/' LC3T») 

DD 3 59 := 1,5 

•’=ir e (6,357) (ICST (I,J) ,J=1 ,10) 
yR'T F (6,3‘»0) 

?0 9iA-(/» cozr») 

DO 3 71 :» 1 ,5 

wril K (6,357) (CCZT (T, J) ,0=1,10) 
WRT: 2 (6,372) ’’CD 
FOn?.*^(/' UC3=‘ , 2'',576. 1) 

2(5,37 3) MCK 
format {/' UCM=' , 1?5. 1) 

WRIT 3 (6,374) 

FO:-.1A'’(/' UCLC) 

DO 3 •’5 :=1,5 

WS:-^E (o,3-’5) (rJCLC (I , J ) , 0 = 1 , 10) 
rornAT (5>;, lOfS.i) 


TIME T.OOP 


.... "a 12 

DC 9''C2 ':=7’'L,::a ::n 

K 1 =•; ♦ 1 

if (•: .'’F.ITL) 30 TO 206 1 


.. CO/.P’JTE NCPA 
L3C= •'A”LSC 
LS=;i A'^LS 
Xr'.AT =LS*N3*1 
nma:’= (I3C + L0) *::3*i 
DO 9 313 : 10= 1 , 10 
DO 9 312 125=1 ,25 

N'HV;? S’! (T10 , 125) ='!0.?i:N’ (110,125) -MOC?(I 13 , N) 

Ci E;;:12‘- (II 3,12 5) .LT . 0) 'JEW SEN (110, 125) =0 
C3N^ I Nil 2 

DC 9 314 12 5= VI ,25 
.T25= 25 + N1-I25 

NEWSSN (11 3,.I25) =-!!E«SiN (1 1 0 ,J 2 5 ) - N EW32 N (I 10,J25-1) 
IF (’• FWSEV (I10,J25) .L". 0) 'iFWSEN (I 10, J 25) =0 
CO I N l]Y 
CO NT IN7F 

CAL!. NOSA: (.VAT7II,INDFXP.,INDSXC,.1.VAT, NMAT) 

I? (! AXF.N g. 1) 30 TO 5902 

y? ir E (5,3 32) N 
70,'^VA"(//* N* Sl"*^) 

W'^:' F. (5,3 34) (NrC?(I,‘») ,1=1,10) 

70=..': AT ( ' NC 3?- ' , 1013) 
yt c (6, 3 29) 

FOF:". ATC .'.’hMSEN') 

DO 3 .30 := 1 , 10 

21’ IT E (6, .3 11) (NEWSEN (I , J) , J = 1 ,25) 
f 0?V A'’' (3**:, 2513) 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IB POOR 


WF IT E (6,3 3 3) KSAT 
333 F03:i.\T(' KSAr= ‘,5:3) 
5902 DO 9 0''5 LSC» 1 ,.rAXL.'C 
?JOS.\ (LSC, V) =KSAT (ISC) 
90'*5 CCNTIV'Jr 

c 

w 

GO TC 13 
2061 CONTINUE 


DC 

2 062 L 

SC»1 

f •»£% 

XLSC 



V M 

=1 AX7. (L 

3C) 





DO 

2 06 3 - 

= 1,1 

HI 




IP 

(IS (ISC 

,^:) . 

r* 

' M t 

0) G.** 

TO 

206 3 

17 

(:<CSA (L 

sc, N 

-1) 

.LZ.O) 

GO 

TO 2062 

IS 

(LSC,«) 

= 1 





DO 

2064 1 

^' = 1, 

V » • 
• « , 

IS 



Tt> 

(I SXIF ( 

LS,L 

SC) 

. LE. *') 

a C 

TO 2064 

I3S(L‘^,LS 

:,'i) 

= 1 





2 06U 

:iOSA (ISC, 1-1) =NCSA (L30,M-1)-1 
lYTA 3 (ISC ,M) =:' 

206 3 COVriN’IE 

w’ rrv C770) isc,’;,:<Ax? 

A9770 rop-'. .\"C All S/C USED - I.£C, N , 3* ',3*6) 

206 2 CONII'njF 

c 

C .8 15 

C 

DO 2 0*0 LSC=1,:i;:'L:5C 
AV?! (LSC) 

DO 2071 Xsl,M' 

IT (*S (LSC,.'I) . 2Q. C) 30 TO 207 1 
IiDYS (IVJ,MY) 

IF (-l-r.",^. ??3 (II-E-v (15,-^) ,y-:vEA?. (LSC, 'll + 1 ) ) 30 TO 207c 

C B'IS FAILS 

IS (7. 3C,.r) =0 
DO 2 77 2 LS = 1,r.A>:iS 
2072 ISS( LS,LTC,:1) =0 
GO 7 1 207 1 

207 8 CO.'JIINfJF 

C ? 16 

DO 2051 LS=1,XAXLS 

IF (isit::(ls,l«c) .LL.O) go to 2091 

I? (ISS (LS ,L5C,M) . NT. 1) GO TO 2082 
NDY5 (UUi-MY) 

IF (7 N. LF. ?r? (L£, N-IY’^ A3 (LSC, K) ♦!) ) 3CTC 2092 

N'OOP (L3,N) =rv*)np { LrS, N) +1 
GO TO 20-1 

2082 «i=‘AlL (LSC ,“) =Y?AIL (LSC,«) +1 
I7S( LS,Lr.C,K) =0 

206 1 CO NT IN Ur 

17 (:;OS5C(LSC) .GT.NrATl(LSC,F.) ) GO TO 209 9 
TS (L SC,F.) =0 
20S9 co‘'^iru? 

207 1 INUF. 

2070 continue 

C 

C CJ7.PUTE NOS\ (LSC,.'l) 


! 



n o o n n o o f* o no 


L^C’^.A/LHC 
L5=.'- AXLS 
MMAr=LS»::3*i 
N-ATa (LSC»LS) *M3 ♦! 
no 9 3S3 110=1, 10 
DO 3 352 125=1,25 

NZ »3 E»; C 1 : , T 2 5) = N-0.0 7N C 1 0 , 12 5) -NOO E ( I 1 0 , N ) 

I? (N EV'S’I'J (110,125) .L-.O) NEWSEH(I 10, I25)=0 
9352 CO'.riNUK 

DO 9 316 : 25= N 1,2 5 
J25*25f*n -125 

EK (II :,J25) = .‘:ZV5EV (I10,J25) -'J£WS5M(I 10,J25-1) 
I? (NEWSEN (riO,J2 5).LT.O) S EW 3E N (1 1 0, J 25) =0 
9316 co!:riyuE 
9 353 CO N’T IS UZ 

CALL sosAT (:-:A'’ri:'’,::rDE>:? ,i!JD£xc,.'i.'iAT, :;mat) 

I? (M A<?.!!£. 1) GC Tj 5901 
ypiIE (6, 332) S 
'W^IT E (6,3 2 9) 

WK IT E (6,3 3 3) (*I00P(:,N) ,1=1,10) 

DO 335 1=1,10 

335 wr.ITZ(6,331) ( .v £W ?E 1 (-' , ‘J) ,0 = 1,25) 

'rf^ITS (6,3 3 3) K3AT 
5901 DO 9 0T6 ISC= 1 , ;-'A ELEC 
SOSA (LSC, M) =KSAT (L3C) 

9006 CO’JIIirUF. 


13 COMTINfJ'; 

DO 5 3*^3 L5C=1 , r'.AXLSC 
PC83 y03( LSC,N) =SC53A (I,SC,M) 


SPACECTArT LOOP 
LSC= 1 

207 5 com*: UK 

I? (SOS (L?C,K) . TE .0) GO TO 20^0 
2075 LSC=L0C+1 

IF (L SC. E2 . KA. IL5C) GO "'0 2076 
10 TO 2079 


STA?. T 0? LAUNCH SEgJF NCL 

!• 20 
2080 CONTINUE 
3020 II=L :/C’»S (ISC, ’J) 

ILA( LSC,M) =I1.A (LSC,V) +1 
S’.!=? N’DYS (DUIHY) 

IF (7S. G£. EPS (II) ) 10 CO 2050 

C 300 STEF SUCCESS 

8S=S NDY5 (DMyr.V) 

I? {= S.GF. PDE3 (IT) ) GO TO 209 1 

C 9 0CSTKP FECOVESPD 

:p::( LSC,!I) = 13? (LSC,N') +1 
GO TO 2092 

209 1 IP (LSC,N) =1? (LEC,N) f1 
20?2 8t:=? MDY5 (DUrf-Y) 

IF (E N. PCS (II) ) GO TO 3023 
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, n AaCM' 

GO ro 300 1 

209 0 ?•:*?. HUY'S (nnM.-^Y) 

r? (? H. :J. i?Jr- (Tl) ) GO To 3000 
IR (L 9C,H) alR (L3C,r) ♦ 1 
GO TO 3''0 1 

IlOJS'^ER "FCO’/tFiD 

3000 ■"??( L3C,N'1 =I3K (ISC,N) ♦ 1 
300 1 3'-=9 f;DY5(D'»y.lY) 

TF ri :i. 3T. ?ChOF c -) ) ;r ro 3002 

IO«(LTC,F) =IOr (LSC:,.^) +1 
GO i’O 30 2 

3002 IPULFCm-.') =IP!.(L'^C,H) 0 
:o (L3C/:) =10 m 

I? {0 HBC.\P (LCT) .f W. 1) ;0 TO 3020 

:r.'i(Lsc,:i) =:?'i (l.sc,*:) ♦ i 
GO - 0 302 0 

302 3 COV"TN?tr 
. . . . . ^ 23 “ 

I? P ?3CA? (IFC) .F'". 1) 3C TO 3034 
?.!.’=5 MD Y5 (OUT 1Y) 

I? (•. N’. 3F. ?F3 (LSC.n iO TO 3030 

TAYLCAO OK 

’'.‘I=r NDYO (DOF 'Y) 

I: (? S'. ?CF (:•") ) GO TO 3031 

AiCC'/iiSFU 

io3(lsc,:t)’=io'’ (lgc^h) +i 
GO TO 303 2 

303 1 "0 (L3C,r») =10 (LSC,H) *1 
303 2 :PI(I.SC,”) =: *’L (L3C,N') +1 

'O TC 3039 

3030 SS's-’ M0Y5(D*.J''; -.V) 

IF (a 3. L". ?C? CD ) GO TC 303o 
(LSC.H) =10 (L5C,li) M 
i?L(L3c,;d =i?L(Iic,'I) -n 
10 302 '* 

chairs? -TCCVf.r.iD 

30 3 6 I0?(1TC,::) =I 0 F (lSc,N) ♦•I 
:pl= (lsc,w) =:dl? (Lsc,:n ♦ i 
GO 1 0 202 0 
3034 COST TV 'JF 


. . 2 4 

?‘J-=0 :.'DY5(D'f“.^Y) 

TF (^. H.GS. ??Fr. (II) ) GO TO 3040 

. . PM O'K 

?*;=? :.DY5 (MJ.\“Y) 

IF (9N.L". ??r? (IT) 1 TO 3046 

..?i Do;-:r 'jot diacf payload :i; 

RN='' NDy9 ( D'IKF Y) 

(F :j. LT. pjFx (II) ) 30 ro 3041 

OFbITF? TCI ?3CCV£r£D 

TO (LSC,M =:o (LFC ,N) ♦ 1 
GO ”0 3C43 


P FFOCVt’FKD 
T) =10? (I SC,! 
=i ?L(ISC,’ 
~z?^ (L5C,: 

20 


3041 


un 


b I 

ff 


IOr( LF 

c. 

304 3 

:^L( L.9 

^ t 


^ n V ( ^c: 

c, 


GO :o 

30 

3040 

con: in 

HE 


DCS 12 ID OtiUir 


1 - 


1 


1 
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C ?•* D'^'ES •n'T CHECK 01!" OK 

nils’? MI)Y5 (DJI'-IY) 

IF (3 !i. I.: . ?CH (TI) ) GO TO 30 u8 

c OFPiTc*. no: «£ccve:?ed 

10 (Lsc,;:) *10 (!.‘?c,N) ♦! 

GO ■'*0 30M 
C 

30<«3 COIlTI'inP 

C OarilFF PiCCVFRL’D 

IO?(LrC,!i) =:0F (LSC,N) 

IP. ''.3 (L£C, :J)=I?.'1F (L3C,:.)*1 

;o ro 1020 

C ?* PLACEli PAYLCAf) IN DFSIFF.D OP.niT 

C ■» 20 

c 

30‘*o CONTINUE 

r:i=? vdy*' (curiHY) 

IF r.r. LT. (L£C) ) 30 TC 3057 

C PAVLCAD DOE? MOT FU'1C"':0N 

?.’»=?. yDY5(DUiriY) 

K=0 



IF (’J 

N.OE. P? 

LEV (I 

1 1 

= 1 



F.'<=?. 

NDY5 (DUMMY) 





IF (? 

N. 3- . ?r 

.Y?P (I 

I) ) 1 

C T 

C 305 


PN=H 

*iDY5(DU 

!!MY) 





If (5 

N . 3 E . PC 

■» (ID 

) 30 

'”0 

3051 

c 

.OFUI 

TER .’’ECOVLrE 

n 




ICK( 

LSC,N) = 

TO; (1 

sc,.;) 

♦ 1 



IF (K 

. 3'-.'^) 

30 '’’0 

3052 



^ 

. . . 5 

L KA? R 

E.hCQU 

I-’EP 




I?LF 

(ISC,';) 

= I PI E 

(LSC , 

.' ■) 

1 


30TD 

3053 





3052 

:?L( 

LSC,;--) = 

I?L (I 

SC,'I) 



3053 

irr- 

(LSC,;:) 

=if;’F 

(LSC, 

;•) * 

1 


■30 r 

0 3020 





lOl** 

CON"" 

I *’ U E 





C 

PM J 

0 E S NOT 

RETU 

CN ”0 

O.R 

?IT£? 


?.M-7 NDYj (OrirMY) 

IF n N. IT. PC3 (II) ) GO TO 3055 
3051 I? (I SC, =10 (L.<C , M) + 1 
10 TO 304 3 

3055 If':{ L£C,’.) =10.-’ (L£C,") +1 
GO TO 3043 
C 

C PAYLOAD ?U Motion 3 OK 

C i 27 

3057 COVTi:i-JE 

FM*.’' NDY5 (DUMMY) 

IF (~ N.GE. PFMP? (II) ) GC TO 306C 

c ?M detjftss to OPDITEk 

nj;=S NDY5 (DU-^MY) 

IF (~ N.GE. ?rt (XX) ) -o TO 3062 

C opaiTE.l ®ECOVt?3D 

IO^(LSC,M)=IOB (l?C,:i) f 1 
IP.^A (ISC, N) = IPMF. (LSC, S) ♦ 1 
!?-.( L.<^C,N) =IPL (LDC,N) ♦ 1 
3035 CONTINUE 

NOS(LFC, .') =MOS (LSC,N) -1 
C 

C END OF SUCCESSFUL LAUNCH SEQUENCE 
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C* ..P'1 DDES Mn'i' CHECK Otl" OK 

rU=?, MI1Y5 (DliriY) 

I? (3N. I.T. PCH (ri) ) GO TO 30U8 
C. . . 4 . , . OHBTTU". !iO-r 1=£C0VERED 
10 (L SC ,V) =10 (LSC 4 N) ♦ 1 
GO -PQ 30 U 3 

3043 COilT-I'ItlE 

C,, ORBIIEE PECOVERL’D 

10? (LSC,li) =70E (LSC,N) +1 
iPER (Lsc,H)=i?:irML 3 c,:.) + i 
30 TO 3020 

4 4^ F:1 PLACES PAYLOAD IN DESIRED OHRIT 

4 4 . 4 4 ■ . 7! 2C 


30:^0 CON” IN UE 

RI1=?. NDY5 (DUrUlY) 

I? (R N. LT. ??S (LSC) ) GO L'C 3057 

PAYLOAD DOES NOT ?yNC'”IOM 

?.fI=E NDY5 (DUHMY) 

K=0 

IF N.GE. PPLRP (II) ) K = 1 
F.v=?.NDY5 (DUMMY) 

I? (? N. 3E. PPM?.? (II) ) 30 TO 3050 

: PN=H NUY5 (DUMMY) 

IF (E N. 3E. PC? (II) ) 30 TO 3051 
. 4 4 4 4 OR UTTER KECOVLrED 

I0R(LSC,N) =10:. (LSC,;I) +1 
I?(K.3''.0) GO T’O 30 52 
....... 'L KA? REACQUI'=ED 

IPLR (ISC, ’I) =TPLP. (L3C , N) + 1 
GOTO ■ 3053 

3052 I?L( LSC,N) =I?L (LSC,N) +1 

3053 IPM? (L5C,N) =.tF:!F (L5iC, N) +1 
GO TO 3020 

305? CO’” IN Ur: 

DOES NOT RST'dlN "’0 ORBITER 
?.l:-=R NDY5 (DUMMY) 

IF (RN. IT. PGR (II) ) GO TO 3 055 
30 51 IT (L SC/N) -IO;(LSC ,N) +1 
30 TO 304 3 

305 5 IC:(LSC,!;)=r02(LSG,N) +1 
GO TO 304 3 


.PAYLOAD rUNCTIONG OK 

444444 . 02 / 

3057 CONTINUE 

?N=? NDY5 (DUMMY) 

IF (?.N.GE.??f1R?(II)) GO TO 306C 
...4, ..?M returns to OPDITER 

r;II=R NDY5 (DUMMY) 

IF (? N. GE. PCI (II) ) GO TO 3062 

ORB ITER RECOVERED 

IO':(LSC,M)=IOP (LSC,:i) <-1 
IPMR (LSC, N) = I?MF. (LSC , N) + 1 
irL(LSC,N) =IPL (L3C,N) +1 
3035 COM INUE - 

NOS(LSC/:;) -NOS (LSC,M) -1 

END 0? snCCRSSRUL LAUNCH SEQUENCE 
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) 






I 




It 




} 


C 

GO ro 207 9 
3 OC 1 O CONTIfiUE 

c Pn DOES '50? ESTUEN* TO 0 P.D I TEH 

HN=R t!DY5 

IF {? »r. LT. POP (II) ) GO TO 3063 

3062 10 (LSC/?,I) =10 (LSC,-N) +1 
■;o ro 306 4 

3063 ICP( LSC,in =lOr. (LSC,N) +1 

3064 T?f1 ( L3C,4) =!?•-! (I.GC,T) + 1 
i?L(isc,':)=T?L(LSc,:r ) +1 
GO TO 30 3 5 

C 

C 

C RMD 0? I.UIKCH .SEQUBTJCH 

f 

2076 C0 2;7INUF 
C 

c • » » » t ^ 2 1 , 

c 

FAC= FLO.s.T 1) **LCHR 

DO 3 070 L3C=1 A XI.SC 

C(1, L3C,N ) =XLA(LSC,N) ♦LOST (LSC, N) 

c (2,, tsc,‘: ) =irL (LSc,N) ^secsT (lsc,:i) +iplh (lsc,io 

2 sees? (DSC, :r) *?AC 

C(3,LSC,S' )=C(1,LSC,N )4-C{2,L3C,N ) 

DO 3071 K=1,3 

17(2 (K,LSC,N ) .Gr.AM>:(K^LSC) ) A«:<(K,LSC)=C(K,LSC,I'I ) 

I" c (K ,LSC,K ) .L7. A^::: (K,LSC) ) A/N (K,L3C) =C (K, L3C, ) 

S'JF.C (K,L3C,F) -SUMC(K,L?C,F)+C (K/L.SC, M ) 

SU'iSC {K,LSC,N)=S'Ji'!.:.C(K,L7C,M) +C(.K,ISC,N' ) **2 
307 1 CONr-INGS : 

3070 COB'^I'.’fjZ 
r 

DO 3 090 K = 1, 3 
A=0. 

DO 3 09 1 L.3C= 1 , :1A /[LSC 
3091 A=A+ C(K,L3C>>r) 

SnSOT (K,N) =«U.‘!CT (K,?:) +A 
3090 SUMSC^(K,M)=SUi-;sCT(K,j-:)+A**2 
C 

DO 3 07 2 K = 1, 3 
GT (K , V )=0, 

DO 3 073 L3C=1 ,'-!AILSC 

CT(K,N )=CT(K,N )+C(K,LSC,N ) 

307 3; CO'ITINUS 

) .GT.A-IXT (K) ) a:iv^(K) =cT(K,^ ) 

IF(GT(K/;r ) .IT . A3T (K) ) A'IT (K) =CT (K, H ) 

307 2 co;;tt?^uh 
c 

c' v;'.;,:;:-'- ■ ' ■ • ' . ■ : , ■ 

C» . .V. . . 32 3 

3038 CON’TINrjE 

DO 4 000 J=1 , 9 

ir =0 

DO 4 001 L3C=1,f1AXLSC 
4001 H=I I + laOST (LSC, M,vT) 

(J^N) =£U:1J {J ,N) +II 
S'JflJ J (J,!,') ^Srj.’1SJ (J, M) +11*^2 
IF (II.GT. 'lAFI) II = KAXI 


Ir* (I I. 3T. J\Z) -1*1 : 

aooo c^MririiiK 

DO 4 009 LSC*1, -ik /L5C 
II*: LA (LSC,K) 

SMU A (LSC,!?) »5-Jf»LA (LSC,'I) til 
SU.’ISLA (LSC,f») =3U:iSlA (ISC, N) ♦11**2 
I7(II.3T,.1AXI) II*MAX: 

I j =: I ♦ 1 

?LA( IJ,LSC,*:) »PL A CJ,L?C, N) ♦ 1 
IF (II.G-^. MA (LSC,N) ) * LA {LSC , :i) =I I 

II*: PL (L3C,::) 

suMPL (Lrc,K) (L3C,*;) ♦:! 

SU«3 PL (L3C,N) *£U-S?L (LSC,r) ♦II**2 

If (II. IT. .*Axi) :r=MA"i 

ij*:i*i 

PPL(T.],L3':,H) *??L CJ^LSC, J) + 1 
IF (ii.GT. -rL(L£!c,ir) ) -rL(LSc,::) *n 
4 009 CCSriVUE 

Aljpy »ANrv +CT (i,S) *P7? (!J) 

9002 COFTIN'UF 

c 

C 2*!D C? riKF LOO? 

w 

SJ.*? V = SUM?V + Ai:PY 
t? (A ;;9V . GT. PVX) PV»A»I?V 

ir (ANPV .LT.PVV) P'/r.-^AIirV 

O’JXS ?V*S’J.*'.SPV + ANPV**2 

C*****0'J7?rjr A7.RAYS TC DISK AXD ClilAR SI0SAG2 

w 

W2I? E (2) C, Ai:pv 
9001 co::“’ 

rewind 2 

c . 

C END OF '10*:TE CARLO 

c 

c 

c 

00 4 010 K = 1 , 1 

cs'^i r (K) = (AF.::-' (k) (k) )/»:oi:.' r 

DO 4 010 I.3C=1 ,.'iA:{L3C 

cst: (k,l5 0* (A^;;^ (k,lsc) -axn(k, lsc) )/scint 

4010 co:rr\r»i 

Fv:= (Pvx-?vj;)/NOi:;r 

n 

C LOOP BACK OVSP DATA STOF.ED ON D"^SK 

C 

DO 5 000 KP = 1,*!Ay? 

. RRAD (2) C, CT, A Vpv 
c 
c 

DO 500 1 K = 1,3 
C 

IF (CSTT"^ (K) . r,r. . 001) HO TO 4032 
00 4 03 3 !I = 1 ,i'!A::K 
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U033 ?? (K , 1 ,*•) 

;0 T'' U '>3 7 
4032 DO 4 035 1, 

L*(CT (K,N- ) - A'*T (K) )/C3?I? (K) ♦. 3?'»9 

Ifd.^O.O) L»1 
Ir* (I .GT.*. ?T!T) I»MO:!IT 
PV (K,I ,>I) 

4035 CO'J'I^UE 
4034 ccN-i:;u2 
4037 CDI.TITJE 
C 

O') 5 '>0 2 lSC»1,.M/.'-'13C 
CSTI K«C.'iri (K ,l«C) 

A'ir::»Ay;i(K,LsJ) 

ir (CSTTK. 31. .ooc 1) 30 TO 500 J 
Dv*' 4 02 3 

4023 ? (K, irc, 1 ^^y' 

GO :0 50" 2 
500 3 DO 

L = (C (;<,L?C, -J ) -A'lW (K, LSC) )/C3T: (K, L3C) ♦.9)39 
I?(L.S;>.0) L=1 
!*■ (L . L»Kc:r'r 

P(K, L5C,L,r)^.MK,15C,L,:J) ♦! 

500 4 :o*r:r:rE 
500 2 coin ■105 
5C01 co‘j::*>rjii 

I’ (?v:.l: . . 001) 30 to uo4 9 

L-»(A*.TV -?VK)/?V:^.9?'5 

r?(L.?o.:) L=i 

ir (•- .GT.-:oi:iT) 

??V(T.) (1) ♦! 

10 TO 4''42 
4049 c"'*.: I'l'jE 
?^7( 1) 

4(^42 c'.»rr^::uz 
503 0 


33 


AMA*: A". ?. 

T? (DL7C.N7. 1) GO TO 4 0 50 
D<'' 4 051 K = 1 , 3 
no 4 051 LGC= 1 , u;;lsc 
DO 4 051 = 

r.Tsr (K,Lsc,K) =ii;.'-c(:<,iST,:!)/AMAy®. 

«Tik:s^ (K, loc, ':) = d:-:v (gi;:-:sc ( ti , l5C,'I) ,.3CST (k,lsc, n) ) 
4051 crtN-^ :'’ur 
4050 CCIiriEOT 

DO 4 060 3 = 

DO 4052 J= I , 1 

y.E (0 , 'I ) =7 U y J ( J , N ) /A y A X ? 

5VDM (J,V) ttl«FV (3UKSJ (J , N) 

405 2 CO N'T iron 

DO 4 05 3 K = 1, 3 

« ,N) =?'!VCT (r, N’l/ATAXH 
3TDT C<,N') =DZ'/ (EU''SCr (K,E) , hT (K,N) 1 
4 05 3 TO'.T T’”jr. 

DO 4 054 L.5C= 1 ,::A 'LEC 
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» 




I 


I 


«LA(L£C,:J) »S1I»»LA (LSC, M)/A?»AX5 

srn A (LfC,N) *J?V (P'J'lJIAtL^C,*!) ,rLA (LSC,H) ) 

nrL(LSC,N) «SUSFL (L:JC, S)/Ar.AX? 

S^'D? L (l?C , ll) * JKV (S't/.S?L (L5C, '!) ^'IPL (L3C, S) ) 
4054 CO-'"I'»U£ 

4050 CONfl’IU^ 

p?v»s!i’*rv/AXA>:r. 

STD?V*Dt:V (£Ur.S?V ,.' 1 ?V) 


*♦****♦*• ?fc?0?TS 


CALL Oy^P'jT 

src? 


K>v'T3 


S'lB/. OU*:!”." OUTPUT 


3 TAL KPV 

C 0 S'* ON/') II T 3/ K : , ? p , A 

1 riTDCrT,Pl\,rA.<K,/. 

2 ,PT,PPL, ^.?L,STD?L, 


vr,r?,*:TDE,AhT,:i 0:?:T,C£T:7, r,«T,STDT,r.C3T 
1 A , « 1 0 L A , ? V M , P V I , ? ? V , K ? 7 , S T D ? V , 1 1 L , D f . T C 

cst: , rAXLSc, of 


0 


r*«TiSiJ?F 4 ->£TC 


DiK£Ms:o:; 3,^) , lai>:l (3,8) 

H£AL AIN (3,5) ,CS7I (3,5) 

H2AL r-i- (2 5, ^,25) ,.'*£(?, 25) ,srOE (9,2 5) ,AXT (3) ,CS TIT (3) ,? (3,5,4 0,25) 
S3AI XT(3,25) ,£T0"(3,25) , «C?T ( 3, 5 , 25) ,37DC£T (3 , 5, 25) ,ELA (25,5,25) 
PSAL ?:?L(5,25) ,S IDFL (5,25) ,??V (40) 

?*:AI, r.LA (5,25) ,?TDLA (5,25) ,r: (3,40,25) , PPL (25, 5,25) 

AXF.*r.Ay?. 

.'ll + i 

DO 511 LSC»1 , XA/LS: 

0) 5 11 Na1,.v.\XN' 

DO 5 11 r* 1 ,'*! 

611 PPLf.r ,L3C , .V) »?Pl ( %LSC ,K) / Avp 
511 PLA('!,I5C,;i) *2IA (M,ISC,N)/ AXr. 

Dr 1000 '<=1,3 

?2AD (5, 1001) (1A2 2L(K,I) ,2=1, a) 

52A0 (5,1001) (LAPrL(K, I) ,2 = 1,8) 

10'>i AT (9A4) 


DO 4 44 :i=i,.xa:'.n 

DO 4 uu MO = 1 , NOirT 
DO 5 02 L‘^c=1 , XAVLSC 
502 P(X, LTC,:;0,T)=?CK,LSC,NT,':)/ AX?. 
444 (K ,*:C,'.) =r'T (K, <0,N') / AZ? 

1000 CON"IM'!»: 

DO 2 I’1 , XII 
DD 2 J=1,9 
DC 2 ‘;=:tl,.xa>:m 
2 Pi c ,J,<) =?t C,J,M)/ .UP 


c 

DO 1 S = I7L,.XAx:J 
WHITE (6,999) 

999 F9F.1 a:* (i:-: 1) 

W^ir E (6,70C) 

700 FORM A” (29.", ' E5TAEL2S!!.XFN7 AKD .XA2fi:SSA:JCE ' ) 

M-^irr ( 6,75 1 ) 

7C1 POriAT(15M,'OF SYSTEK C" EA5TH CKDITIU3 3PAC2.CP A FT ' ) 
U^’ITF ( 5, 1 002 ) 

100 2 F 3 K:iAT (/1 

wr.n E (6,7 02) '• 

WHITE (6, 1002) 

702 ror-AT ( 6 oy,» yean : • ,ix, i: 2 ) 


I 


o n 
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WRITS (6, ■^03) 

703 AT (20!t, ' PROP ABILITY OF INDICATED O’JA^TITY') 

wp.r E(#>,T19) 

DO 3 ii»i,nii 
I*ST -TI^ 1 
J*I» 1 

3 Wti:E(6,7 0U)I, (t F(J,K,i:) ,K«1,9) 

70 4 ?OR^ AT (lA,i:3,9y.,9F11 . 3) 

WRr »: (6,705) 

705 pory AT (17>: , ' LAUr JIP ,05y, » ADDIT, ' , 05X , • ADDII. • , 05M, • ADDIT . ' , 05)C, 

1 • SPACE IE A FI' , Ol • BOOST LEi' ,03/.,' 0P3 IISES' ,0 3.{ , ' ?X3» , OBX , ' SFACt 

2CEA.-T') 

W?ITS(6,706) 

706 FOR'1 AT (1'.'*’', ' ATTFPPT3' ,03!’, 'BOCSTTiiS' , 03X, 'ORBIIiPS' ,03>:, ' PXS ' ,03:<, 
r, '?. saniF LD' ,03/, ' EFFOaiJ. ' ,04X ,» R2.»U3D. ' , 04X , ' F RFU 5 . ' ,C4X, • F*FUr& 
4 . • ) 

WRITE (6,7 07) 

707 F0.3.'' AT (3'':<, ' r.Egy IP.2D' , 03/ , ' E2QUIE ID' , 03X , ' IvSQniF.ED ' ) 

W?TI K (6,70«J) 

70S ?0:;iAT (OV',' EXPFCTED' ) 

WRIT F, (6,709) (1S(T,:0 ,I«1,'') 

709 fOP.RAT (Or/., ' ,5’',9?1 1. 2) 

1 WFi:z(6,710) (5TDE(I,N) ,I'1,9) 

710 TOP/ at (0i:-:,' 3T0. DFV. • ,02y,9?1 1.2) 


DO 4 K*1 , 3 
(K) 

A17K=A“.T (K) 
wri.? E (6,?99) 

WRT^K (6, "CO) 
wr-.rrn (6,701) 

W?IT 2(6, 1 002) 

WRITE (6,7 1 1) (LA3.-:L (K,T) ,I«1,?) 

WRIT E (6, 1 002) 

711 rC?y A" (2C •:, ' 7ROEABILIIY OF INDICATED ',9A4) 

W.-':'K (6,7 12) 

712 ■’0?>: A? (051, ' COST '/05>: , 'F.A'lGE'/,04y. ,' (MILlIOKS'/,OaX, »C? DOLLARS)') 
DO 4 4 r:c=i,NOiVT 

i=j- 1 

A=Ar " K*I*CK 

D* a;i tk+j*ck 

4 4 WRIT K (6,7 13) A, B, (P^(K,J,H) ,:i=1,:iAXN) 

WFi: Z (6, 1002) 

WRI'E (6,921) ri,N = 1,.nA':N) 

921 foe:'. A" (02:v,' YEA? • ,13X, 1019) 

W? II £ (6, 1002) 

71 3 70?.:: a: (l::, 1F7.2, , 1?7.2,6X, nFS. 3) 

WSrF(6,71U) (FT (K,N) TAX*:) 

714 POP-! AT (05 •*, ' F/ ?tCT'-:0 CC.STC ' , 0 IX, 10? 9,1) 

WRIT £ (6,7 15) (STD- (F.,*!) ,N=1 ,.-AXN) 

715 70PF, AT (05::,' STD. DC V . ' , 9X , 10 78 . 1 ) 

4 coi;?i::uE 

c 

c 
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I 


I 


m 


-0. 


I 

I 


f -4 ) 

I 


IP (D ETC. L‘Q.0) GC TO 5 
DO 500 K»1,3 
DO 5 00 L?C»1 ^.lAyiSC 
CK«TS'^r (K , L«5 :) 

AK*A (K, LSC) 

W?TT Z (6,999) 

Wf'^T’; (6,7 00) 

W?2T S (r.,‘»0 1) 
wr:t?((,, 1 002) 

wo:ri(6,oii) aAb-L(K,i) ,:=i,‘^) ,lsc 
911 FC?.'*. AT (20/ ,' .’.-OE ability OF TJ.’DICATED • , 3 A4 , 112 , • ) ' ) 

E (6,1002) 

W?i:E(6,712) 

DU 5''! :.’0* 1 , MOT”"' 

j*:.'0iNT-rro-n 

:=j- 1 

A-= A< ♦! *CK 

B=AK ♦J^CK 

50 1 W?i:2 (6,713) A,B, (P( K , LSC , J , :i) , :i= 1 ,M Ax:;) 

W?.:T E (6, 1 002) 
wf::e(6,921) (:j,:;*i,.-iaxn) 

(6, 10 02) 

WEI- H (6, 7 14) (''C57 (K,I SC,.'i) ,N*=1 ,.iAIK) 

V?ITr.(6,7l5) (ST0CF:(K,13C,J.) ,:I=1, ^AMN) 

50C coir?3'nr 
5 corriNyE 


DO S 10 LS C-»1 , MAY LSC 
W?!’' E (6,999) 

W?I^ L (6,7 0 0) 

WPr E (6,7 01) 

-FITE (6, 1 002) 

'-r.ITE (6,7 16) LSC 
V3i: r (6, 1 002) 

716 FO-'! (40 >;, • .S?ACL'C?A?T TYPE :',1I3) 

WFITF. (6,7 16) 

716 FO'"!AT (10T, ' PHCBAniLITT 0: IIILICAriD LAUNCH AIIL/iPTS') 
W^ir F (6, 1 002) 

( 6 , 7 10 ) ' 

7 19 FO’! A" (CO:*, ' QUANTITY*' ) 

DO 5 12 Y.= 1,XI1 
I=Nr-N+1 
J=I+ 1 

WrlTE (6,720) I, (PLA (J,1SC,:.') ,N=1,«AXN) 

720 FOPXAT (03X, 1I2,10X, 10F6. 3) 

512 coNTiNrir: 

WEI?E(6, 1002) 

WBT'E (6,7 21) (N,K = 1,;iAXN) 

72 1 FOF.yAT (02 X, ' YEAR • ,06X, 1016) 

W?i: E (6, 1 002) 

WFirr (6,7 2 2) (UA (LSC, N) , E - 1, N A XN) 

722 ’C?YA' (31, 'EXI-£C^EO'/3X, • NUXflEK' , 5;<, 10F6.2 ) 

VSITk (6,723) (STDLA (L3C,M) ,N=l,f1AMN) 

723 ?0?N A"" (IX, 'STD. DEV . * , 2X , 1 0?6 . 2) 

510 con:in!ie 

c 

r 


reproducibility of thb 
ORIGINAL PAGE 18 POOB 


f") n n r> 


DD 6 10 LSC»1 ,.1AXLSC 
• Fir E (C,99?) 

WPir 2 (6, ■’ 00 ) 

WFIIE (6,7 01) 

WFirS(6, 100?) 

WRITE (6,7 16) ISC 
WRITE (6, 1 0C2) 

WRITE (6,7?5) 

725 FOr."^ AT (10"' ,’ ??.:2 ABILITY 0? INDICATED SPACECRAFT REQUIRED*) 

W?:r 2 (6, 1 00?) 

WRITE (6, 7 19) 

DO 6 12 r- 1 ,“'1 
I = .'?I-!i + 1 
1 

V?1TF.(6,*’20)I, (Eri(J,LSC,::) ,N=1,rtAX!l) 

612 COMTIN'UE 

WRIT E (6, 1'702) 

WRirt (6,7 21) (V,N = 1, ’IAIN) 

WRITE (6, 10C2) 

WFI7 E (6,722) ('’.PI (ISC, S) ,N = 1, MAXN) 

WRITE (6,723) (FTDPL (LSC,N) ,N=1, KAXN) 

610 CONTINUE 


DO 3 O'' N0 = 1,NCT»’T 
500 ?'>V( NO) =PPV (NO) / A 'R 
W»ir E (6,0 99) 
wriT-E (6,7 00) 

WRITE (6,7? 1) 

••--<:te (6, n? 2 ) 

w ‘7 1’ !•: (6,726) 

726 at (10'*, ' ’’REE 'ENT V A lUE ’, 1 3F 30b A BI L ITY OF INDICATED') 

Wr r E (6,7 2"') 

027 FCR/. at (10X, * (''II LION? Or DOHA PS) ' ,0 5X, ' PP.ES’EN T VALUE 0? COST') 
WRITE (6,72?) 07 

729 F071 at (36/, • (DISCOUNT ?.ATZ=' , 1 ’’u . 1 , ' ) ' ) 

WPi: ^ (6, 1 002) 

DO ■’2 9 N'' = 1 , 

J-=P0IN'7-V0 
I = J+ 1 

A=PVN + J<’?71 


n=?v N*i*?v: 

WRITE ( 6 , 709 ) A,!3, PPV (I) 
7 2 9 CO N7 I N U E 


799 

790 

791 


wot; E ( 6, 1 002) 

?OHXAT(10/,1F'i.2,'--' , IFS.2,16*', 1 F5 . 3) 
WFi: E (6,790) N?V 

FO?.X at (10**:, ' E/.PECTSD ?V',19X, 1F3.2) 

WRIT n (6, 79 1) ::ropv 

=‘OR-l AT (10.:, • .RID. DEV. ' ,?2X, 1F8. 2) 


REIN PN 
END 

FUUaCUTINE A LIST 1 (SUN PL, SU KS PL , P , PPL , PT , PF V, AN P V, S U.UP V , S U.NS ? V , 

1 1A-; LS,'^:\;*’N, .''.a:CKB,.NAXLSC ,?FS ,FF3,SUI1LA,SUFiSLA, MLA, 

2 PL , .9 UMC , SM N.NC , PLA , S UX J , S UO S J , P E) 


FF.AL 


SU.NPL (5,25) , GUNS PL (5,?:.) , P (3 , 5 , 60 , 25) , PPL (25, 5, 25 


I 


> 

i 




f 


r 

i 


i9 ■ 


I? 


I 

i 

] 


i 


1) , n (3,40,25) ,?PV(40) 

RFAL ('3,2 5) , SU;^SJ (9,25) , PR (25,9, 2 5) 

REAL r>FS (10,25), ?73 (5 , 25) , 5UIUA ( 5, 2 5) 

PEAL Sl’r.i^LA (5,25) ,5U^C (3,5,25) , Sn.ISC ( 3, 5 , 25) ,PL A (25, 5, 25) ,^.LA (5, 2 5 
3 ),rPI(5,25) 

DO S''0 LO*1, 10 
DO fl O'^ !:=» 1 ,25 
300 ?75(L.P,‘:) =0, 0 

DO 301 \’5=1,5 
DO 30 1 N» 1,2 5 
801 ?F3(S3,:n=0. 

DO 3 02 M»1 ,25 
00 3 0.3 l;*c=i, 5 
S»If<? L (LSC , m) = 0 
PL (La - ,N) =0 
SO ''LA (L=IC,'’) =0 
S'J''.o LA (LSC ,S) *0 
MLA( L*^C,’i) =0 
f1PL(LPC,*i) =0 

DO 304 K= 1 , 3 
DO 2013 L=1,40 
2015 ??V(L)=0 
2015 P' (•■< ,L,‘l) =0 

2013 P (K, L?C,L , V) =0 

so::: (k,lsc,n) =('■. 

804 (K,LSC, :>) =C 

DO 3 05 TO*!,. 3 5 

2014 p?l(i:. ,L.JC,N’) =0 
305 ?I \(IL,LSC,:.) =0 

90 3 COX'" IN Mi 

ro 3 05 j= 1 , ? 

5:u.'!j (J,‘.’) =0. 

J (.J,V) 

DC .3 06 11 = 1,20 
90 6 ?E (I :,.j,::) =0 
80 2 O'! MIN OF 

.a;i ?v =0. 

SU.''P v = c. 

Sll.''.5 PV = 0. 

SET” PM 
F.MD 

SUUPOUTINE ALIOT2(COLT, 

1 SAM LSC, .S y.IS , KA ML , IM , SCC?T, LC 52 , I LA , IB , lU IPL ,10 , lOP, , I ?M , I PLR , 

2 IP-S, SrCP, NFAIL,IS,I3S) 

PE.AL IC5T (5, 25) ,CO'’T (5,25) 

TNTEO-r IX (5) ,ILA(S,25) ,10 (5,25) , I 3S (5,25) , I?L (5,25) ,10 (5,25) , 
nor, (5, 25) ,iry (5, 25) , IPLF (5,2 5) , r?;iP(5,25) , NOOP (10,25) ,.;faII (5,20) , 
2 ir. (5,20) ,I3S (10,5,20) 

PEAL 3CC0T(5,25) 

DO 4 20 LFC=1 ,5 
420 I>:(LSC) = 1 


I 

I 

1 

s 

j 

I 

I 

m 

I 

X 

1 

1 

1 

I 

"i 

1 

I 
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4 




r 


» 


5 




c 

c 


2008 


DO 1 OS9 ?: = 1,25 

DO aooe LFc=i,5 
CC2T (ISC, i:) = 1. tic 
sees? C3C,N) »1.E10 
LOST (LSC,*:) =1,210 
IIA(!,.?C,V)=0 
IP (lSC,':) *0 
:hf( L sc,") ^0 
IPI(13C,.\) =0 
10 r.sc,:;) =o 
IOF.(L?C,:!) =0 
i??:(L5r,‘ ) =r 
ipi.^ (I sc, *0 
iPf'?. (L3C, •:) = : 

CO NT Ivor 


2012 

109® 


DO 2 012 L3-=1 , 10 
NOO? (1S,>:) =0 


CO N? I’.-’TF 
DO 2010 M=1,20 
DO 2010 : FC= 1 , 5 
Nf AT L (L£C , 7.) =C 
I? (L SC,'1) =0 
DO 2011 1S=1 , 10 
2011 :S£( LS.LSC,'*.) =0 

2010 comvu" 

r-.OTL' PI* 

END 

SnaTOOTINE D?L"'T''f 
1*’OS2N, TSMIX, 13:^1 ■,A”T5?, 


calc,;. 


r 


AN I) L , A .n 2 AM L , A T 8 S , 3 7 A ND £ , A M E Al.’ S , A A a C , 

?. r» _ ^ n A T 55 A 'iJ ? «.• M A k. r- »• u ♦ «« • 


, Ai'ilMu , I f?3 , CALS, A E.Ayi,A^; i::7 , IPC.X, CAL3 , 


701 
700 
70 2 


2 

3 A"! A .. IC, n VIC, T 7~1C , CALLc, L NCiS , JR BCA? » ?PS , p 8 S , r hH£ , PBRF, PCS , 

3 ?7 FOE , p .■> r r r , PC p ! c ,??.•:??, ? ?lf r , i a v* , n , dp. , d etc ) 

IV-CG?^ ::0Fr::(ic,?5) ,:s.Y-v{io,S) , ib;:: (5 ) , :pt ( 1 0) ,:?C3(5), 

1 i^':LC(r>, 1C) ,li:ch5 (5,25) (B) 

?EAI CALC (1C) 

PC’iL ??;-?(10) 

KE AL (10 , «;) , TTANDL ( 10 ,5) , A'lTANL ( 10 , 5) , AIITBS ( 10) , STANDS (10) 

’^LAL A''C.\::S (10) , AVAXC (10) , A.'iINC ( 10) , A PAX 8 (5) ,AXTM8 (5) ,CALD (5) 
vFAL ;XA::LC(5,1C),A.NIVLC(5,1C) ,CALLC{5,1C) ! P?S (5) , PDS ( 1 ( 

1 ) »P D.S ( 1 n) ,. DPr ( 10) , ?03 (1C) / r08O?(10),PORy?(10),?OP. (10) ,PF«CMC) 
INTEGCr; CF=C.\?(5) ' /» 

REAL PP!<5 ? (10) ,rri.FS ( 1C) 

DO 100 LS= 1,10 
DO ■’Cl N = 1 , 2 S 
moss:.* (l.3, :i) =0 
DO 7 00 LSC= 1 , 5 
Ts.vi (Ls, L5C) =0 
DO 7 02 LSC = 1 ,5 
(L“C) ~1 

DO 701 18=1,10 
DO ■» 0.1 N?.= 1 , s 
AXTOF n 3, ML) =5. 

S^AN DL (L3, NO) =.5 
AXLA NL (I.T), »:in =5. C 
DO 7 on ir= 1 , n 
AF"3S (I.S) =5.0 


7G3 


STA:i (LS) *^>.0 
A.'IFASS (LS) '5.0 
A!«.A:'C (LS) =0 

Artl'iC (LS) «') 

IPT( LS) '1 6 
c»i': (LS) =?o. 

DO 1 05 <0=1,5 
h»\y. ti (!t3) *0. 

AKTS 3 (!:3) =0. 

IPT3 (< B) = 16 
CAL3 (?:l) =90. 

CONI I” OS 
AS xy K'O. 

ASIN .^'0. 

IPTS = 16 
CAL? ='^0. 

DC "»07 LSC = 1 ,5 

DO - 0*’ 12=1, 10 

A<A/LC (1SC,TS)=C. 

AS.IV 1C (LHC ,I i) =0 . 

TPCLC (ISC,IF) =16 
CALLC (LSC ,’F.) =00 . 

7 CCS? I '’0 2 

DC 70? LSC = 1 ,5 
DC " 0? S= 1 ,25 
9 L”CHS (L32,S) =1 
DO 7 10 L3C = 1 , 5 
0?3C P.o (LSC) =0 
0 PPS (LSC) = 1 .00 
DO “'ll IF= 1 ,10 
?3r(I?)=l . 
pBPs'rr) = 1 . 

?n?? (iL) = i. 

POf (IS)=1 • 

PO?CF (It) =1. 

?0?3? (IE) =1. 

PO?(I.L)=l • 

PP10 ("E) = 1 . 

Pf MP (IE) = 1 . 

?py.5 P (IS) = !• 

)1 ??I??(IE)=1. 

DO 7 12 LEC'l ,5 
12 MAM'*. (L SC) =10 
KM =5 
D’r* = l 0 . 

OETC=0 

?rTUP.!’ 

SND 

^ q f5 ^5) .\rtI3F(10,5) ,STA:JDL(10,5) ,A.1EA1:L(10,5) 

SFAL F^?(10,25) , PELS (10,25) , AS rSS (10) ^SIAM E5 (1 0) , AKEAKS ( 1C) 

Rvl"'= SC»I{^*0) 

DO icoi n=i,.*ia>;k 
F=N 

DO 1002 M3=1,SAXMD 
A= 1. 


n o o — non 


* 
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DO 1 003 LD=1 ,MAXLD 

R (?,AXTP? (LE,N3) , A.>1SANL(Lo ,N9) ,SrA:iDL (I3,WD) ) 

1C03 A*A*RFL 

pPl.U (NB,K) =A 

??D( N0,N) =1. -F.ELD (NB, M) 

If (i: .ro.1 ) 00 TO 1002 

=1. 

IF p.ELB(VB,M-1) . JT. l.E-50) FFti (MD,I.) = 1 . - PELB (N D , N) /B tL3 ( SB, :i - 1 ) 
100 2 C0STISU2 

DO 1011 LS = 1,:iAXL5 

RELS (LS,::) ="» (f , AMT33 (I.E) , ASEAWS (L 3) , STANDS (LS) ) 

P?r‘(LS,‘*’) = l.-rELS(LS, .1) 

IF . ZO. 1 ) JO "0 10 1 1 
PfS(LS,*:) =1. 

'I (?E^ S (1.3,11-1) .G'. l.E-50) ?FS(LS,H) *1 . -5KLS (L 3, N) /P.ELS (LS,K-1) 

ion co*’:i?:uf 
1001 coririKUS 
P.ETOfv 

£»: J3 

sun O'lTIS Z N03 AT (KATE IX , I'.iDErCF. , IN’ D£XC , K KAT , SEA T ) 


C 

c 

c 

c 

c 


c 

c 

r* 

c 

c 


i*:t3GL~ sat,df.m 

PEAL *5 Ki T-.I>; (.'■'KAT, SKAT) 

i»;TiGEF uebog 

i:!T£G’fR + 2 isnsy.c (Ki-iAT) , INDEME ( KMAT) 

:0 "" 0 N/ IT n A 5/CC S I ( 5 , 2 5) , S S S , il S A T , N Y i A R , S T I X E , 

S AT (10, 5) ,DEK (10,25) 

?OK'10>:/IA'Jf-0R/’*tIi^SAr (5) ,TTCST 
C0.X1 OK’/®? I/DEBUG 

CO‘‘F'JTES ?I;e KOKBLE CJ satellites 0? TYPE J TO LA'JNTH 
IN THE GIVEN INPUT YEAP 

INP'’T 

COST:: TOTAL COST 0? S/C OF TYPE I IN YLhP J 

nsfn:: NU*’n~? of sen’s^'.® types 

NSAT:: NUKIiEF CF S/C TYPES 

MYEA?;: NU'IBFR OF YIAES TO LOOX AHEAD 
‘rTIK!’:: INDEX FOE THE CURRENT YEAR 
SA";: •.'UKB-? OF SENSOF 0? rYPc T ON S/C CF lYPE J 
DE‘':: DEKAND FOR SENSC?. OF TYPE I IN YEAR J 
OUT? U*” 

NUMSAT:: N’UKBEr. OF 5/C OF TYPE I TO LAUNCH FOE THE CUrrENT YEAR 
TCOST::CCST CF LAntJCHIFG THE TOTAL SUKPEi. Cr 3/C 

MVAR=V'IAT-1 
NTI'I ER=N?irF-1 

ZAP .KATSIX 


D01IZAP=1 ,NHA'’' 

DO 10 7.A?=1 , NNAT 

riATR IX (IZAP, JZAP) =0. ODO 

INITIALIZE INDEX FOR POW AND COLUflN VARIABLES 


D0 2JC0L»1 ,NKA? 
INDEXC (OCOL) =JCOL 


151 


* 


» 


I 


I 


I 


I 


» 


I 


DO 31 R0W=1 , 

MVAr:=>’VA? ♦ 1 

IN OE X" (I?OW) =N’VA 3 

GENERATE SI:!PLEV HAT3IZ 

DO 10 C”YRAF.*1 , VYEAR 
HT I.i EE*”"’*! M2 ? ♦ 1 

lYMa : yea?- 1 
JCCL 1 = IYM<=NSA7 
NCOL=JCOL 1 

JCOL 2*IY '’^NSEN + NSAT”'!: YEAH 
JCCJL 3=JCOL2-M:l£N 
I\CV 1 = IYM*N'S S” 

DO 90 IS • ** - 1 , VSii; 

IHOW 1-IEOW If 1 
JCOL l^JCC'L 1 1 
JCOr. 2= ^CU2+ 1 


C 

C 

c 

c 

c 

(7 

10 


POSITIVE IDErniY matrix 
MAT?. ly {IRO W1 , JCCL2) = 1 . ODD 
17 (T Y.M) 20,20, r: 

NEGATIVE IDENTITY MATFI.X 
JCCL 3 = JC0L.3+ 1 

MAT? II (I.IOW1 , JCCL 3) =-1 .ODC 


C 

C 

20 

C 


c 

c 

70 

C 

80 

90 

100 

c 

c 

c 


DEMAND FOR SENSORS 

mat:. IX (T’owi,nmat) =-d-.m (isen, ntime?) 

JCCL U*NCJL 
DO 80 ISC=1 ,;:.S.A7 
JCOL U = Jcr L^f 1 

— SENSOR F ELP.TICNSHI? 

MAT?. IX (IPO n , JCOL4) =-EAT (ISDN, ISC) 

IF (I3EK-1) 70,7C, SO 

COST EOF. .S/C 

MAT? l:-: JC0L4) *-CCST ( ISC , K TI.M ER) 

CO .NT 1 1; ” M 
CO NTINUE 
COST IN US 

RUN SIMPLEX ALGORITHM 


200 CONTINUE 

CA LL 7 V3F LS ( .M AT® 1 1 , r M A T, N N at , I RO 4 , JCO L) 
IF (I ROW. Eg. 0.0?. JCCL, IQ. 0) GOrC300 
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C 

C 

c 


300 

310 


330 

C 


CALL PIVOT3 (.1ATF.IX, 

GOTO 200 '* 


N ’.’ lAT , IMD £ Xa , 


I ? iDi ; xC , IHOW , JCOL ) 


D S~Sr..'!:VE TilE 

determine the 


N’U:i3E? CF S/C 
total LAL’HCH 


TO LAUNCH 0? EACH '’'/PE 
COSTS 


D031 ‘)I2A?=1,I?5A"' 
NUM5 A-" (12 A P) ^0 
KHfl=MHAT-l 
TCOS T = 0.0 


DO 33 OlCHKsI, MMM 
I*'D= I-IDE/. = (ICHK) 


TF (I NO. GT. VFAT) GOTOJ30 

J '- iAl ) ♦.POOCO 

-oar ^ 


I He) 


RETU r.M 


flOO no 8 01 1 = 1,5 
90 1 lilJMS A~ (I) = 1 
9 ETU ?*' 

END 

?V8FL9 :S8C- 

HSTNG *,15 

STr lU, 12, 12(13) 


i ’ iAT , M , N , i ? o <<, jccL ) 

^ S-AL DOOSLE-??'’''! c-r,»j 
. r.e.'a..,/. Ve.T^. 0.1 5 A n)C 1y75 bY R. HECKSTRCTH 


* 


* 


♦ 

* 


P EGISTM PS : 


2 

3 

a 

5 

6 

7 

8 
7 

1 0 
1 1 


PAsr; 705 (m,jm:n) (.H-i)=e 
BASF ^0^ CMI'I,:;) ( n - 1 )* m *8 


Mir ? 0 «, T.iZN BAS :. 70 F 
HINCCL, BASE ?0F ( '^t;; 

LOOP inoe:-:, 'jcol'"’ ’* 
HOW IVCPSHENT (=8) 

HMO LIMIT (=(«-2)*5, 
COLUMN increment (■'.♦8) 

VMO LT-ii: ((H*V-;i)*8 = 


(IF.OW, JMIN) 
r JCOL) 


M ( M - 2 ) * 3 ) 


'ISO* • 


FLOATING-POINT REGISTERS; 

0 MI!J 
2 MAX 

-1.D-30, THEN l.D-30 
6 VALUE 


2,4 ( 0 , 1 ) 
3,3 (0, 1) 

4.0 ( 0 , 1 ) 

2.0 ( 0 , 2 ) 


add? (M) 

A DDR (N) 

A DDR (MAT) 


0010 

0C20 

0030 

0040 

0050 

0060 

0070 

OOtiO 

0090 

0100 

0110 

0120 

0130 

0140 

0150 

OlbO 

0170 

0180 

0190 

0200 

0210 

0220 

02JO 

0240 

0250 

0260 

0270 

0280 

0290 

0300 

0310 

0320 

0330 

0340 

0350 

0360 

0370 
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..... 

L 

3, 0(0, 3) 



BCTF 

3,0 



SP 

7,7 



LR 

11,3 



IK 

10,2 



L 

fl,«r» ft» 



LE 

3,11 



S3 

11,2 



3LL 

11,3 

> 


IF. 

10,2 



SP 

5,5 



SLL 

10,3 



3CTR 

2,0 



SR 

6,6 



LP 

9,2 



SLL 

2,3 

» 


S LL 

3,3 



A?. 

2,4 



A P 

3,4 



OCTr 

9,0 



LD 

0, = D* 9. 999F50 



SIL 

9,3 


* 



1 

* 





CNC? 

0,9 


DO 10 

ID 

0,0 (7,3) 



BL 

END10 



L F 

5,7 



ld 

0,0(7, 3) 


END 10 

B !CLS 

■» , S , D 0 1 0 

1 

« 




* 





ID 

C,=D' O' 



BNL 

KAV 



L 

12,12 (0,1) 



LP 

7,5 



SF.L 

7,3 

1 


•LA 

7,1(0,?) 




7,0(0,12) 



LD 

4,=D' -1.6-30' 



L D 

0,=D' 9.999L50 



S t 

6,6 



S R 

7,7 



AS 

5,4 

¥ 


SR 

7,10 


* 




* 





I NO? 

0,9 


D020 

CD 

4,0 (6,5) 



BL 

ENO20 



LD 

6,0 (6 ,2) 

» 


DD 

6, 0 (6 , 5) 



CDS 

6,0 



B ML 

SND20 



LR 

7,6 



L DP 

0,6 


END20 

a XI.E 

6,10 , 3020 



L 

12,16 (0,1) 

I 


A R 

7, 10 



SR 

6,6 


N 

N-1 

lNDi:X 

N-1 

X (N-1) 

(N-1) ‘I 

(N-1)*1-N= (N-2) *X 
Ji* (N-2) *8 
K 

HINSOW 
X *9 
«-1 

KINCOL 

«-1 

(M-1) *3 

(N-1 ) *3 *.‘l 

PASS 

PASS 

M-2 

XIN 

(X-2) *3 


XIN VS. «A? (ININ,N) 


IF (XIS.GS.O) 

30 TO 'l.iX; LLS2 DO 2^^ LOG? 6 SilUHN. 
ADD (IF. On) 

CCN7E9T XINFOW TJ FCHTP.AN FOKM. 


mow 

KIN 
INDSX 
' JCOL' 

BASE FOR KAT (IPOW,JKIN) 
- (K*S) = FOP.TF.AM ZERO 


-1.D-3^ VS. /.AT (IROW,J«IN) 

NAT (K,JKrN) 

/:iAT (I?OW,JMIN) 


•JCOL* = JMIIl. STORED BELOW 
KIN = VALfJE 

ADDR (JCOL) 

CONVERT MCCL‘ TO FORTRAN FORK. 


0380 
0390 
0400 
0410 
0420 
0430 
0440 
0450 
0460 
0470 
0480 
04^0 
0500 
0510 
0520 
0530 
0540 
0550 
0560 
0570 
0580 
0590 
0600 
0610 
0620 
0630 
0640 
0650 
0660 
0670 
0680 
0o90 
0700 
0710 
0720 
0730 
0740 
0750 
0760 
0770 
0780 
0790 
0800 
08 10 
0820 
0630 
0840 
0650 
0860 
0870 
0880 
0690 
0900 
09 10 
0920 
0930 
0940 
0950 
0960 
0970 


» 
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DP 

6r 10 

JCOL/ («*3) 

0980 


ST 

7,0(0,12) 

JCOL 

0990 

* 



. 

1000 


L K 

14, 12, 12(13) 


1010 


BB 

14 

RETURN 

1020 

* 




103C 

* 




1040 

SAX 

L D 

2,*=D* -9.99«E5C 


1050 


S 3 

7,7 

M7IJCOL 

1060 


SR 

6,6 

INDEX 

1070 

* 




lOdC 

« 




1090 


CN'O? 

o,e 


1100 

D0 11 C 

0 D 

2,0 (6,2) 

r.AX VS. hAT (E,JXIN) 

1 1 10 


3Ii 

EMD110 


1120 


L P 

7,6 

KIKCOL 

1130 


LD 

2, 0(6, 2) 

KAX=KAT (.'1, j;i:s) 

1 140 

ESDI 10 

3XLE 

6 , 10, 001 in 


1 150 

* 




1160 

♦ 




1170 


CD 

2,=D' 0» 

MAX VS. 0 

1130 


3H 

'!0?F 

GREATER: GO OK 

1190 


S P 

0,0 

ELSE I?. OV, JCOL=0 AND RETURN 

1200 



2,12(0,1) 

ADD? (IRO^) 

1210 


T 

3,1t(0,1) 

A DDR (J CCL) 

1220 


ST 

0,0 (0,2) 

IROW=0 

1230 


ST 

0,0 (0,3) 

JC0L=0 

1240 


LM 

14,12,12(13) 


1,^50 


3R 

14 


12o0 

* 




1270 

« 




1280 

KORE 

S T 

7,hINC01 


129C 


A R 

7, 10 

CONVERT MIMCOL TO FORTLAN FORM. 

1300 


5 ?. 

6,6 


131 C 


OF. 

6, 1C 


1320 



12, 1(> (0,1) 

ADD? (JCOL) 

133C 


1 D 

0,=D' 9.999F50' 

'll!; 

1340 


L D 

4,=D' 1.E-30' 


1350 


S T 

7,0 (0,12) 

JCOL 

1360 


S R 

7,7 

INDEX 

137C 


S ? 

5,5 

' IPOW 

1380 


f 

6, «INC0L 


1390 


s s 

5,3 

-6 = FORTRAN ZKHO 

1400 


AP. 

6,4 

MIHCOL ♦ ADD? (MAT) DASS 120 LOOP 

’ 1410 

00120 

C 0 

4,0 (7,6) 

1.E-30 VS. MAT (I MIN, JCOL) 

1420 


31! 

SMD120 


143C 


L D 

6, 0(7, 3) 

MAT (IMIK,K) 

1440 


or 

6,0 C ,b) 

/ MAT (I KIN, JCOL) 

1450 


C OF 

6,0 

VALUE VS. Mi:i 

146C 


3HL 

EMD"20 


1470 


L P. 

5,7 

•IP04‘ 

1480 


LD? 

0,6 

HIN=VALUE 

1490 

ZNDl 20 

3XLE 

7,9,DC120 


1500 


S 3L 

5,3 

CONVERT Ir.OW TO FORTPAN FORM 

1510 


LA 

5, 1 (0,5) 


1520 


L 

12,12 (0,1) 

ADD? (I ROW) 

1530 


ST 

5,0(0, 12) 

I RO W 

1540 


LM 

14,12,12(13) 


1550 


BR 

14 


1560 

* 




1570 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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I 


I 






i 


MIKCOL DS 
PND 

PIVOT3 CSFCT 
USING 
ST.'* 


IP 


*,15 

1U, 12, 12 (13) 


* 

* 

* 

* 

* 

* 

* 

* 


* 


* 

* 

* 

* 

* 


* 

* 


* 


* 


* 


SUUPO'JTIN E PIVCT« (.MATPIX , 'i, :i ^INDiXa , IS DEZC,IF C W, JCOL) 
lAT'lX 13 AH BY *i HEAL I'CU 31 F.- ??. LC I3ION ARKAY. 

INnP.Xr (:*.) A.ND I'lO^XC (H) APt HALf'WOKD INThGE-l ARRAYS. 
OP.IOIVALLY W“ITTEt! IV. POETRA’I 3Y HHILI? A3R.AK 4/4/75 
p?.is::»:r vf.^gich 3Y r. iecksthct^ u/3i/75 
rsgi.st':k 5 .gufisg sxpcutiom or loops: 

1 I HD*:*; OF i?cw,J?iv 

2 'fiGi! IMKFF LOO? IHCrZMf.NT 

3 H IG !! H H 2 ?. LOOP 1 1 X IT 

4 ADDI'CSS OF HATP.IX 

5 :hde>: plgiste'! ro? i?iv,jcol 

6 IHDE/. FrilSIER 70? IPrV,J?"V 

7 ouTF.r LOOP isde:: 

e. OUTER LCO? IHCRE12ST 

I DU'E? LOO? LI HI': (C:-*AHGED DSTWEES T’.»0 LOOPS) 

10 I SHE'’ LOO? I?.'C?JKEST 

II I »';*•» LOO? Lo:v lint? 

12 THH''P LOOP I!!DS*: 

13-15 I IHKAC.H 

FLO A'’IHG F.LGISTE.’.G : 

0 r. ECI? 

5 r?.j" (.'*r.J*5£CI?) 


» 


f 


I 


I 


♦ 

12.20 (0.1) 

A DDR (IP.CW) 

L 

12.0(0.12) 

IR02 

BCT?. 

12.0 

IROW-1 

L R 

2,12 


S LL 

2.1 

(IRO --1 ) *2 

I 

3.12('>,1) 

ADD? (IHDiX?) 

LH 

0.0 (2,3) 

ITH?»ISDEXR (IP.OW) 

L 

4.24(0.1) 

• ADD? (JCCL) 

L 

4,0 (0,4) 

JCOL 

RCT:-i 

4.0 

JCOL-1 

LR 

9.4 

JCOL- 1 

SLL 

4. 1 

(JCOL-1) *2 

L 

5,16(0,1) 

A DDR (ISDEXC) 

I.U 

6,0 (4,5) 

IHDEXC (JCOL) 

S TH 

6. 0(2. 3) 

ISDliXR (IROW) sISDiiXC (JCOL) 

3TH 

0,0 (4,5) 

INOEXC (JCOL)*I?H? 

L 

2.4 (0, 1) 

A DOR (M) 

L 

2 . 0 (0 , 2) 

:1 

LP 

7.2 

n 

SLL 

1.3 

.N*8 OUTER LCO? INCPEHEST 

LR 

5.9 

JCOL-1 

A P 

4.7 

(JCOL-1) ♦'1*8 

ST 

5, JCOLDASS 


L 

3.3(0, 1) 

ADD?. (S) 

L 

3. 0(0, 3) 

S 

ST 

3.N 



1580 
1590 
IbOO 
0010 
0020 
0030 
0040 
0050 
0060 
0070 
0080 
0090 
0100 
01 10 
0120 
0130 
0140 
0150 
0160 
0170 
0160 
0190 
0200 
0210 
0220 
0230 
0240 
0250 
0260 
0270 
0280 
0290 
0300 
0310 
0320 
0330 
0340 
0350 
03b0 
037C 
0380 
0390 
0400 
0410 
0420 
0430 
0440 
0450 
0460 
0470 
0480 
0490 
0500 
0510 
0520 
0530 
0540 
0550 
0560 
0570 
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LR 

5,2 

S 



KS 

1,3 

.S*N 



SR 

5,2 




SLL 

5,3 

(fl*N-S) *3 SECOND OUTER LOOP LI KIT 



ST 

5,N?WLI.'<IT 




BCTR 

9,0 

JCOL - 2 



'I? 

e,2 

(jCOL-2) • « 



I.?. 

3,2 

fl 

f 


SLL 

9,3 

(JCOL-2) •'1*0 FIRST OUTER LOOP LifilT 



IH 

8,7 

OUTER LOC? iSCaZfttST 



SR 

7,7 

OUTER LOOP lilDEX 



LR 

11,12 

IRJM-1 



I P 

2,12 

. 



SLL 

2,3 

(IP.WO-1) *3 



SY 

2,IuCWBASt: 

(USED AT END IN PIVOT ASSI3NKENT. ) 

w 


h 

2,JCCLUASF 

IROW.JCOL INDEX 



3CT" 

11, C 

I ROW -2 



SLL 

11,3 

(IROW-2) *9 LOW iriiER LOOP LIMIT 



1 

in,=F' 3' 

INKER LOC? INCRL.XENT 



SLL 

3,3 

M*8 



SP 

3,10 

m*8-3 * (K-1)*, HloH i:;k3R loop LIM 

f 



i,0(n,l) 

A DDR (K A TP IX) 



LD 

0,*D' 1. ' 




DD 

0,0 (2,4) 

1 ./‘l.M?.IX(I?OW,JCOL) EECIP 



t 

1 ,:rcwi3A.sn 

IR0W,JPIV INDEX REGISTEF 



L 

2,»F' 6* 

INKER LOO? INCREMENT 



S F. 

6,6 

IPIV,JPIV REGISTCf. 



C 

6,0 COL BASE 

13 (JCOL-1) a 0? THEN IT'S AN EDGE 

• 

♦ 

BE 

2D1FC 



« 

crop 

0,5 



0UTE31 

T 

5,JCOLB.\SE 




SE 

12, 12 

IMN2E LOO? INDEX IlGISTER 



LD 

6,D (1 ,4) 

MP.J * .MATRIX (I?.DW,J?IV) 



CD 

6,=D' 1.L-3C' 



• B H 

KF.jnNIOr. 

SIGNITICANTLX POSITIVE 



CD 

6,»n* -1. 2-30' 




3SL 

SOSFF*=‘Cr 

NOT^IGNIFICANTLY NEGATIVE. 



OR MD® 

6,0 

MUD* EECI? 



w 

12,IFOWBASn 

IS (l?0li-l)=0? THEN EDGE 



BE 

ED3E1 


1 

* 

* 

CN'OP 

0,8 




L D 

2, 0(6, 4) 

MATRIX (I?IV,J?IV) 



LD 

4,0 (5,4) 

K ATP. IX (I? IV, JCOL) 



K D?. 

1,6 

♦MRJR 



S DP. 

2,1 

- lERM IN P.EJISTEF. U‘ 

1 


STD 

2,0 (6,4) 

INTO MATRIX (IrIV,J?IV) 


.\R 

5,2 

iucrembn: (i?iv,jcol) 


” * 

A F 

6,2 

INCREMENT (IPIV,JPIV) 



b>:le 

12,10, i:r;sR i 



2D3E1 

LA 

12,p (C, 12) 




LA 

5,S(C,5) 




LA 

6,R (0,6) 


1 

* 

STD 

6,0(1 ,4) 

MATRIX (IF.OW,JPIV) = MPJF. 


* 

LP. 

2, 10 



0530 
0590 
0600 
0610 
Ob 2 0 

063 C 

06 uO 
0650 
0660 
0670 
066 0 
0690 
0700 
0710 
0720 

07 30 
07a0 
0750 
0760 
0770 
076C 
07 5 0 
0800 
0610 
0320 
0630 

064 0 
0650 
0?60 
0670 
086 0 
069 0 
0900 
091 0 
0920 
0530 
0940 
0950 
09 60 
097 0 
0950 
0950 
1COO 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1060 
1090 
1 100 
me 
1120 
1130 
1140 
1150 
1 16C 
1170 



i 
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I 




CNOP 

lVMEa 2 LO 
LO 
I D® 
3D? 
STD 

AP 

DXIE 

b 

NOEr^ECT A* 
EP*’FCT AP 

Bale 


b,s 


2, 0(6,4) 

“#0 (5,tt) 

2,4 


2 , : (f. , 4) 

5.2 

6.2 

12,2, INNR?2 
EFFECT 

6, e 
i,« 

7, »l,0nTEP 1 


aAiaix (iPiv,jFiv) 
PAT?::< (i?iv,jcoL) 

- TBnx IN 3Ei3IST2?. 4 

:sro -ATsix 

incpsksnt (i?iv,jc:l) 

I NCKENENT (IPIV,Ji>IV) 


6 TKCREMfNTSD IK LOOP 
HASN'T been done Dy INNER LOOP 
iNCRSHEtiT (iaow,jr:v) 






I 


I 


I 


EDGEO 

A? 

6,S 


AP 

7,8 


A P 

1,6 

* 

* 

L 

9 , N EW 

• SECOND 

* 

C'lTEr 

LOOP: 


SINCE R6 HAS HISSED A WHOLE 
IMCPEKSNT OUT25 LOO? INDEX 

i.'Icp.chent (I2gw,j?iv) index 

NSW ODTEN LCO? 11X17 


?CW 


TKOP 

00TKF2 L 
SE 
LD 
CD 
3H 
C D 
BI.’L 

HISTSajp HD‘> 

w 

3E 


0,8 

5, JCCLBA3E 

12 , 12 

6,0 (1 ,4) 

6, =D' 1.E-30' 
HIS?E?.J^ 

6, =n» -1 . E-30* 
NOEFECT2 
6,0 

1 2 , 1? 0 W ? A S E 
EDGE? 


IKNL?. LOO? INDEX PEGISIER 
NSJ * NAT? IX (Ir.OW,JFlV) ‘ 

SIGNIFICANTLY POSITIVE 

not SIO NI7TCANILY 
N?J * RiCIP 

IS (I?0*W-1)30? THEi; ED3E 


C VO o 

I N N E ? .3 L D 
LD 
ND*? 

s n? 
STO 
A S 
A K 

BXLE 

ED1E2 LA 
LA 
LA 
STD 


L P 

0 VOP 

IKNZR4 LD 
LD 
N D? 

S D® 
STD 


0,3 

2,0 (b, 4) 

4, 0(5, 4) 

4,6 

2,4 

2 , 0 ( 6 , 4 ) 

5.2 

6.2 

12,10, IVNE? 3 

12 , 8 ( 0 , 12 ) 

5 , B (0 , 5) 

6 ,8 (0 ,5) 
6,0(1, 4) 


2,10 

o,e 

2,0 (6 ,4) 
‘*r0(5,4) 
4,6 
2,4 

2,0 (6,4) 


KATPIX (I?:V,J?IV) 
NATP.IX (I?IV,JCOL) 
*MRJ? 


- tern Ii; REJISTED 4 
INTO HArpix (ipiv,jr:v) 
increncnt (i?iy,jccL) 
increment (IPIV,JFIV) 


NATPIX (I30W,JPIV) a mRJS 


HATIIX (I?TV,JnV) 
KATRIX (I?IV,JCOL) 

*.VRJR 

“ term in RLGISTEF 4 
INTO MATRIX (IPIV,J?IV) 


1180 
1190 
1200 
121C 
1220 
1230 
124C 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
14 10 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1 5o0 
1570 
1580 
1590 
•1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1630 
16 90 
1700 
1710 
1720 
173 0 
1740 
1750 
1760 
1770 
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A B 

5,2 

IKCPEflEUT (I?:V,JCOL) 

1780 

AB 

6,2 

INCPEME.'IT (I?IV,J?:V) 

1790 

b:cie 

12,2,1NNF.?4 


1600 

9 

FP£CT2 

REGISTER 6 ISCRtXEHTBD IB LOOP 

1810 

NOEFECT2 AP 

6,8 

HASN'T BEEN DOSE BY ISNFH LOOP 

1820 

2FECT2 A». 

1,8 

increment (IBO-.i,J?IV) 

1830 

BXLE 

7,8,OUTFS2 


1840 

* 



1850 

* 



1860 

* LOO? TO DO 

WKAT ^OUTRAN STATFKENT 40 DID: 

1870 

* 



1680 

SB 

12, 12 

FORTY LOC? INDEX .REGISTER 

1890 

L 

S, JCCLBASH 


1900 

FORTY LO 

4, 0(5, 4) 

MATRIX (IPIV,JCOL) 

1910 

LCD?. 

4,4 

- MATRIX (I?:v,JCOL) 

1920 

.10P 

4,0 

♦ P.ECI? 

1930 

STD 

4,0 (5,4) 


19u0 

LA 

5,4 (C,5) 

INCREMENTS (IFIV,JCOL) 

1950 

byls 

12,2, FORTY 


1960 

L 

6,JC0LBASS 


1970 

A 

5,IrOyBA3E 


1980 

STD 

0,0(5, 4) 

MATRIX (IROW,JCOL) « RECI? 

1990 

LM 

14,12,12 (13) 


2000 

BP 

14 


2010 

* 



2020 

* 



2030 

TROWDAFE DS 

IF 


2040 

M DS 

IF 


2050 

NEWLIF.r" DS 

1 » 


2060 

JCOLDASE DS 

1 ? 


’ 2070 

K?!D 

• 


2080 

RFNPU7 




MAVN»ir, 




^1AxLSC= 4, '1AnLS«5, MA CK 0*2, .'5AXLr*»4, 

MA.vIE = 5, ;1AX3»10, 



1,2,2,3/ 3,3, 

KSK!J^ 2=1, 1,2,2 ,3, 3, 3, 3,3, 3, 
HSEN0 3=0, 0, 1, 1 ,?,2,2, 2,2,2, 
M3K'I0 4 = C, 0,0, 1 ,2, 2,2, 2,2,2, 
HSE*lO5=0, 0,0, 1 ,1,2,2, 2,2,2, 
1=1, 0,0,0, 
iKi::c 2 =i, 1 ,c,o, 

I KI '20 3=0, 1,0,0, 

I MI '-:o 4=0, 0,1,0 
i*;ix05=c, 0, 1,1 , 

IB«IX = 1,1 ,2,2, 

AKTBF1=20 . ,20. , 
aktb?2= i: . , 10. , 

A-"T3F3=10. ,10. , 

AF.T3 74s2 0. ,20 . , 

STAND 1= 1. , 1 ., 

STAND2=1. , 1. , .. . . . ■ 

STA!ID3=1. ,1., 

STAMD4=1, ,1., 

AFFAS 1 = 20 . ,20. , 

AKCA:J2= 10 . , 10. , 

AMF.AN 3 = 7. ,5. , 

AKSAF4=20 . ,20. , 

AMTB5=5., 7. ,7. ,5. ,S., 
STA»IDS*1. ,4 = 1. , 

AKEAHS=5* 5. , 


A 




4 . 
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XJ!AXC*15» 

AMINC»10. ,5. »3. , JO. , 

CALC*5* 100. , 

,\^AX3»15. , 30. , 

A-.i::3=lO. ,2C. , 

CAL5= ICO. , 100. , 

AXAMM-15. , AKI!.'.1=10. ,CAL3=100. , 

AXYLCI’ 12 . .20. ,22 . ,?0. ,20. , 

AiiNLCI* 10 ., 15. , 17. , 15. ,15. , 

C.MLC 1*5* 130., 

A XX 1C 2=0. , 20., 2 2. ,20. ,20. , 

A-NLC2 = 0. , 15. , 17. , 15. , 15. , 

CAL1C2=5* 100. , 

ISCaS 1= 3* 1 , 3*3 , u, 4,5, 5, 

LXCV:S2=3*2, 3*3 ,4,4,5, 5, 

0 ?acxp=c, 0 , ?F3 = . 45, .«? , 

?3-:=. 55 ,. 95,?3HS»4*‘’>. , .5, PB??=5*0. , POS = .95, . 9 5, 

^0:*'0=‘ = 2*'' . , ,PC’ = 2*0.,?L-'10 = 1. , 1. ,.75, .95, .95, 

??X?= .95, .9 2, . 95, .95, .95, 

??X?P=3*0 . , . 9, . 95, 

P?LS?: = U*3 . , .9, 

, 0r=i c. ,Dt;rc=i 
XAX2 = 1 , 

DEE'JG = 0, 

XAXL 3= 1 , 

A X7 Pr 1 = 1 0 . , 10., 

AX3A;n=4. 0,4.0, 

A XL\ K9 = 5. 0 , 5 . 0 , 5 . 7 , 5 . 0 , 3 . 0 , 

Axr:c=l5. 7.0, ■'.0,5. 0,4 0.0, 

A ilN P=15. 0 , 30. C, 

AX XL C1= 10 .C, 15.0 , 1''.0 , 15.0,15.0, 

AXXLC2=0. , 15. , 17. ,15. , 15. , 

AX::LC3 = 0. , 15. , 17 . , 15. , 15. , 

AX;<LC4=0. , 15. , 17 . , 15. , 15. , 

AX1:LC3 = 0. , 15. , 17 . , 15. , 15. , 

A*'KLC4=C. , 15. ,1*'., 15., 15., 

LVC:!S3=2, 2, 2, 3, 3, 3, 4, 4, 5, 5, 

LNC‘*S4=2, 2, 2, 3, 3, 3,4, 4,5, o, 
h=i5. , 

XA VO = 1 , 

XAZ?=1C0, 

AXT3 '•1=999 ., 9® '3. ,37A'! D1-=. 31, .01, Ax 73 5 = 5* 999., STANDS* 5*. 01, 
?P5= 4* 1. , ?3S = 5*1 . ,?3r.3 = 4*C . , 1 . ,?0S=5* 1. , P?XC=5* 1 . , ??X?=5* 1. , 
PFX? P = 3*0 . , 2*1 . , ?PL0? = 4*0 . ,1 . , 

5EVD 

TOTAL T.AL7. C” COSTS 
LAn*'CH COOTS (LSC= 

:0~AL S?ACHC?A?T C7ST5 

sPAcscpA-r ccsrs (isc= 

TOTAL LAI’S' CH & 3PACEC3APT COSTS 
LAUNCH 0 SPACEC2AF7 COSTS (LSC= 
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